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16031 IMPACT OF SOLID. MISSILES ON ONCRETE BARRIERS 
KEY WORDS Barriers; Concrete; Impact; Missiles; Perforating; Structural 


ABSTRACT: formulas are developed to predict scabbing and 


_ perforation thicknesses for reinforced concrete panels subjected to impacts by sald 
_ steel missiles. The scabbing thickness is defined as the thickness that is just enough to” 

_ prevent the peeling off of the back face of the panel opposite to the face of impact; 

_ perforation thickness is defined as the thickness that is just enough to allow a missile x 
_ to pass through the panel without exit velocity. In the development of these formulas, 

_ simplified models for scabbing and for perforation were assumed; engineering 

_ mechanics principles were applied in the formulation, and the Bayesian statistics were 
utilized to determine coefficients in these formulas. The Bayesian approach was used _ 
_ because it can handle the estimation of parameters better than the classicial statistical | a 
_ method when suitable data are limited, as is presently the case in low velocity missile 
impact on concret how hu 4 
_ REFERENCE: Chang, “Wen S., “Impact of Solid Missiles on a Barriers,” 
_ Journal of the Structural Division, ASCE, Vol. 107, No. ST2, Prec. fetes 16031, 
February, 1981, pp. 257-271 


16034 SHEAR STRENGTH OF SLABS IN STUB GIRDERS are 


KEY WORDS: Beams; Buildings; Codes; Composite structures; Concrete 
(reinforced); Design; Girders; Limit design method; Load transfer; o 
Reinforcement (structures); Shear strength; Slabs; Structural design cles 
_ ABSTRACT: The results of a test of a stub girder which failed by longitudinal shear i in 
_ the slab are presented. The measured force in the slab at failure is compared to the — 
hn theoretical longitudinal shear transfer strength: A procedure is outlined for the — ” 
| investigation of shear strength in slabs of stub girders. A numerical example is given to — 
=z REFERENCE: Buckner, C. Dale, Deville, Danny J., and McKee, Dean C., “Shear _ 
_ | Strength of Slabs in Stub Girders,” Journal of the Structural Division, ASCE, Vol. 107, 7 
OF 1 
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16039 EFFECTIVE WIND LOADS 0 ON FL FLAT ROOFS 


_ KEY WORDS: Aerodynamics; Codes; Correlation; Design; Pressure; Roofs; = L 


ry ‘ABSTRACT: An experimental investigation of wind-induced pressure loads acting on 
_ two square-plan flat-roofed model ‘Structures (SOft and 225ft high in full scale) is 
described. Measurements of mean, peak, and root mean square pressures acting on 
- points and over larger roof areas have been made in turbulent boundary layer flow 
_ simulating wind over urban and open country terrains. Area loads on the roofs have 
| been measured both by a pneumatic-averaging technique, and by using a large flush- 
| diaphragm transducer. Correlation coefficients and spectra of effective pressure loads 
{| on various regions of the roofs are also presented. A significant overestimation of peak 
| _ loads may occur for larger areas if full correlation of the worst point pressures acting 
A REFERENCE: Stathopoulos, Theodore, oi David, and Davenport, Alan G., 
_ “Effective Wind Loads on Flat Roofs,” Journal of the Structural Division, ASCE, Vol. 
107, No. ST2, Proc. Paper 16039, 1961, 200-28 
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16041 SKEWED BEAM AND SLAB BRIDGE DECKS 

KEY WORDS: Beams (supports); Bridges (concrete); Deflection; Design; 

*) Matrices (mathematics); Models; Skewness; Slabs; Stiffness; Stresses; Tests 

_ ABSTRACT: A stiffness method for the analysis of monolithic beam and slab skew 

| — is developed, the deck geometry being modeled with two-node beam and ine 


4 
Yq node skew plate finite elements. The T-beam action of the deck is effected by joining 
4 


_ the centroidal surface of the beam and plate systems at the corresponding nodes by 

_ infinitely stiff connections. This form of representation implies that, at a beam- il 

_ connection, the displacements of the two elements may be related to each other so_ 

a that at such a connection there exists only one independent node. Thus, the number of — 

| equations necessary for an analysis is only of the same order as that required for a flat 

| grillage solution. However, the necessity of prior calculation of section properties, as is 

REFERENCE: Just, David J., “Behavior of Skewed Slab Bridge Decks, 
Journal of the Structural Division, ASCE, Vol. 107, No. ST2, Proc. Paper 16041, 


February, 1981, pp. 299-316 a 


3 160 rerum Py 
KEY WORDS: Bending; Box girders; Bridges; Correlation; ao Field 
_ ABSTRACT: The induced static and dynamic response of a two-span continuous 
_ composite steel curved box-girder bridge, located in Seoul, Korea, has been examined. 
The deformations and strains were induced by a calibrated test vehicle. The 
experimental results, deflections and strains have been compared to results from three 
_ analytical techniques. The comparisons show excellent agreement between tests and 
- theory. Because of the high torsional rigidity of the bridge, basic beam theory is a 
REFERENCE: Heins, Conrad and Lee, Woo H., “Curved Box-Girder Bridge Test: 
Field Test,” Journal of the Structural D Division, ASCE, Vol. 107, No. ‘ST2, Proc. ‘Paper ; 
16083, 1, February, 1981, pp. 317-327 


16056 EDGE ST STIFFENERS FOR THIN-WALLED MEMBERS Lis es 
WORDS: Buckling; Longitudinal stability; Plates (structural members); 
- Steel plates; Stiffeners; Structural stability; Thin-wall structures hi 
ABSTRACT: An effective width approach is presented for predicting ultimate 
of thin-walled compression elements with edge stiffeners. In conjunction with 4 
procedures for predicting ultimate strengths, a requirement for adequate stiffener — 
rigidity is given. Stiffener adequacy is assessed as that stiffener rigidity for which the 
ultimate strength of an edge-stiffened element equals that of an element of similar 
dimensions and material properties but supported by a web at the stiffener location. 
Critical and post-critical behavior of the assembly are studied analytically and 
experimentally; however, the stiffener requirement is based primarily on the 
_ experimental results. Simpie formulations for assessing the performance of the = 


REFERENCE: Desmond, Thomas Pekoz, Teoman, and Winter, George, | “Edge 
Stiffeners for Thin-Walled Members,” "Journal of the Structural Division, 
107, No. ST2, Proc. Paper 16056, February, 1981, pp. 329-353 NTS ee ae 
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160s? DISTORTIONAL BUCKLING OF -BEAMS 


KEY WORDS: Beam columns; Buckling; Distortion (structural); Elastic 


distortional buckling of doubly symmetric I-section beam-columns. The method allows _ 
an analysis to be made of elastic uniform members under various conditions of loading, 
end support, and restraint. A computer program based on this finite element method is 7 
_ described, and the convergence and accuracy of the buckling solutions are 
_ demonstrated. ae astic is then used to investigate the significance of different end 
Zz conditions on the elastic buckling of I-beams under uniform moment. 
REFERENCE: Bradford, Mark A., and Trahair, Nicholas S., “Distortional Buckling of i 
4 I-Beams,” Journal of the Structural Division, — Vol. 107, No. ST2, Proc. Paper 


16055 LOAD DISTRIBUTION AND DOUBLE SKIN WALL 


~ KEY WORDS: Beams; ‘Load distribution; Skin (structural member); Stress “¢ 


H distribution; Structural engineering; Walls; Ww ind eve 


_ ABSTRACT: The stress , distribution j in a double skin wall system. is found to be very . 


- sensitive to the load distribution between the two skins. The understanding of wind Fi 
load distribution between the two skins is essential to justify the design adequacy of a 
double skin wall system. Wind engineering research is urgently needed in this area in- 

_ order to eliminate the design confusion and legal litigations arising from the dispute of 
the assumed wind load distribution used in the design. An example design analysis 
_ indicated that an 85% difference in the computed allowable load resulted from the two — 

- extreme wind load distribution assumptions (i.e., total load on the interior skin versus 
- total load on the exterior skin). Using an extreme wind load distribution assumption in 
- the design would lead to either an unsafe structure or an overconservative design. he 
_ REFERENCE: Ting, Raymond M.L., “Load Distribution and Double Skin Wall,’ - 
_ Journal of the Structural Division, ASCE, Vol. 107, No. ST2, Proc. eae 


16054 TOPOLOGY AND OPTIMALITY OF TRUSSES 


Iteration; 


results of Maxwell’ s and Michell’s theorems on on structural 


optimization are briefly reviewed. A special characteristic of a certain class of common > 
truss is then presented as an extension. The condition for trusses to have a 
_ charateristic is then derived, and its implications on topology, optimality and 'Stress- 
ratio design are examined. Simple and more complicated examples are presented. © Pre atl 
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| Finite elements; Ground shock waves (direct); Model studies; Motion; Shear 
iH ABSTRACT: The formulation is presented of an analytical model using the — 
element technique to compute the nonlinear dynamic behavior of reinforced concrete 
planar structures under earthquakes. The model is based on a lumped-mass approach, _ 
and utilizes the Newmark B-method to perform the step-by-step integration of = 
' equations of motion. This model is applied to trace the response of a three story shear- 
to the El Centro of The results computed 


Mot ior No. ST2, Proc. Paper 16061, 1981, pp. 395- 411 


_ 16088 TORSION AND ULTIMATE LOADS OF H- COLUMNS 


method; Structural ral stability; Torsion; Ultimate loads 
ABSTRACT: Ultimate moment capacities of H-columns with or without flexural end 7 
restraints have been computed by two methods—one including the effects of torsion © 
and warping, the other ignoring them. Comparisons are made for both the ultimate 
‘és moments and for deflections and failure modes. Some cases are indentified where 
torsional effects appear to have only a small influence on the moment capacity or 
failure mode. These include: the relatively short columns (weak axis slenderness ratio 
up to about 45) with no end restraints, and the longer columns (weak axis slenderness 
ratio up to about 75) with only minor axis flexural restraints. In some other cases, 
Serious errors can result if torsional effects are ignored. 
REFERENCE: Virdi, Kuldeep S., and Sen, Tapa K., “Torsion and Computed Ultimate 
Loads of H-Columns,” Journal of the Structural Division, ASCE, Vol. 107, No. ST2, 
Proc. 16058, February, 1981, pp. 413-426 ay 
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STRUCTURAL DIVISION © 


" - Theo objective of this paper is to establish formulas to predict scabbing thicknes 
and perforation thickness for concrete panels” subjected to solid steel 1 missile — 
impacts. The ‘scabbing thickness is defined as ‘the panel thickness that is just 
enough to prevent peeling off of the back face of the panel opposite to the — 
_ face of impact, and the perforation thickness as the panel thickness that is 
> just enough to allow a missile to pass through the panel without exit velocity. " 
The formulas will be derived using simplified models, principles of mechanics, ‘ 
Some early studies on missile impact problems (8,12) included the Modified 
_ Petry formula which was originally developed in 1910, the Army Corps of 
Engineers formula, the Modified National Defense Research Committee one? 
_ formulas, the Ammann and Whitney formula, and the Ballistic Research Labora- __ 
} 7 tory (BRL) formula. These formulas were primarily developed empirically for 
nm high velocity projectiles, ‘say over 500 ft/sec, and are questionable in their 
applicability for lower velocity missiles. Recently some impact tests (1,2,5,6, 10, 13) 
were conducted for lower velocity missiles, and several empirical formulas 
’ 7 (3,7,9, 10) were developed from these tests. Sliter (10) has examined some empirical © 
. formulas by comparing the calculated results with 103 recent solid missile impact B 


tests, and has assessed the accuracy of ‘these empirical formulas. 
in this Paper a are shown i in n Eqs. l and 2 respectively: 


i. Engrg. Specialist, Bechtel Power Corp., 520 South Post Oak Road, P.O. Box 2166, 
_ Note.—Discussion open until July 1, 1981. To extend the closing date one ~ol 
a written request must be filed with the Manager of Technical and Professional Publications 
_ ASCE. Manuscript was submitted for review for possible publication on March 5, 1980. 
This paper is part of the Journal of the Structural Division, Proceedings of the American — 
Society of Civil Engineers, ©ASCE, Vol. ST2, 1981. 0044- 
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IMPACT OF SOLID MISSILES ON CONCRETE HABBIERS 
g 


e= scabbing and perforation thicknesses, 
missile; v = the impact velocity; d = the missile Paeniier: rf! the concrete 
_ compressive strength; u = areference velocity; and aand 6 = numerical constants. 
5a i The coefficients Y and Z in Eqs. | and 2 are assumed to be normal aa 
__- variables. The parameters s characterizing the normal distributions will be estimated _ 
7 _ by statistical methods. The classical method of statistics to estimate the m mean Oe! 
a and the variance is to calculate the arithmetic mean of samples and the sample — g 
z variance. Accurate estimates require large amounts of samples. When test data 
are limited, as is presently the case in missile impact of concrete barriers, — a 
the Statistical estimates can be handled better by the Bayesian approach. This 
or has the advantage of being able to combine, in a rational and systematic “ 
all available pertinent information including sample data. hy 
a The Bayesian method treats the unknown parameters, such as mean and 
a Berets oy as random variables. The joint probability density function of these 
_ random variables, called a prior joint distribution, is to be constructed from 
all of the information that happens to be available, such as sample data or 
engineering jodgment, The prior joint distribution represents our prior uncertainty 
in probabilistic terms in order to make an inferential statement or to “make ee 
a decision. The prior joint distribution is combined with subsequently obtained © 
sample data to produce a revised joint distribution for these random variables, 
called a joint posterior distribution. The estimation of parameters is made from — 
the joint posterior distribution. When additional data become available, the joint — 
_- posterior distribution is treated as a prior joint distribution and a revised jolt i 
= “ posterior distribution can be obtained by repeating the same procedure. Thus 
the values of Y and Z in Eqs. | and 2 for concrete scabbing and perforation A 
wil be updated as new data is acquired until sufficient refinement is niche 


"Bayesian Anatysis of SAM FROM Distrisution (4) 
‘The coefficients Y and Z in Eqs. | and 2 are assumed as normally entins’ 
random variables. The mean and variance of each random variable will be — 
ie determined. Because of the limited amount of suitable test data available, the — B 
4 ‘Bayesian approach, rather than the classical statistics method, will be used. ai 
To begin this approach, the Normal-Gamma distribution is assigned as the prior 
g joint distribution to the mean and the inverse of the variance of the normal 
random variable. Four parameters, a,, Bo, and 
distribution, ane by the following 


a 


te CONCRETE BARRIERS 


in which ) represents the expectation; Var ( the variance; the 
value; and the variance the random variable. The values on 
the left of these four equations are estimated on the b basi of sample data, 


+ 


in which ae the of ‘samples; x, = the sample value; x, = the: 

_ sample mean. The revised mean and variance of the normal random variable ~ a 

can be computed by Eqs. 3 through 6 by using new parameters with subscript = : 

_ 1 obtained from Egs. 7 through 10. The same revision procedure can be repeated » 

_ whenever new data are available until the results have been refined ey. 


en a missile strikes a ‘the barrier is deflected that 4 
the maximum negative moment occurs at the periphery of a circle of radius 
OR where te oo bending stress is created i in the rear face of the barrier 


missile. “When the kinetic energy absorbed by ‘the bending i is large enough and 

exceeds the capacity of the strain energy along the periphery of the circle, 
the barrier starts cracking and scabbing. The scabbing thickness is calculated 
by equating the kinetic energy and the strain bie © capacity. 7 The capacity - 


in which SE = = the capacity; R= “the of the affected 
~ circular area; M, _ = the ultimate moment capacity per unit of length; and 7, 
= the rotational capacity of the barrier which is taken as the minimum angle _~ 
_ corresponding to the ultimate moment in the idealized moment- curvature relation- 
ship. Assume R is proportional to’ the product of V ds and (v y / u)* ‘in v which 
= the diameter of the cylindrical missile; s = the A gay gee. thickness; v= 
“the: missile velocity; u= the and “a” =a numerical constant 
m 


in which C, = arandom variable. 1 The ultimate moment can expressed 


. ng those four parameters of the joint distribution as follows: __ — 
| 
| 
| a 


in which f/ = the compressive strength of the concrete barrier; and -- 


a random variable. The barrier rotational capacity r, is relatively a constant. 
_ The portion of the kinetic energy of the missile used in , bending process 


in which C, = a random variable less t than 0. 5; end 1 m= the mass of the | 


PS. 


Face 


= 


oe! 1 —Missile on Concrete Barrier 

‘The - scabbing thickness of a barrier is pone by equating SE of Eq. 


u and KE of Eq. 14, and 12 and 13 into Eq. 11: 


its distribution will be we the Bayesian ‘method. 
reference velocity = was chosen as 200 ft/sec or 60.96 m/s for convenience, 


66,99 


and the constant “‘a’’ was chosen equal to 0.325 by comparing the calculated 
values from from test data. Therefore Eq. AS — 


— | d/2| | |d/2 SC 
| 
- andom variable Y is assumed to be a > a 


: ‘tee e.g., variations in steel reinforcement occur from a barrier to a barrier, 
1 the actual concrete strength varies from the nominal strength specified, the 
Bi of missile impact can occur on or between reinforcing bars, and missile — 
_ orientation varies frontal face on or slightly off normal, although both « cases 
sidered as normal impacts. These and other uncontrollable parameters — 
will af fect the impact results. The chance that all worst conditions or all favorable _ 
- conditions for the impact take place at the same time is small and this observation - q 
leads to the assumption of normality for 
The two parameters (mean and variance o”) used to characterize the 
distribution | are considered as random variables in Bayesian : Statistics. The F pl 
distribution for and is “assumed as Normal-Gamma distribution. The 
prior information for the joint distribution is estimated as E(p) = 1.8, Var : 
(@) = 0.01, E(o*) = 0.04, and Var (o”) = 1.44 x 10~*. The corresponding — 
parameters of the prior distribution are Seg by Eqs. 3 through 6 which > 
= 13.11, By = 0.48, Ho = 1.80, and 
posterior joint ‘distribution will be by combining the prior 
distribution and sample data. The test data shown in Table 1 were selected © 
as sample data from recent impact tests. The criteria of the selection are that ; 
the test result should yield minor scabbing, i.e., the secondary concrete particles 
, emerging | from the rear ir face of the panel should be small, have very | low velocities, - 


solid steel with a cylindrical impact ar area. a. For ‘these ‘selected am 
X, = 1.865, and (x, — X,)’ = 0.1505. The parameters of the 
joint has been computed by Egs. 7 through 10 to yielda, = 17.11, 
 -B, = 0.56, p, = 1.84, and A, = 12. By using Eqs. 3 and 5, the two parameters al 
of the random 1 variable are computed: E(p) = 84; and E(o’ )= = 0. 0348. 
_ Therefore, the mean scabbing thickness of a concrete barrier am 8 to a 


cylindrical solid steel can be calculated by 


= 200 or 6, 096 the units” used in the 


The safety margin factor corresponding to a confidence level can be assigned - 

r to 5, of Eq. 17 by utilizing the mean and the variance obtained herein. For 
instance, if the safety margin factor was chosen to be equal to 0.2 or y equals — 
the confidence level is 97.6% which is evaluated by = 2.21|n = 


When a missile strikes a concrete barrier, the barrier develops” resistance 
to the missile’s penetration. Assume the resistance function is maximum when | 
the missile just strikes the barrier and decreases ih to zero when — 
e the missile meus the barrier thickness as expressed in Eq. 1 


= 
| 
| 
q 
(18) 
3 
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which F(x) function; f = the resistance when x= x 
the distance along the concrete thickness measured from the impact face; 


Note: W = missile weight, v = missile velocity, d = missile anata. f! = concrete ‘ 
compressive strength, = concrete panel thickness. C = (200/v)°'* (mv’)° */ 


light scabbing © 
‘moderate 


moderate scab- 


{ 


light scabbing 
2.033 | very small wa 


chips of con 
crete dis- 
lodged from 


periphery 


‘Tear face 
fractured 
1.2 cu in. parti- 
traveled 1 


2 cu in. pert 


cle traveled 1 
ftatlfps 
0.6 cu in. parti- 


 cletraveled 
«ftatl 


y =t/C, 1 lb = 453.6 gm, 1 in. = 2.54 cm, 3,300 psi = 22,700 kN/m?. — vale ce’ 


 ande = the p perforation thickness of the barrier. The energy required ‘for Poa 
_ missile to just perforate the barrier is evaluated by Eq. 19: 


The re resistance f is further assumed to be: 


n which a the diameter of the cylinder of length e to resist the missile’s 
penetration and assumed to be proportional to the product of the missile diameter 


(20) 


ds ‘second ‘ches | inch | ches| ches 
= 213 | 1 | 5,910 | 6 | 3.476 |1.726 | 
‘se | | 1 | 6 | hum | 
| 4810 | 6 | 2.985 
| | 18 | 8.853 
== | 
4,300 | 6 | 3.429 |1.750 
| 
| 


barrier; and c = a proportional constant. 
Match 


The resistance f can be rewritten as in Eq. 21: 


‘in which =a random variable; and b = a numerical constant. 


ye Let the kinetic energy of the missile consumed in the perforation process 


‘in Z= (GBC, Z is assumed to be a normal 
random variable following the same reasoning as for scabbing. The determination 
of Z will follow the Bayesian approach by assuming the Normal-Gamma joint 
distribution to the mean and the inverse. of the variance of Z. The numerical — 


i was chosen to be equal to 0.5. Thus, mc 23 is aes as in 


Test where no ‘scabbing 
@ Test where scabbing ocowred 


“AIG. 2.—Test Barrier Thickness /, in inches versus sus Scabbing Thickness s, in inches - 
Colouleted SA Formula for Seabbing tt 


4 
| 
| 
: 
20, 
4 
| 


0.25 | 
AY 


a the test data in Ref. 10 and judgment is estimated as follow : E(w) : = 

Var (uw) = 0.0025; E(o”?) = 0.01; and Var (o’ 10~°. The 
‘parameters for the prior joint distribution are computed by Eqs. 3 through ' 
6 which yield: a» = 12;B, = 0.11; 


y When additional sample data become available, the parameters of the posterior 9 _ 


_ joint distribution can be obtained, and thus the mean and the variance of the 


beans 3 z. At this time it is ; recommended that the mean perforation hd 


Test where scabbing occurred / 


3 tests 
‘Sapo 


FIG. 3.—Test Barrier Thickness /, in inches versus iit Thickness s, in inches 


Calculated by the Modified NDRC Formula for Scabbing 


(e,) of a concrete barrier to a cylindrical solid steel missile 


in which u = 200 fps or 6,096 cm/s an upon the units used in the 
equation. Note that a similar empirical | formula for perforation thickness i 


4 


Comparison OF Mooirieo NDRC Formutas- ULAS 


__ The modified NDRC formulas (see Appendix II for the complete set of formulas) = 


are Perhaps: the most widely known in the field of missile aed ener 


a) | 
j 
q 
: 
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= CONCRETE BARRIERS 
It can be that predictions by these for: scabbing thickness 
ie are in good agreement with the test data available at the present time. However, 
predictions by the NDRC formulas for perforation thicknesses are not in equally — 


close agreement with the test data but provide conservative estimates. Moreover, 
# further examination of these formulas reveals the probable existence of a hidden 


i’. po modified NDRC and semianalytical (SA) formulas are studied and a 
se Fig. 2 shows the plot of all the pertinent test data for scabbing described i 
4 in Ref. 10. The horizontal axis of the plot represents the calculated — 


Test where no perforation occurred 


“FIG. 4 —Test Barrier Thickness ¢, in inches versus Perforation Thickness, vated in aaaedll 
thickness using , the SA formula, and the vertical axis represents the thickness 
of the tested concrete barrier. It should be noted that each test data which © 
includes five parameters is able to be plotted on the two-dimensional figure 
as a single point and that due to coincidence a point on the plot may represent 
_ more than one test data. Ideally, all points falling below the 45° line on the 
plot. ‘should be tes test samples s where scabbing occurred and points above the 45° 
line should be test samples where no scabbing occurred. It can be seen from 
Fig. 2 that the predictions of the SAF for scabbing are in very good agreement 
- with the test data. A similar plot using the Modified NDRC formulas for scabbing 7 
7 is presented in Fig. 3. It can be seen that these formulas en almost oquaty 
in dividing the scabbing and test data. 


| 
a 


a Test where perforation 


FIG. 5. 5.—Test ‘Barrier Thickness t, in inches versus Perforation ‘Thickness e, in ine 


Calculated by the Modified NDRC Formula for Perforation 


SAF 


v = 295 fps 
fo= psi 


aid 


J FIG. 6. &—Seebbing Thickness s, in inches and Perforation Thickness e, in we 
Calculated by the SA Formulas versus Missile Area a, 


«| 


10 using the SA formula for perforation and the Modified NDRC 
for perforation respectively. From these plots it can be seen that the SAF 


_ Fig. 6 shows a plot of calculated scabbing and perforation thicknesses against 

pe os impact area for both the NDRC and the SA formulas. In this plot 2 a 

turbine missile having a weight of 500 Ibs with a velocity of 295 fps has been 
assumed. The concrete barrier has been assumed to have a strength of 5, 500 
psi. Because the orientation of the missile can vary as it strikes the barrier, 

_ the impact area will vary. This area has been assumed to range from 40 * 
a  in.-150 sq in. It is seen that the NDRC formulas predict a greater scabbing “4 
thickness as the impact area increases. This is contrary to’ the rational expectation — 

that greater impact areas, all other parameters remaining constant, should lead 

we. decreased barrier thickness requirements. To a lesser extent, a similar 3 

a is exhibited in the plot for the NDRC perforation formulas. This = 
contradictory prediction by the Modified NDRC formulas is the deficiency 

referred to herein. Note, however, predictions by the Modified NDRC formulas 4 
do: not t invariably exhibit this deficiency but in some cases only. The SA formulas b. 
are shown to behave in conformity with the rational expectation. __ > a <a 

This study has developed two semianalytical formulas, Eqs. 17 and 25 

respectively, for concrete scabbing and concrete perforation for concrete barriers 

subjected to impact by y cylindrical solid steel missiles. Unlike previously developed 

- empirical formulas, these two formulas are unit consistent, since basic principles 
_ of engineering mechanics have been applied and consistent units were maintained — 
in each equation in the formula derivation, 
_ The undetermined random variables Y - Z in these two formeles were 

; evaluated by Bayesian statistics. One of the advantages of the Bayesian approach _ 

is that all other available information in addition to sample data can be utilized Bs 
‘d in the determination of Y and Z. This is particularly useful at the present time 7 

because suitable data for this determination are limited. The Bayesian approach — ; 

also allows the revision of Y and Z by Eqs. 7 through 10 systematically as 

new data become and formulas can be refined progres- 


_ the Bayesian estimates of mean and variance of Y and Z, and using the table 
of the standard normal distribution | to obtain the desired results. 
These fe ormulas are applicable to current solid steel tornado and turbine — 
problems in nuclear power plant design. The values of Y and Z for each ie 
_ datum listed in Refs. 1, 5, and 10 were calculated based on these formulas ps 
g Be and barrier thickness, and they appear to predict with reasonable accuracy the 
s of barrier damage. These test data cover missile velocities ranging from 
55 fps (16.73 m/s)-1, 023 fps (311.81 m/s), missile weights from 0.24 Ibs (0.11 
__kg)-756 Ib (343.64 kg), missile diameters from 0.79 in. (2.01 cm)-12 in. (30.48 


Fe The developed formulas also provide a rational way to determine a ay 
roin f: according to elected idence level his isd hy uiating 
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batrier from 2 in. (5.08 cm)-24 in. (60.96 cm). 
4 hers a final note, since test concrete barriers used and referenced for the 
development of these formulas are normally reinforced concrete (0.3%-1. 5% 
each way), it may be conservative when these formulas are to be applied a 
heavily reinforced barriers, especially for perforation. 
Appenoix |.—Bavesian ANALysis or SampLes FROM NorMaL DistasuTioN ae 
Suppose ‘that the variables X,, form a random sample from a 
distribution for which both the mean and the variance are unknown. 
the Normal-Gamma distribution is assigned as the joint distribution to the 
_ parameters y and T, where T is called the precision and defined as T= ‘Ae, 
_ it was proved that the posterior distribution that is derived from the observed 
values in the sample is again a Normal-Gamma distribution. The subsequent — 
theorem describes the conjugate property of the prior distribution. The theorem | 
¥ is stated and the statistical | relationship between the variance o *and the precision 
Theorem (4). —Suppose that , ..» X, form a random sample from a normal " 
‘distribution for which both the mean p and the precision T are unknown (—o 
s 4. < o and T > 0). Suppose also that the joint prior distribution of » and hr 
_ T is as follows. The conditional distribution of » for any given value of T 
is a normal distribution with mean and precision Ay < pp, < and 
‘ee _> 0), and the marginal distribution of T is a gamma distribution with parameters © 


a, and B, (a, > 0 and > Then the joint posterior distribution of » 


and T, given that X,= x,(i = 1, ..., n), is as follows. The conditional distribution 
for any given value of Tisa distribution with mean , and precision 


in which z, = = : the s; sample mean. Also, the ma marginal distribution of T is a ‘gamma . 


= Bo +— > 


From the marginal distribution of the mean it is found: 


= practical problems, the Prior i information was given ‘in terms of the variance 


| 


4 
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rather than the precision T. The f between 
and T is derived i in the following: 
kay 


fi 


1 


l — 2}. 


in which x = 
FS equal to 0.72 for flat nosed bodies, 0. '84 for blunt and bodies, 1.00 for average . 3 a 
. bullet nose (spherical end), and 1.14 for very sharp nosed bodies; W = the 
: he missile weight, in pounds; d = the ‘missile diameter, in inches; and v = ed 
a missile velocity, in feet per second. The factor K isa concrete 


factor which was pro to be: 


a 
Bs 
— nine the concrete ed 
The formulas 


K= 

which fi= 

‘The scabbing and thicknesses are predicted by 4B and 44 


i ‘tid were developed by regression analysis in the specified x/dranges: 


132+1.24—; for 135s —<135.. 


be in 1 which s and e = scabbing and perforation thicknesses, in inches, respectively. 
_ For the small values of s of x/d the ome Sere curve- fitting extrapolation equations 


were developed: 


with the result predicted by formulas developed i in this p paper (SAF). 
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StRENGTH OF SLANE Srp 
The following symbols are used in this pi paper: = 


numerical constant, missile impact area; 


e= proportional constant; 


«=D = diameter of cylinder of concrete resisting ileal s penetration 


missile diameter; 


energy required to O just perforate concrete panel; 


concrete compressive strength; 


= energy 0 of missile consumed i in scabbing process; 


M, ultimate moment capacity per unit of length of concrete panel; 

the n = number of sample data; Toa 


= radius of affected circular area upon missile strike (see Fig. 1); 
concrete rotational capacity; 
strain energy along the periphery of circle of radius 


oncrete scabbing thickness; 


thickness of tested concrete panel; 


reference ‘velocity; 
variance; 


missile impact 


= distance along thickness measured | impact act face; : 
y = random variate of Y; 
Z = random variable for e; bat he 
= parameters of prior joint distribution; we 
= parameters of posterior joint distribution; 
= mean of random variable; 


variance of random va 
Vi 


tree 
| 
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of Steel Construction (AISC) Specification 
(8) for the design of shear studs in composite steel-and-concrete beams does 
not consider the longitudinal shear strength of the concrete slab. The specification - 
establishes the allowable horizontal shear load per connector ba ed upon the 
type of connector, yr, the strength an and unit weight of concrete, and the aracteristics 
of the decking used to form the slab. Limitations for the transverse and ee 
spacing of the connectors are also specified. ake 

a Longitudinal shear failure in a composite beam is charactorined by splitting — 


of the concrete slab parallel to the longitudinal axis of the member. Shear 


- stresses develop along the cross section to to transfer the compressive fore e from z 
“the: slab into the shear connectors as shown in Fig. 1. 
_ The importance of transverse reinforcing steel on the longitudinal shear capacity © 
of the slab in a composite beam has been covered by Johnson (4). The British 
specification composite design (2) requires a minimum area of trans-— 


- Recent trends in composite construction increase, in the writers’ opinion, 
the possibility of longitudinal shear failure in certain composite girders. These 
. Use of “stub girders” 
2. Use of lightweight c concrete which generally result in thinner -— A 
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SHEAR STRENGTH OF SLABS IN STUB GIRDERS 
5 
| 
7 
| 
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3. Use of compos | metal decks which, “since the are usually 
parallel to the girders, reduce the concrete shear — area. ee = 


test of a stub girder which was at 
_ University (3) indicates that the longitudinal shear strength of the slab may 
bet the weakest link in the shear transfer mechanism for such members. The 
purposes of this paper are (1) To provide a brief description of 


4 iA compnession IN SLAB 


 THEADED sTuDS @ 5%" (TYPICAL) m 


test; (2) to to caution designers of stub girders. of possibility of a 
shear failure; and (3) to call attention to the use of transverse slab reinforcement z 


simply sapported composite st stub 26 ft (. 92 2m) was loaded 


TuUctea Of a lOWRT CHOT Wi A WC st 5, 


of Testing and Materials (ASTM) A36 steel, and 2 a solid 48-i 
=o by 3.5-in. (122-cm by 8.9-cm) sand-lightweight concrete slab with unit weight — 
of 110 pef (17.25 kN /m’ ) and 3,500-psi (24,150-kN /m’) compressive stren ath. 
The slab contained transverse reinforcement consisting of #3 rebars with = 
— 60,000-psi (413,400- kN /m* ) yield strength spaced | 12 in. (30.5 cm) on centers” 
in the middle of the : Slab. The stub | elements | contained 36-5/8-in. (1.6 cm) | 


14-3/4" A325-xX @ 4"0.C. EACH SIDE 


(wom FIG. 3.—Section A-A of Test Girder 


in concrete at failure was 298 kips (1, 326 kN), hadi is below the “atimate 
_ value, C, based on the smaller of the compressive strength of the slab oped 


A. = 500 kips (2,225 kN) whet 


| 
——— 
[| 
i 
| 
ares ‘FIG. 4.—Slab Failure Region in Test 
—— diam by 2.5-in. (6.4-cm) long headed studs spaced within AISC limitations. - a 
Details of the test specimen are shown in Fig. 3. » 
The failure mode for this girder was by longitudinal shear and occurred at 7 
4 


ort the ‘strength of the shear connectors 


computed in accordance with the AISC 


_ The failure region of the test girder is shown in Fig. 4. iE ingore nae 


Stus Girvers in Compcere Fioon System ow 
The behavior of a a stub girder in in an cont: floor system n may differ fer significantly — : 


from that of the test girder. In a a complete floor system transverse membrane - ' 


forces can develop in the slab which increase its longitudinal shear strength. a 


On the other hand, the slab over the web voids may be subjected to additional 


bending moments caused by the floor loading which could reduce this strength. — i 


The flexural strength of a conventional composite prea is hon 30% greater 
than the flexural strength of its noncomposite steel section. In a stub girder 
the flexural strength is typically more than 100% greater than for the steel 
section comprising its bottom chord. Thus, loss of composite action is potentially — 
more catastrophic | in a stub girder than in a _ conventional ‘composite girder. 

From a qualitative standpoint girders which ‘would be most. susceptible 
longitudinal shear failure would be in terminal spans (where the slab is fre ~ 
along one transverse edge), 
_ Since the floor slab generally spans parallel to the longitudinal axis of a 
girder, negative flexural reinforcement is not required in the slab over the girder. 

A nominal amount of transverse reinforcement may be provided over the girders — 
to serve as crack control and temperature and shrinkage steel, but this steel — 


more than likely will not be to develop the 

sof the slab. 


An analysis taking account of these factors has not — be at ie 


Tentative Desicn 
In the absence of the e writers: | that 


The longitudinal shear strength should be sufficient to develop the capacity 


‘connectors as indicated in Fig. 5. The ultimate longitudinal shear force, V,, 


2. The critical planes for are to rows 

> 
to be transferred across these Planes can be approximated by the — _ ae 


Q, = the ultimate horizontal shear capacity of the connectors; 


| 
4 | 


can be assumed to ‘be uniform (7). ‘Thus, for the case of a slab 
sides of stub nominal shear stress can be 


in which = the length o of the stub under consideration. 


ield st of the reinforcement. The limits for sand-lightweight concrete are ay ae 


= [0.8p/, + 200 
Wig: 
4. At least half the reinforcement required by ee 5-7 should be 


in the bottom of the slab (4). 


‘fox FIG. Critica Longitudinal Shear 
5. The longitudinal shear stress should be te” vary ‘linearly 
its critical value, given by Eq. 4, to zero at the extreme edges of the effective 
_ width of slab. The area of transverse reinforcement can be reduced accordingly. 
Lonerrupinat Shear Stress in ‘Test , ang 
a Consideration of horizontal equilibrium that the resultant c: compressive 
- force in the slab of a stub girder must equal the tensile force in the lower 
chord at sections through the web openings. These forces were determined 
fer the test girder from strain gage readings taken on the lower chord, and 


is 
4 
: 3. The shear stress computed by Eq. 4 should not exceed the limits propose ) 
— os by Mattock, et al. (5,6). For normal weight concrete these limits are a 

a q 

a 
| a 
i | = 

a 
| 
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The critical longitudinal shear planes for ‘the test are by | 
distance of 5.6 in. (14.2 cm). Thus, the critical longitudinal shear force in the 


test girder can be approximated as 


which a pry over ‘to 60-in. ase. stub 
= 570 psi (3,930 kN/m*) 
te the ‘slab of the test girder, the criteria of expression 6 yields a limiting 
value for v, of 375 psi (2,580 kN/m’). This criteria was developed for the — 
case of an initial crack across the shear transfer plane, such as might be caused 
by shrinkage or slab flexure in an actual floor system. The longitudinal shear 
planes of the test girder were initially uncracked, and were capable of developing ~ 
higher shear strengths than predicted by this criteria. From the data obtained © 
by Mattock (6), the shear t transfer strength of an initially uncracked shear plane, 
in sand-lightweight concrete, is approx 200 psi (1,380 kN/m’) higher than for a 
an initially cracked plane. Thus, there is good agreement between the nominal 
stress computed for the test girder and the shear transfer strength obtained 


To illustrate the design criteria the required stub length and area of transverse _ 
reinforcement will be computed for a stub girder with the following characteristics: a 
36 5 /8-i -in. (1.6-cm) diam headed studs arranged in three rows at ‘minimum 
transverse spacing, ¢ = 5.25 in. (13.34 cm), f’ = 4,000 psi (27,560 kN/m’), | 
longitudinal slab atateoniene consisting of seven #6 bars arranged such that - 
one bar is contained within the critical shear planes, f, = 60,000 psi (413,400 
_kN/m?), and sand- -lightweight concrete with a unit weight of 110 pcef (17.25. 

_kN/m? ). The AISC Specification will gov ern connector strength and spacing. - 
_ The ultimate shear capacity of a 5/8-i -in. (1.59-cm) diam stud is taken as 
twice the working load value tabulated in Table 1.11.4 of the AISC Specification. ‘ ; 
This value is reduced by the lightweight concrete coefficient, 0.83, in accordance © 
with Table 1.11.4A of the specification. The minimum transverse spacing of 7 
the connectors is four times the stud diameter, therefore, the critical longitudinal | ~ 
shear planes are 5.6 in. (14.2 cm) apart. Thus, sali longitudinal shear — 


= (36 x 0.83 x 2 x 9.2) - (0.85 x 4.0 x 5.6 x 5.25) - 4 


800 psi (5,520 kN/m’). Thus, Eq. 4 can be solved to obtain the minimum 


> criteria of | expression 6 limits the ultimate longitudinal shear stress —_ 


1, = 423,000 + (2 x 5.25 x 800) 


(11) 
which is is adequate to accomodate the shear cc connectors without violating the 


| 
| 
= 


_ AIBC minimum spa spacing six times mes the stud 
ss The required percentage of transverse reinforcement corresponding to the 
50-in. (127-cm) minimum stub length is obtained from. a first criterion of | 


= (8 800 250) + 8 x 60,000) = 0.0114. 
“which is ante idl six times greater than would normally be furnished . 
temperature and shrinkage requirements in a slab of Sin 
Sun Summary AND Conctusions ‘ 
One of the pa in the design of a stub girder is to create web openings * 
which are large enough to accomodate mechanical ducts, electrical conduits, ‘_) 
= intersecting floor beams. In order to make the web openings as large as 7 
> arrange the 1 required number of shear a 
connectors s at the minimum spacings permitted by the design specification. Such — a 
an arrangement can result ina longitudinal shear weakness in the slab, as evidenced 
— longitudinal shear capacity of the slab of astub seem should be investigated ; 


in "part, by a Graduate 
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m ._ in this paper should provide an adequate investigation of this failure mode. = 
_ The testing of the stub girder was performed in the Structural Laboratory §& 
i 


following symbols are used in this pape 

A; = area of concrete flange in composite member; Bo 
Aa ‘ad area of longitudinal slab reinforcement between longitudinal shea S.) 
resultant compressive force in concrete flange; 

specified compression strength of ave 


Se 


= length of segment containing ‘shear « connecto 
= ultimate horizontal shear strength of group of connectors; jectors; Sh 
thickness of concrete in flange of composite member; ‘Sat of 
= ultimate longitudinal shear force in flange of composite member; sr 


sq 


Mie 


P 
q 
1 
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EFFECTIVE Winp Loaps ON FLat Roors 


By Theodore Stathopoulos," M. ASCE, David Surry, 
and and Alan G. G. Davenport,” M. ASCE 


At any location on ‘i exterior + of a structure, the wind- induced wan 


_ (both positive and negative) are likely to be highly unsteady, and also to — 7 
_s significantly from point to point. This is due to turbulence in the flow, and 
turbulence caused by flow separation from the sharp edges of the building. 
The scale of the resulting pressure fluctuations must then depend on both can 
building size and the size of the eddies in the oncoming wind. Since these 
_ pressure variations are not necessarily well-correlated either in time or in space, 
peak loads affecting large areas would be expected to be less than the 
sum of the peak loads acting at individual points within such areas. Point loads 
_ are significant for the design of fasteners; however many other structural elements — 
have tributary areas of significant size and thus their design should recognize 
this potential alleviation of load rather than conservatively assume full correlation — 
 betweenthe peakloads§ 
Numerous wind tunnel experiments on models have been carried out in a 
Bs ‘past to determine appropriate design pressures, but until recently almost all 
of them have been performed under uniform steady- cofeng 
: ——- often yield results considerably different from those observed in 
full scale. Modern boundary layer flow simulations have led to results in 
% satisfactory agreement with full-scale measurements (see, e.g., Refs. 1 and 2); 
however even among these recent studies, many have been mainly confined 4 


Asst. Prof., Centre for Building Studies, Concordia Univ., Montreal, Quebec, Canada * 
a ? Assoc. Research Dir., Boundary Layer wae Tunnel Lab., The Univ. of Western 


*Dir., Boundary Layer Wind Tunnel Lab., the Univ. of Western il London, ~~ 
 Note.—Discussion open until July 1, 1981. To extend the closing date one ‘month, 
‘a written request must be filed with the Manager of Technical and Professional Publications, _ 
_ ASCE. Manuscript was submitted for review for possible publication on June 22, 1979. i 
a This paper is part of the Journal of the Structural Division, Proceedings of the American — 
of Civil Engineers, -@ASCE, , Vol. 107, No. ST2, February, 1981. 0044. 
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Roofs i in are vulnerable to wind action, especially local areas’ 
near windward edges and corners. Near these locations, large local suctions — 
- (Caegative pressures) may occur which tend to lift roofing elements. The dynamic 
>: characteristics of these wind forces and their effects on roofs have not received - “ 
much attention, perhaps due to difficulties i in the 7 
an This paper describes wind tunnel of p pressures on ‘points and 
on larger areas of flat-roofed buildings, and compares some of the results SoG * 
current practice. Special emphasis has been given to the dynamic characteristics _ 
of the loading, as well as to the combined wind forces acting together on different _ ; 
roof regions. A novel application of pneumatically averaging the unsteady "4 
_- pressures acting on ‘Specific: model areas has been used to define peak area ap 
x loads for different roof regions. Parallel measurements using large a 
i transducers set flush with the flat roof have verified this technique. Conventional — 
point pressure measurements have also been made. ‘The comparison of _ 


+ 


4h 


FIG. 1 —Ghetsh of Models: (a) With a Dense Grid of Pressure Holes; ; and (6) With 
‘Openings for Flush- Diaphragm Transducer 


aii coefficients and spectra for different roof areas s help to give a ‘picture 
of the development of the flow over the building, and therefore to determine 
_some further characteristics of the random pressure loads. This paper is based 
on work reported in more detail in Ref. 10. Me 
Ug All experiments were carried out at ‘the: University of Western ‘Ontario’ . 
te == Layer Wind Tunnel Laboratory yi flap The wind tunnel has a 


j 


ROOFS 


roof height ; averaging 7 ft ( (2 m). Details of the wind are given 


by the third writer and Isyumov (4). 
_ Two plexiglas models were used as sketched in Fig. 1. ‘The first model has 
a regular d distribution of pressure t taps on the top for conventional measurements | 
using pressure transducers with | or without pneumatic averaging. The second - 
_ model has a large-diameter, flush-diaphragm pressure transducer as part of the 
_Toof surface. The diaphragm area approximates 1/16 of the roof surface 
\ corresponding to the tributary area of nine pressure holes on the first model. 
_ The general arrangement of the measuring points was made so that, for wind a 
_ perpendicular or at an angle of 45° to an edge, all roof areas could be investigated - > 
by using symmetry of flow and models. When using the flush-diaphragm a 
transducer, only the inner roof region could be investigated. Each model is 


made of two parts. The top Part includes the roof surface and can be ood 


“together with the top part, Tepresents a t taller building. 


- 

6 

2 


FIG. Profiles of Mean intensity for Two 10 Exposures 


_ The models were mounted on the wind tunnel turntable which is located 
- 80 ft (24 m) downstream from the belimouth and 64 ft (19.5 m) away from 
‘. the beginning of a random homogeneous roughness distribution. Two types of = 7 
terrain were simulated. Open country terrain (the “amooth” case) was 
a carpet. The urban environment (the ‘‘built-up” case) was Tepresented 
by a surface of approx 2-in. (50-mm) cubes distributed randomly. Typical nea / 7 
profiles are shown in Fig. 2. The free stream or gradient | velocity, V, | used 
in the wind tunnel was approx 45 fps(13.7m/s). = 
cm Point-pressure instrumentation includes four | /2-in. (13-mm) diaphragm strain- _ 
gage, differential-pressure transducers (Statham model PM 131TC) set in scani- 4 
valves under the turntable of the wind tunnel. Tubes [24 in. (610 mm) long 
and 1/16 in. (1.6 mm) in inside diameter] lead from the surface taps to the og 
scanivalves. These tubes include damping constructions to keep the frequency — 
response flat. The resulting system responds to pressure fluctuations on the 
model of up to about 100 Hz with negligible attenuation or distortion. Higher 
suffer increasing attenuation, although some response is obtained 


j 
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"parallel weleg several scanivalves, each of which i is is sampled a at a rate > of about 

used 
uy to: measure sare spatially-aver -averaged pressures acting on roof areas, by interconnecting. ah 
_ pneumatically an array of pressure taps through special multi-input manifolds. ey 
= detailed study of the response of these manifolds has been carried out (11). a 

The inherent frequency response of the manifolded systems is similar to that 
_ of the pe ete tubes, as presented in Ref. I. For direct area measurements, i 


Le “that i its response up to 100 Hz is essentially flat. In all cases, differences + ol 
_ the pressure exerted on roof points or areas and the static pressure of a 
_ Digitization of the pressure signals and analysis of the data using a PDP-8/I — 
“ digital computer 0 on-line, gave the maximum, minimum, mean, and root a 
square (RMS) pressures over a 60-sec period. The dynamic pressure of the 
flow above the boundary layer (1/2 pv? ~) was measured similarly, and used uff 
to determine pressure coefficients. “Appropriate filtering and recording of the 
pressure signals coming from the various roof segments were camtied out in a 
order to allow correlation and spectral analysis of pressure loads, as detailed 


For measurements of pressure loads on a rigid model, correct wih requires 

‘the characteristics of the velocity profiles to be similar to full scale. Gus 
H 4 this, the ratio of any length from the boundary layer | to a characteristic length pe 
7 : of the model determines the length scale. The velocity profiles of Fig. 2, taken a 

7 at the location of the models without the models in place, yield a number — 
of boundary layer characteristics for consideration. An exponential law V/ 
= (Z/Z,)*, in which V = the mean velocity, at height Z; and a an 

determined constant, t, fits very well with most of the 


= von 
constant to 0.4; =a surface drag defined in 
terms surface shear stress Tt, as = Vi: ; and Z, = the roughness 
length, i.e., a quantity which characterizes the surface elements, 
> is fener considerably smaller than their physical size. Typical values of — 


: _ these pe parameters for different full-scale terrains, taken from Ref. 5, are compared 4 


to those obtained from the measured velocity profiles in the wind tunnel in 
‘Table 1. The data indicate that a 1:500-length scale is appropriate. The caret 
ponding full-scale dimensions of the models are shown in Table 2. Jensen and 
_Franck’s results (7) indicate that the experimental coefficients may well be 
for other scales | as well, as long as h/Z, is maintained 


There’ is some flexibility in the sstestion of the velocity scale and eessiiaiate 
the time scale since te/ty = (Lp/LuVu Vp ), in which M and P signify 
model and prototype, respectively. The wind tunnel speed is chosen to a 


acceptable signal to noise ratios. The velocity and time scales are then based 


# 


on the ratio of this speed to full-scale wind of interest. In this case, 


Wwe 


a. this leads to velocity and time scales of about 3:10 and 1:150, respectively. 
Consequently, the sampling rate corresponds to about 7 Sam-— 


gare to about 0.7 Hz in full scale can be detected without distortion or attenuation. — 


_ The choice of sampling period is determined to be sufficiently long to provide _ 
: Statistically s stable estimates of mean and RMS pressures and to ensure that 
measured maximum and minimum pressures provide ‘representative re estimates 
of peaks encountered during a full-scale interval of approx 1 h. The hourly "1 
_ average wind speed is a convenient, statistically-stable value which can straight- . 
forwardly be determined i in full The 60-sec used satisfies 
these ese requirements. 


TABLE 1 Velocity Pre Profile ter Full- Scale = Simulated Terrains 


ay 


Scale 
Open country 
Forest, suburban areas 


Smooth 
Built- -up wad 


Note: in. = 25.4 mm; and | ft = 305 


e 


TAME 2.—Dimensions of Buildings | in Wind Tunnel a and in Full Scal 


Height | Width Height Length | Width Built 

250 =| 250 a 


Note: 1 bd 25.4 mm; and | ft = 305 on. eee +) 


_ Fastest-mile wind speeds are also often used as the reference speed for building 
_ codes. The relationship between fastest mile and hourly speeds can be determined | 
¢ a function of the fastest-mile speed. For example Hollister (6) indicates | 
: that, for open country conditions the fastest-mile wind speed ~1.29 (hourly 
mean speed) for wind speeds of design interest. Moth 
‘The m maximum, minimum, mean, and RMS pressure coefficients for each | 
5: individual roof point [see Fig. 1(a)] have been measured for the two different 
exposures described, two model buildings (low and high), and two wind directions _ 
(normal and oblique). The Pressure coefficients are as 


pies /sec/Channei in fu scaie, an pressure uctuations wit requencies up 
j 
_ 
: 
1,700 40-200 | 040 | 0.04 
(183) | 0.01 | 016 | 0.042 
333 | 04s | 037 0.048 
Building 


4 in which p,,,x = the maximum instantaneous pressure measured over the sampling 


period; Puan = the Pressure measured over the sampling 
period; p = = the temporal n mean pressure; Puss * = the RMS pressure = =V(p-p)’; 


and q, = 1/2 p Vi, = the dynamic head associated with the hourly mean 


velocity at roof height 7. All are of with respect 


FIG. 3. —Local Mean Presoure Coefficients 


coefficients with respect to the at ‘gradient height. 
conversion factors g,/q, have been determined from the velocity profiles of rl 


-Resutts: Locat Pressure Coerricients ey 
- Contours of mean, peak, and RMS pressure coefficients for all parameter 
_ combinations considered are shown in Figs. 3, 4, and 5. Considering the symmetry 
of flow, contours are shown for each wind direction on half the roof only. 
ry However, it should be noted that the contour pattern is not absolutely symmetrical _ a 
 ohthes of the random nature of the surface roughness, some local nonuniformity = 


of the flow, and the natural statistical variability associated with such measur a 
» _ Regarding mean pressure coefficients (see Fig. 3), it can be observed that 
es is almost always suction all over the roof. Regardless of exposure, the 


q 
>. 
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‘FLATROOFS) 


there with incr increasing Geant There i is some tendency fi for slightly higher suction 
-— coefficients for the smooth exposure in the case of 45° oblique wind | flow. 
_ This has also been reported by Jensen and Franck (7) and others. sit ‘bata 
4 The comparison of mean pressure contours for the normal wind a 
high and low roofs under smooth exposure shows that reattachment of the 
_ ‘mean flow occurs further from the leading edge for the higher roof. The result © 
: is that the front part of the lower roof supports a steeper pressure gradient — 
_ and the rear half is under no substantial mean pressure load. The flow reattachment 
length is approximately a constant number ‘of building heights. The built- “up 
_ exposure results are similar, although local flow disturbances cloud the comparison 


igh 


corner. The edges” near the corner (but not at the v very y comers), are more 

heavily loaded than in the case of normal wind. This is consistent with e 
_ formation of vortices along the leading edges, which has been commonly observed 

__ The contours of uncorrelated negative peak instantaneous pressure coefficients — 
are shown in Fig. 4. The general pattern appears disorganized, particularly for 

a the built-up exposure. _ Although not shown, the peak positive pressure coefficients 
present a similar pattern, with magnitudes of the order of 0.6-0.9 for the smooth — 
and 1.0-2.0 for the built-up exposure. The broad range of the peak fluctuations, — 

2 more pronounced on the low roof for both wind directions, is very characteristic. ne 
Again, the most severe suctions are caused by the 45° oblique wind. Generally 
speaking, the extreme pressure coefficients on the high roof are either equal 
to or less severe than those on the low roof. (Keeping in mind, of course a 
reference velocity Pressures are higher for the higher 


IDS 


wind acting on a cube-shaped building a little different from the high model 
described here. Their results are similar to to the peak instantaneous suction 

in the two cases are strictly comparable in terms of the reference pressure 


a 


‘BUILT-UP Al * AVERAGE OF ¥ 


6—P —Peek ‘Suction Coefficients from Different Techniques 
Fig. 5 shows the RMS pressure coefficient contours which provide a measure 
_ of the spread of the pressure fluctuations around the mean value. As aan al 
the RMS pressure coefficients have much larger values for the built-up exposure. 


The larger RMS values occur close to the windward edges of the roof a 


nts for a 45° 
| 
Bae 
wa 
>, 
——~ 
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FLAT ROOFS 


In order to study the influence of the tributary area and the nonsimultaneous 4 

- effect of the local peak pressures, experiments with flush-diaphragm transducers 

"covering specific roof regions were carried out. Measurements were also made : 

_ by pneumatically averaging the pressures from taps corresponding to these _ 
regions. For these group measurements the roof was divided into the 16 equal “4 
_ Square regions shown in Fig. l(a). Each contains nine pressure taps. . The round 

flush-diaphragm transducer covers an area corresponding virtually to the circle 

es in each square area. The fact that the corners of each region are 
of not covered by the sensitive part of the instrument is not expected to oe 

significant errors, particularly for the interior regions where it was used. es 

_ Although not shown here, mean ‘pressure coefficients obtained from the | 

manifolding technique agree very well with those derived by averaging the mean 
¥ values from the individual pressure taps (10). For os areas loads, Fig. 6 


4 MEASUREMENT VIA PNEUMATIC AVERAGER 
MEAGUREMENT FLUSH DIAPHRAGM TRANSDUCER 


FIG. 7.—Comparison of Spectra from ‘Different Measurement | Techniques (Smooth 
Exposure, Low Model, Oblique Wind, Inner Region) 16 
compares minimum pressure coefficients for inner roof regions obtained from 7 
_ flush-diaphragm transducer measurements, the average of the individual local ’ 
peak measurements (assuming complete correlation), and the pneumatically- 
averaged results. This comparison indicates that the pneumatically-averaged — 
a compare well with those of the direct area measurements of the fush- 
diaphragm transducer although they a are always high. This is ; due in part to 
‘the overestimation of the high-frequency components, as shown by the compara- 
_ tive spectra of Fig. 7. The overestimation physically arises when the area over - 
~ which the pressure fluctuations are coherent becomes much less than the tributary L 
¥ area associated with a single tap in the grid. This phenomenon has been examined 
by a simple heuristic analysis described in Ref. 10, and is subject to analytical — 
i correction if required. Generally, the resulting overestimation is conservative. 
‘Note that the differences apparent in Fig. 6 may not be solely | due to this | 


9 igh frequency effect. — may arise from the fact that the data for the transducer 


| 
d hae 
| 
| 


- ont those for the local and pneumatically- -averaged technique were taken from 


_ experiments carried out several months apart under nominally the same conditions. a 
| Furthermore, different weighting functions are inherent in the two methods, | 
the pneumatic technique being uniform, the -transduc cer being weighted 


MEASUREMENTS WITHIN ZONE 


‘PRESSURES 


PRESSURE COEFFICIENT (REFERENCED TO q, 


«FAG. 8.—Comparison of Peak Suction Costticionts Derived Local and Area 


Measurements (Built-up Exposure) / 
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9.—Spectra of Pressure Loads on All Root Fenriater aye Exposure, Normal | 


yaa 


Wind 


both wind directions, and both models for all roof regions. Typical results are _ 
a shown in Fig. 8 for an oblique wind acting on both the high and low roof. _ 
Both the local pressure coefficients, a1 and the effective pressure coefficients for 
= the component areas shown, have been averaged over the four zones as indicated, 
in the inset. Similar results are available (10) for the other wind direction /exposure — tJ 
combinations. In almost all cases, both the average and extreme limits of the — 
7 noninstantaneous peak values given by the individual pressure taps are much | 
om more severe than the true spatially- averaged | peaks. . In some cases, the whole - 7 
‘range . of the s spatially- averaged peaks is less than the range of individual extremes. _ 
Thus, a significant overestimation of the peak total load occurs when point 
pressures alone are used. The amount of overestimation increases with the 
proximity of the zone to the edge and appears to depend on the location of 
the roof area with regard to the wind 
_ It is interesting to note that Kim and Mehta (8) measured peak, mean, ~* 
RMS loads on a flat roof in a full-scale experiment and determined a probabilistic 
) 7 - model for predicting the peak loads from the measured mean and RMS values. 
_ They applied this model successfully using the results of the present study 
(10); i.e., their measured peak pressure coefficients are predicted from the 
— mean and RMS pressure coefficients of the current study, despite — 


= effect of the peak point pressures. Measurements were made for both pened 


the different geometries of the present models and the full-scale building. 


an AND ConneLation Coerricients For Roor Component Areas 

_ Since this study was primarily aimed at ana characteristics of the effective — 

-. | pressures acting on roof areas rather than points, the pneumatic- averaging — 
technique has also been used to determine spectra of pressure loads affecting — 
the 16 equal-sized regions. Representative results are presented in Fig. 9. Each > | 
power spectral density function of pressure load is drawn in the appropriate _ 

Be . All spectra are normalized by the variance, the values of which are also’ na 
indicated, in coefficient form, normalized by the square of the dynamic pressure 
at roof height. spectra can be classified ‘into two. ‘different, although 
< clearly separated regions, i.e., the v windward region and the wake region. 
_ The pressure spectra in the windward region of the roof are similar to Ps 
spectrum of longitudinal turbulence which can be described by the universal 
ie suggested by the third writer (3). According to this expression, the | 
_ maximum energy of the wind turbulence in full scale ‘is associated with the 
_lowerv wave numbers (peak at 4.6 x 10~* cycles /ft = 1. 5 x 107 > cycles/m), 
which corresponds to the low frequency peak, at about f/V = 1.4 xO il 
. cycles /m, of the measured spectra in equivalent full-scale units. This accumula- 3 
3 tion of energy at low frequencies is clearly shown in the windward regions 
for all cases. It is interesting to compare these results with the information 
given by Kramer and Gerhardt (9) for power spectral densities of pressure 
loads on a flat roof. Although their geometry was not the same, the —- 7 


ag obtained for the windward regions in this study, the main difference being 


that Kramer ond Gerhardt’s peak occurs at t higher onan 
equivalent full-scale value of f/V ~ 3 x 10~° cycles/m. Note, however, as 
the spectra presented here correspond to area loading, whereas Kramer and 


Gerhardt’ s are for point pressures. Furthermore, Kramer and Gerhardt’ 's — 


FIG. 10 preseure Correlation Costticients for Various Root Regions 


‘two other energy peaks become more pronounced ‘the downwind 
side of the roof, reaching their maximum value in the wake region. The highest 


_~40 Hz (fh/ p ~ 0. 45) for the high model (model scale frequencies). Whereas 

the response to the gustiness of the wind can be seen over the entire roof 

on the low model, it only affects the first windward zone of the high model. _ 


| a for the apparent shift in the energy to higher frequencies. Fig. © also indicates __ 
F 


= Since the effective pressures acting on the subregions of the roof retain a 
 eemidinaite amount of the dynamic nature of the point pressures, it is necessary — 

>. define the cross correlations or cross spectra between different roof areas 


in order to determine the unsteady load acting on still larger parts of the structural 
_ system of the roof extending over more than 1/16 of the area. The effectiveness _ 
of the dynamic wind loads experienced by large structural elements is thus 
_ further reduced. The loss of correlation between two points, which is beneficial " 
_ for the roof loading, depends on their position and the distance between them. pa ie 
Correlation coefficients between the area loads on the different 1 regions of 
the low and high roof are shown in Fig. 10 for different exposures and normal Jf 
‘ wind. Each zone of Fig. 10 gives the correlation coefficients between the > 
associated shaded area and all the others. The following observations can be > 


“factors in . the development of high or low correlations. Negative correlation 
coefficients have been found only in the case of the high model for the smooth 
_ exposure wind and are small. Generally, the correlation coefficients aa : 
me rapidly for the higher roof and more gradually for the built-up exposure. 


_ 2. Although the correlation coefficients are expected to decay with j increasing 
- sepeaation between the various subareas (the rate of decrease depending on 
the position of the two subareas, as related to the wind flow direction), this 


_ does not always happen. There are cases in which the correlation of the p1 pressure 
z is higher with a more remote zone of the roof than with a neighboring» 
“region. No clear explanation is available; however in some cases it appears 
e be associated with the reattachment of the separated shear layer. od Sai 
3. With respect to the wind flow direction, the correlation coefficients (10) 
appear higher for regions located in the trace of the vortices dovelaged over 
for the 45° oblique wind direction. prs ult lap 
s Results of cross-spectral analysis show that, in most cases, the quad- spectra 
are insignificant; this implies that the out-of-phase components of the pressure 
_ load fluctations are negligible. The co-spectra given in Ref. 10 indicate that 
_ the pressure load fluctuations are strongly « correlated for ‘small separations; thus 
“a the co-spectral densities appear fairly comparable to the p power spectral densities 
for adjacent regions, particularly at low frequencies. 
ASSESSMENT OF Errective Wino Loaps ON LarGe Roor ‘Sussecrions baw 


used to derive s spectra of the total fluctuating wind force acting 0 on n larger sections 

of the roofs. Results were obtained for sections made up of multiples of the 
_ 1/16-square areas under the various conditions of exposure, wind direction, 

and height cain (10). Typical cases for force spectra acting on a roof quarter 


5 The snectra of prescire acting on the 6 sections he roo Were 
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- In general, the variance coefficients are smaller for larger areas, a 
_ they may be higher for the larger area than for some of their component areas _ 
_ if the larger area includes components on the windward side of the roof. It 
4 is also interesting to note that larger spectral peaks appear for the lower roof. 


No clear explanation of this fact is available at present, but the shoorvation 
indicates that lower roofs may be more susceptible to resonant effects, — 

Comparison oF Loca. Pressure Coerricient Resutts with Coves ‘AND STANDARDS 


_ Since pressure coefficients are determined with respect to different reference _ 
speeds in various countries, it is very difficult to make comparisons directly 
_ between coefficients. Instead, a simple illustrative comparison is made in terms — 
« pressure acting on the roofs of a low and a higher building, similar to those ae 
(wood in this study, according to the specifications of the 1977 Commentary 
_ to the National Building Code of Canada (NBCC) (simple method), the American =, 
_ National Standard Institute (ANSI AS58.1-1972), and the present results. The 4 
comparison in summarized in Fig. 12 which shows both the actual pressures, 
and ratios normalized by the current experimental results. Ratios greater than 
t 0 imply code conservatism. The wind speed considered corresponds to an — 
hourly average of 60 mph at 30 ft over an open country exposure. Considering — 


a similar storm system, i.e. » the s same ‘gradient height speed over a 


aa be ‘only 23 mph. ‘Since the simple ‘method of the NBCC does not recognize ” 

_ any terrain differences, the same values are specified for both exposures. = 
= ANSI refers all pressure coefficients to the fastest mile speed at roof | 
' height. Conversion of the hourly average of 60 mph to a fastest mile of 7 
mph was made by using Hollister’s charts (6). Different values of velocity 

_ pressures (which include gust effects) are specified for different exposures. = 
To determine the pressures from the present study, the wind speed at roof , 
= height was calculated using the experimental velocity profiles of Fig. 2. The 
pressure coefficients used from the experimental study are the worst peak suction ~ 

coefficients measured at any point in the area considered, thus do “not 


- adopted because the a areas used in these experiments do. not conveniently aline 
with the divisions of roof areas considered in the NBCC or ANSI standard. _ 4 
_ The comparison shows that the NBCC significantly overestimates the pressure 
loads, in particular for the corner region, the built-up terrain, and the lower 

building. The results ‘Suggest some adjustment should be made for the edge 
and corner specified coefficients which are very high, and that allowance for 

_ @ rougher exposure in the simple method of the Canadian Code might be 7 
worthwhile. The ANSI also significantly overestimates the negative pressures. “d 
on the perimeter zones; however the ANSI appears to underestimate some 

_of the local loads in the interior of the roof, particularly for the higher building _ 


and the terrain. This is to spatial averaging 


for an oblique wind and different combinations of height and exposure are 7 
presented in Fig. 11. Also shown are mean square effective pressure coefficients, § § 
a 
q 
| 


¥ Examination 1 of the measured pressures reveals that ¢ generally the dynamic 

component predominates. Codes usually overestimate the mean loads and 

underestimate the gust factors. For example, the NBCC considerably overesti- 

_ mates the mean (by factors of 1.5 and up), whereas the gust factors found a 
ic es were typically in the range 2-7.5 compared to the value of 2.5 F 


EXPOSURE MODEL COEFFICIENT 


SMOOTH LOW  .0264 


FIG. 11 of Pressure Loads on a Roof Quarter 


4 ft. =305Smm under the provisions of the 
— 1980 NBCC, for taller build- 


ings the coefficients specified 7 
oh tT for C are reduced by a factor — 
of 0.67 and for W by a factor 


FIG. 1 12.— —NBCC, ANSI, and Comparative Experimental Values 
specified by the Code. These large experimental gust factors may be particularly 
Significant if the structural response of all or part of the roof is susceptible - 
_ to dynamic effects. It should be noted that since the experimental coefficients — 
a used in the aforementioned comparisons are local values, spatial averaging will 


_ further reduce the loads appropriate for overall structural considerations. This 


| 
a 
ar | _psf__| ratio | psf | ratio | [ratio | psf [ratio | psf | ratio | psf | ratio | 
nsec | Low | -26.7 | 1.17 | -53.2 | 1.87 | -79.7 | 2.80 | -26.7 | 268 | -53.2| 302 | -79.7| 453 
High | | 0.98 | 1.57¢} -109.07| 1.967 | -36.2 | 1.52 | -72.57 | -1090t} 3.21t 
anst | Low |  -189 | 083 | 64.8 | 2.27 |-135.0 | 4.74 | -70] 0.71 | 240] 1.36 | soo] 284] 
High -25.2 | 0.68 | -86.4 | 1.87 |-180.0 | 3.24 | -13.3 | 0.56 | 45.6 | 1.68 | -95.0] 2.79 
Present | Low 722.8 1.00 | -28.5 | 1.00 28.5 99 | 1.00 | -176]| 1.00 | -17.6]| 1.00 
Study | High | -37.0 | 1.00 | 46.2 | 1.00 55.5 23.8 | 1.00 | -272 | 1.00 | -34.0 | 1.00 


is recognized in the NBCC for the design | of primary structural systems only 
(reduction of gust factor from vy to 2.0), a and in the ANSI for ‘areas $ greater 


_ have been undergoing revision. The ANSI revised coefficients have not been Pi 
- examined in detail, although they appear to be somewhat reduced while still 
remaining conservative. Some indication of the 1980 NBCC changes have been - 
added to Fig. 12. For taller buildings, the roof corner and roof edge coefficients | > 
have been reduced so that for a smooth exposure the experimental ratios quoted sf 
are reduced to about 1.2-1.3. For low-rise buildings, the 1980 Code hasintroduced _ 
‘Significant changes. These include: the specification of peak coefficients (i.e., aa 
ha they inherently include a gust factor); revision of the coefficients to more closely — 

a _ reflect the most recent experimental data (see Refs. ‘2, 13, and 14); reduction | 
in load coefficients for i increasing tributary area; and specification of completely _ 
* separate coefficients appropriate to the design of the primary wind-resisting 

= structural system. The effect of terrain roughness has not been introduced into _ 
the simple method of the Code; however the reference wind pressure now 

down to 6 m rather than the previous wee of 10 m. 

_ 1. The wind-induced pressures acting ona flat roof have a major oo 

i & 2. A significant overestimation of the wind load acting on roof areas cums 
when the lack of spatial correlation of gusts is not taken into account. This 
overestimation is higher for areas more heavily loaded and increases with the 


. Even within the aforementioned conservative comparison, based entirely 


on local pressure » coefficients, the 1977 Commentary to the Canadian Code on 
overestimates the suctions acting on a flat roof, in particular for the lower - 
roof and the rougher exposure considered. The ANSI A58.1-1972 document — 
ae some of the interior zone local loads (which may be compensated 


by spatial averaging), but significantly the on roof 


advice provided by their colleagues at the BLWTL. The helpful suggestion 2 
q the Validyne transducer as the most appropriate instrument for the purpose 
_ was made by R. D. Marshall of the National Bureau of Standards, of the United 
— Financial support for for this project was provided by th the National Research 


Council of Canada. 
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power er law exponent; 
3 = square root of surface. drag coefficient; 
pressure coefficient; 


frequency; 


+ length; 
dynamic head at gradient height; 
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By David J. Just” 


INTRODUCTION 


to high traffic speeds road or railway echemes | can seldom be modified 
in order to eliminate the skew of their bridges, and therefore a considerable — 
Sumber of bridge decks are constructed with a greater or lesser degree of skew. 
ool The effect of skew above angles of about 20° in single span decks can have 
a considerable effect on the behavior of the bridge, especially in the short = 
cr: to medium Tange of spans where the span and width are of the same order, a 


_ In addition to exhibiting skew, a large proportion of medium-span decks are 

of beam and slab construction, the beam system composed mainly of longitudinals 
with perhaps transverse beams at the support sections. At present the design» 

. of decks of this description | is usually carried out using a “‘grillage analysis”’ 


- (2). In this technique the deck is idealized as a framework, and equivalent 
section properties ascribed to the members so formed so that the action of 


. as in the downstand grillage (2) or in a lattice approach (2), “the attraction of 
the grillage treatment lies in its simplicity in “conception and use. The method 
of finite elements could be used, and the structure modeled completely with 

4 plate elements, but as with the more complex framework simulations, this becomes 

: impractical in terms of data preparation and computer time. Thus any develop- 
ments in the analysis of such decks should, if they are to be generally used, 


J possess the qualities of the e grillage approach, ‘that is is baste should be both simple 
'Lect., Dept. of Cy. Engrg., Univ. of Aston in Birmingham, Costa Green, Birmingham 

7ET, England. = 

Note. oes open until July 1 1, 1981. To extend the 2 closing date one month, 
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 — ACR. Manuscript was submitted for review for possible publication on March 4, 1980. 
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8001 /81 /0002-0299 / $01.00. 
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paper by the writer (4) and involves the representation of the slab by skew 
plate finite elements and the beams by unidimensional beam elements. In orde 


systems are preee: together at the corresponding nodes by members that are 
considered to be of infinite stiffness, so that no curvature | of these connections 
can occur. This assumption implies that at these nodes, in addition to Smad 


caw 


(6) Beans and plate in same plane 


elit 


sections remaining the ‘the plate and the displacements 


__ of the beam may be related to each other, so that at the plate-beam connections 
there i is only one independent node. Therefore of nodes, 
the same time the sophistication of the finite element ‘odilions t is combined 
= the simplicity of the grillage treatment. It is with the development of this 
method for skew decks of T-beam construction that this paper is concerned. 


| _ The approach described herein 1s an extension of that described in a previous “a 
al 
| 
Beams and plate in different planes ed 
—Beam-Plate Structure Idealizations 


or SLAB E DECKS 


INTERACTION © or Beam AND PLate ELements he: 


_ The stiffening of a plate by a system of beams produces a significant effect 7 
in the behavior of the plate. The effect of beam and plate interaction in one 
plane can be quite simply investigated using standard displacement method 
‘programs incorporating these two classes of element. However, the behavior : 
of structures where the beam and plate systems are in different planes i is quite : 
‘dissimilar, and if standard programs are used, the number of nodes becomes _ 
_ excessive since approaches similar to the ‘‘downstand grillage’’ (2) must be 

used. In the procedure described in this paper, this effect is simulated without 
_ any increase in the number of nodes or structural elements. oN, ae 

Consider the beam- -plate system shown in Fig. 1(a). In the idealized representa- 

tion, the beams are represented by lines along their centroidal axes, and the 
y plate by a lamina in its centroidal plane. Fig. 1(b) shows this idealized representa- 
tion applied to 2 system where the beam and plate systems are in the same 

Plane. If the | beams and plate | are in different _ planes, however, the System 


nee with R Rigid Nodal Connections 


an by connecting members, thus i increasing the number of modes. © 3:1) 46 
If, however, these connecting members are considered to be infinitely stiff, 
then the displacements of, say, a beam node can be related to the displacements 
of the corresponding plate node and therefore there is no increase in either 
_ the number of independent nodes or in the number of elements i in the structure. 
ot is thus necessary to relate the centroidal line displacements at the « ends 
of a beam element and the centroidal plane displacements at the corners of 
i a plate element to the displacements of the single independent node chosen. == 
Ria Ai Relationships between Centroidal Line Displacements at Ends of Beam Element 


and at Nodes. .—Fig. 2 depicts the centroidal line of a beam element = 


_ directions. Thus at end i the three _ connections are p,, dy and r,, while 
at end the are denoted aS Pj, how ? 


ij 
: ; = J. These two connection vectors may be considered to be each com a 
hre ano ian Om nonen one in each o he 


Since the connections are 4 infinitely stiff, the relationship 
between the beam displacements at i and the nodal displacements at I can 
be 7 in terms of p,, and) ras: 


and sila relationships can be written: | for the displacements at the other end 


be 


 islacements [Y,] at the nodes may now be related to the displacements 


the — [X] being in a set of global right hand a ain 


ucture, by the equations 


m, 


in n which the of = values of the various cosines 
by the two sets of axes; and the relationships at the node J may b be similarly 

‘Thus the relationship between the two sets of ‘displacements [Z,] and [X] _ 


and sideline the stiffness. matrix of the beam referred to its local axes Pp. 
q, r by [k,], the global stiffness matrix of the beam with rigid connections 
IK, |= = = 1A, 
The heats peri used in this paper include the effect of shear deformation, ss 
the matrices being described in full elsewhere (3). 
Relationships between Centroidal Plane Displacements at Corners of Skew Plate” 
and Displacements at Nodes. between the corner and nodal 


4 
| 
axes fo 
= [V,] [ seth & 4 


displacements for the skew plate ‘element can be assembled i in a similar meme j 


3. The skew plate finite elements used in 1 this ‘paper | were formulated such = 
that the two in-plane translations are defined along the. p-axis and q-axis and — ae 
_ the three out of plane displacements in the r, and directions. These 
elements are the ‘node, 8 x 8 skew plate. element for | plone 


a 


a 


vith Rigid Nodal Connections 


4 
stress 6) andt thes 4-node, 12 x 12 nonconf orming bending element 
q a that formulated by Dawe (1). Thus, at corner a, the equations relating the 
ae of | the s skew plate to those of the rectangular plate may be expressed — 


in which u, v, w = the displacements i in the Ps r disections, 


7 
| | 
| 
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equations at the remaining corners may ‘te formed. These 


Now consider | the plane of the element with 
comer a displaced from the corresponding node A and joined to it with a ery 
connection vector. Again this vector may be considered to be composed of 
three rectangular Cartesian components, ¢ one in each of the > q; and r directions. 


Similar relationships can be formed between the dieplacements at the rer remaining “— 
corners and those at the corresponding nodes in terms of the components of 4 
the associated rigid « connection vectors, and the complete set ad of equations m — 


to the ] at the corresponding nodes can now be expressed 
4 


Now the nodal displacements can Telated to the 


these relationships being written in full for node A as 
Jo 


| 

in which the elements of = the the cosines formed 


_by the two sets of axes. ‘The were at the remaining three nodes may 
be similarly formed. 


| 
- relationships between the rectangular plate displacements at corner a and the 
displacements at the node A may be writtenas 
‘ 
4 
| | va | | 


the between the [S,] at the 


the skew plate and the a at the nodes is by 


a equations being written dons full for comer a of the cow aie which 


= 


= 


mt form of many mc 


FIG. 5. —Model Skew Beam and Slab site 


| 
Denoting the suliness matrix of the skew plate, referred to its local axes by _ 
y 
[k,], the global stiffness matrix of the plate with rigid connections is thus 
w Beam ano Sis Broce 
— dern beam and slab decks is a thin slab connected _ 
_ tothe top of a senes of longitudinal beams. Transverse beams are placed over 
7" 
tyes 
— 


transversals ‘reduces ‘the load characteristics ‘of the deck, but 
provides a simpler and more economical construction procedure. 
ies _ A small scale simply supported 45° single span skew deck of this form was 
“thus investigated both numerically and experimentally. The geometry of the 
deck is shown in Fig. 5. The model bridge was constructed from perspex for __ 
_ which the values of Young’s modulus and Poisson’ S ratio were found to be ae 
3,176 N/mm ? and 0.335, respectively, for the 6-m mm wide beams and 3,354 _ 
N/mm? and 0.390, respectively, for the 3-mm thick slab. The beams were | 

glued to the underside of the slab with “‘Tensol”’ cement. The deck was considered | 


J 


At these points, both movement in n the direction of the 
axes of the transverse beams were prevented. Horizontal deflection normal 
‘to the transverse beam axes at these “points was, however, allowed. In the 
experimental analysis, these conditions were simulated by supporting the deck 7 J 
at the appropriate points on mild steel rollers whose length and diameter were | 


7 6 mm. The deck was loaded on the edge at midspan. —t™ CTT 


i Numerical analyses were carried out on the deck using two element meshes _ 


4 for the slab, these meshes being shown in Fig. 6. 6. The > coarser I mesh consisted — 
_ of 24 elements and the finer one of 80 equal elements. Results of the analyses _ 
performed on the deck are presented | in Figs. 7, 8 and 9. It can be noticed = 


| 
De supported at ihe cnds of cach of the equally spaced longitudinal beams. 
| 
[Vii 
| 
__FIG_ 6 —Meshes Used in Analyses of Model Skew Deck 
| q 
| 
tha er SUIS Tor defiection Show 4 Very good agreement Wi 


‘those fron rom the numerical analyses : and also that the difference between deflections 

Obtained from the two mesh sizesis quitesmall, 

- Calculation of normal stresses, as shown in Fig. 8, shows that the coarser 
mesh predicts maximum stresses of about 88% of those predicted by the finer 


"mesh, although | elsewhere the differences between these results obtained from 


(a) Deflections along arid line 1 


As 
Loading: 100 N downward on node Cl 
b) Deflections along grid line 


FIG. 
the two analyses are much smaller. ‘Fig. 9 shows typical shearing stress 
distributions in the slab, results which cannot be completely obtained in a standard 
~ In order to assess the effect of skew, a right deck of similar form spanning 
"300 mm was solved. Table 1 shows the comparison between results obtained __ 
using t the 80 ‘mesh from both skew and right d decks, where it is 


| | 
a 


: 


on 

Soffit stresses. or 
(tension) 


he. 


sis 


oso 


———— 8.66 
= (compression) 


Normal stresses on top and soffit of slab along GLC 

ne. 8. terme! ‘Stress in Model Skew Deck, in 


= 


| 
i ~ 
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Note: In both diagrams, stresses on top and soffit 

Of slab are in opposite directions 


(a) Shear stresses in between and parallel to eid lines and 2 


wer 


Parte 


9. —Shearing Stress Distributions from 80 Elmont Mesh ‘Analysis of Model 

- higher : stresses occur in the skew bridge although a a lower sonieneme deflection 


be subjected to distributed loads and also to poi«t loads acting between the nodes, 
_ the most universal form of the first category being that of self weight loading. — — 
In the dis In the method of analysis, these loadings | are represented as 


—— 

he 

du 
4 20 

| 
L 
Biber 
— 


fixed ad acting at the nodes of the and thus the 


_ Solution procedure remains the same whatever type of loading is applied. ps 

_ Considering a uniform skew plate element, of the type described in this paper, 

subjected to a uniformly distributed load, the equivalent out of plane nodal © a 


forces are as shown in Table 2. The manner in which these results are obtained 


is described elsewhere (15). 


TABLE 1. of Deflections and Stresses in Skew and d Right 


pron Deflections at X and Y, in millimeters 
— J 


(b) ) Maximum Stresses on Top of Slab and Soffit of Beam on Grid Line 2, in Newtons 


- Note: All results obtained from 80 element mesh. 1 Ibf = 4.45 N, 1 in. = 25.4 mm, 5 


qj For a uniform beam of length L under a uniformly distributed loading of — 
_W per unit length, the equivalent nodal shearing forces and moments are wl /2 a 


and WL? /12, respectively, providing that shearing deformation is neglected. 
In this paper, however, the effect of shearing deformation in the beams is 
taken into account, so that fairly deep beams can be considered in the analyses. 
In this case the results for the equivalent fixed end forces must be modified. P 
a Consider a uniform propped cantilever of length, —L, and flexural rigidity, 
EI, subjected to a uniform lateral load of W per unit length. By symmetry _ 
the equivalent fixed end shearing forces” are each WL/2 whether shearing 


sig 
| | 
| 


kate 

deformation is into account or not. In fixed 

moments, with shearing deformation ae the rotation, 6, at the propped 


2 


TABLE Out-of-Plane Nodal Forces in Skew Plate 


Plate subjected | to load q per unit area, positive 


(X/6)y 


| 


— . >. 


Center | line 
, upstand beams 


Wee 


FIG. —Vertical Deflection Profiles in Ds Deck, in 


Deflections lity center ot 


a 


in which n = a coefficient dependent on the shape of the cross section; A 
af = the cross-sectional area; and G = the modulus of rigidity. Simplifying this 
expression and replacing G by ‘E/2(1 + the fixed 


For pili sections n is 6/5, Se thus for such a section of depth d, 
the equivalent fixed end moments become had 


(14) 


‘Anatysis oF THREE Sran Continuous ‘Beam AND o Stas Brice Decx 
a _ To demonstrate the effect of self weight in a beam and slab bridge, a sini 
span skew deck of proportions was analyzed, ‘geometry of the 


- 99 ‘Swesses on top 


*AIG. 12.—Normal Stresses on . Top of Slab and Soffit of Beam along Grid Line 2 
_in Continuous Skew Deck due to Self Weight and 1,000-kN Point Load, in Newtons eg 


-@ 
E 


- edge upstands, a feature which can improve the load distribution characteristics 
of the deck (2). The deck was supported | continuously along the lines AA’ 
BB’, CC’, DD’ , these supports “preventing: both vertical displacement and 

__ The bridge was considered to be of concrete construction, the values of Young’ ‘Ss 
modulus and Poisson’s ratio being assumed to be 30 kN /mm/’ and 0.2, respectively. 2 
To quantify the distributed dead load on the deck, the density of the — 
7 was considered to be 2,400 kg/m’, and the surfacing of the bridge assumed 
to weigh 3 kN/m’. In addition to the self weight, a vertical load of (1,000 as 
KN was applied to the deck at the point indicated. 
_ Typical longitudinal and transverse deflection profiles are "depicted wo 
11 for both (1) Dead load; and (2) dead load and 1,000-kN point load, while 
= “the variation of normal stress along the soffit of the beam and the top of 
,% the slab along grid line , % and the normal stress distributions on the top and 


| 


FIG. 13.—Normal Stresses on Top and Soffit of Slab along Center Line of Continuous : 
_ Deck due to Self ey and 1,000-kN Point Load, in Newtons per square millimeter _ 


- bottom of the slab across the center line of the deck, both for the aici - 

case of dead load and 1,000-kN orn load, are shown in Figs. 12 and 13, 


‘The analysis of skew decks of beam and slab construction, although important 

_ inmodern bridge practice, can present intractable problems in design. The method 
described in this paper is an attempt to represent easily and accurately the 
geometry of such decks and to maintain data preparation and computer time 
to proportions comparable with the grillage approach. The preparation of data 
is simpler than that for a grillage analysis, as no prior calculation of section 
properties is required, and the output includes horizontal deflections and horizon- — 
tal slab shearing stresses that cannot be obtained from a grillage representation. 


| 


on bridges, in which the occurrence of corner uplift is a possibility. ad = 


1. “‘Parallelogram Element in ‘the Solution of Rhombic Cantilever Plate 


Problems,”’ Journal of Strain Analysis, No. 3, 1966. 


z. Hambly, E E. Cc. , Bridge Deck Behaviour, John 1 Wiley and Sons, Inc., New York, N.Y 


3. Just, D.J. ‘Deep Plane Frameworks with Joints, Structural 
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the Institution of Civil Engineers, Pt. 2, Vol. 67, Dec., 1979. ) Lid, 
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_ The following symbols « are used in this pa paper: ee F come 


= area of bea beam; 


depth of | beam; vf 
-Young’s modulus; 
- modulus of rigidity; 
second moment of of t beam about neutral axis; 
[K,], [K,] stiffness matrices of beam and plate referred to axes; 
natal [k,].[k,] _ stiffness matrices of beam and plate referred to local axes; 


ar = loc 
or at corner and node N of 
[U, matrices relating [ [Z,] tc to and {z, | to to to 


translations wong p-axis, q-axis, r-axis; 
= local displacements at node N of beam or plate; 
8, displacements at end n of beam or corner n of plate; 
matrices relating [Y,] and [Y,] to [X}; 
matrix relating [S,] to [Z,]; 
= load per unit length of beam; 
- vector of nodal displacements in global axes; 


_ Comparison of numerical and experimental results indicates the viability of the 
| 
t 
| 
| 
7 
q 
g 
| 
— 
| | 
‘Pret: uy 
1 
| 


ow global displacements at Nofbeamorplate;s 

Y,], [Y, i = vectors of nodal displacements of beam and plate in local 

vectors of end displacements of beam and and corner 

ments of plate in local axes; 

,6,,6,, notations about p- axis, “axis, r-axis, s-axis; 
Poisson’s ratio; and he 


* 
angle of skew. 
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ry ana during the past 10 yr (1,2). The use of such structures has ‘forced - 
the engineer to develop appropriate design 21 and analysis criteria to n meet this 
> need (3). In developing such criteria, it is imperative to compare the performance 
of the completed structure to that of the design criteria or analysis techniques. 
‘It is therefore the purpose of this paper to present the experimental results _ 
; - obtained j in October 1977, from the field testing of a two-span continuous- curved 


-composite-steel box girder, located in Seoul, Korea. These results are then 
compared to results from several analytical schemes presently being used in 


the design of such structures. 
me bridge which was tested under : oe > and dynamic live loading is located 
along the Kwang Gyo Viaduct at the Sam II Ro Interchange, in Yeon 
Seoul, Korea. The bridge, as shown in Figs. | and 2, consists of a two-span 
continuous single steel box girder 5.25 ft x 10 ft (1.6 m x 3.0 m) in composite 7 
action with a 7-in. _ (180-mm) reinforced concrete slab. The bridge has a center | 
line radius of 164 ft (50 : m), and span lengths of 112.86 ft-116.6 ft (34.8 m-38. 5 


_ m) respectively. The general plan and cross section of the bridge | is shown 
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Prof. of Civ. Chmn., Ang Univ., Seoul, Korea; also 1979 Visitin 


a written request must be filed with the Manager of Technical and Professional Publications, =) 

ASCE. Manuscript was submitted for review for possible publication on January 23, 

1980. This paper is part of the Journal of the Structural Division, Proceedings of the 
American Society of Civil Engineers, OASCE, Vol. 107, T2, February, 19 1981. 1. ISSN’ 


: 
| 
Curvep Box-GirperR BripGe: Fietp Test 
4 
. 7 


‘The bridge was fabricated it in : five ¥ various s shipping pieces a and nd spliced at four 
locations, the details of which are shown in Fig. 5. 
‘The initial bridge was designed for one lane of traffic. Also, although not 
designed for composite action in the negative moment region, shear studs were — i 
_ provided throughout the entire length of the bridge and spaced at 3. 3 = (1.0 


4 


2.—Interior Support 


is As shown in Fig. 4, the top flange is solid and has stiffening 1 ribs, as do 
= _ the bottom flange. Depending on the section location, the compression flange 
"7 has six longitudinal stiffeners, and the tension flange has three longitudinal > 
stiffeners. Fig. 4 is a typical section for the positive moment region. pe 
i ‘The web of the box also has horizontal stiffeners, located at one fifth “a 


the girder depth . Vertical stiffeners, spaced between 3 ft and 4 ft (0.9 m owe 
1.2 m) were also attached to the web along the span. 


CONCRETE SLAB ert 


: 2.S4cem = n. 
Additional full cross interior transverse diaphragms. spaced at 15 ft (4.5 m) 
and partial diaphragms spaced at 5 ft (1.5 m) were also used to inhibit distortion 


EXPERIMENTAL 
General. —The bridge response was by strain and deflection gages, 
Which w were applied at various locations along the span. The bridge was Gen <>. 


in 


| 
= 
= 


: _ Strain Gages.— There were a total of 42 strain g gages attached at three locations — 
along the box girder. As shown in Fig. 6, six of the gages were rosettes and 
co: the remaining 36 were uniaxial. The gages were attached to the interior of Le 
-. steel box flanges and webs, at ees of each span and 4.6 ft (1.4 m) 


y 
Deflectometers were also attached at midspan- as each span as shown i in Fig. Se 
7 - Two gages were used at each location, in order to evaluate rotations, and ce 
were attached at the exterior parapets on either side of the bridge as also shown 
i Recording Equipment.—The static load test data were recorded on manual 
r bridge / balance equipment. During the dynamic tests, automatic strain / i 


After all of the were ond initial 2 zeros established, the bridge 
was subjected to a known 44 k (20 T) type two-dimensional vehicle, as shown — ; 
in Fig. 8. This vehicle and one similar to it, was utilized to examine the bridge - 
under three static loading conditions. These Sere as shown in Fig. ‘ % — : 

established to induce the following: 
| Loading “i Maximum positive moment truck. 


_ subjected to a known truck loading, positioned at three locations, such that piss 
= j|§ maximum positive and negative stresses were induced. Details relative to the — - 
Loading 3: Maximum negative moment 


The tests were a single truck across the 
a crawl speed of 3 mph (5 km/h)-25 mph (40 km/h). 
; _ General.—The apis of curved box girders can be predicted by various 
P= analytical s schemes. One scheme (4), which has directly incorporated warpi 


LONGITUDINAL STRAIN GAGE 
LOCATIONS 


i 


Rosette Gage 


«FIG. 6. —Strain Gage 

- distortion of the box section has been previously used in examining the — 
actual response of a three-span curved box bridge (5), located in Baltimore, 
Md. However, this bridge did not have a sharp degree of curvature and a 


‘not have a solid flange. Thus the of the CURSGL program 


the actual Tesponse of a bridge. #2 ex 


the | plane frame STRUDL (7) computer program and the SAP finite 
computer program were utilized. 


56.43" (17. 7.2 56. 43'(17.2") | 58.3"(17. 


MEASURED D 


q 
Dial 


} FIG. Rand —Deflection Gage Locations 
ealization of Section.— —In using the CURSGL and STDUDL computer 
rograms, the “section properties must be computed. In particular, ‘CURSGL 
equires the bending inertia /,, torsional constant K , and warping constant 
.In these constants, the program we uniquely evaluate 


| 
4 
| 


such constants when the plate sizes are Also warping distortional 
considers are considered. This i is not the case in using STRUDL, 


— 

in feet | in inches* K,, in inches* Section 


407,922 585,402 


390,651 


_ Note: The section lengths are measured sequentially from left to right along the span; 


_ Sections 1 are located at the end of the spans, section 3 over the support and — 
is between section land3. = 


| 
“4 # 


d 

‘that t the FE program SAP was as employed, which would thus f a evaluation a 
of the warping condition. 
In modeling the bridge, the . effects of the stiffening ribs on the engi and a eene ottom 
— were included by summing their —- area. This total area was then 


4 
: 
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wel TABLE 2. —Defiections Data, in inches” 


CURSGL | STRUDL | SAP an. | Inside | Outside 


N 


LOADING 


$8.28" 12.37" 45.92" 
40m) (17. (3. (14.00%) 

, 


100139.921" = 1399.21" 


st. 
- 
1 
(1) 
| at | 038 | 036 | 037 | 033 | 02 | o3 =o 
Loading Positions 


“vied by the box width to obtain an equivalent — thickness. The yen 


CURSGL | STRUDL 


> 


>awroaw 


0.658 


ri 


0.02 (R) 


of 4. —Normal Stress: Loading 2, in kips per square ir inch 


Section | Cl ‘STRUDL Experimental: 


. 


‘The calculated bending and torsional of the bridge are given in 


connectors which we! 
of 
| | -053 | -012  -048 
‘Top flange | | | 
q | 
3 — 
| 
—— 
| -045 | -022 | -092 
7 schemes has resulted in the data given in Table 2. The referenced locations 


in ‘the table are as in Fig. The results show the CURSGL and 
STRUDL programs give excellent correlation with the test data. In general 
_ SAP gives greater deformations, which is probably due to the element size. E 
: Only 100 rectangular elements v were used i in modeling the bridge, , due to a limitation ¥ 
_ Normal Stress.—The mean stresses in the top and bottom flanges. and web, 
a for each case loading, are given in Tables 3-5. All three analytical schemes * 
‘show reasonable correlation with the test data. Although warping was 
be - for, it contributed only 0.5% of the bending : stress. om the stresses = 
ra s the top and bottom flanges: were uniform. — 


ss t 


4 Location ‘Section CURSGL 


g Note: | ksi = 6.9 MPa. 


” Dynamic Response.—The impact factor was calculated using both the experi- 
mental deflection and stress results. The largest deflection impact factor was — 
computed as 0.75. The largest impact, using the stresses was 0.16. 

; a _ The experimental natural frequency of the structure was observed at 2. at 2. 7 siteaaal . 

The computed value, based on the following © ia 

EI, _ 


The testing of a two-span continuous composite steel-box girder has 
: provided information as to th the induced stresses and deformations when subjected = 


to calibrated loads. In g general, correlation between three different analytical  . 


‘Experimental 
Top flange | -0.14 -041 | 
| 
q 
— 
gives a value of 2.6 cps. 


and the test excellent. Due high torsional ‘stiffness « of 
the box, warping was negligible and thus beam idealization is acceptable. —__ 
Examination of the induced stresses and deformations indicates a maximum 
; _ normal stress of 2.6 ksi and deflection of 0.37 in. Comparisons with the American 
Association of State Highway and Transportation Officials =" code 
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ae STIFFENERS FOR THIN-WALLED MEMBERS 
Thomas P. Desmond, M. ASCE, Teoman M. ASCE 
and Wister,” Hon. M. ASCE 


IntRopucTION 


Local buckling is an important and usually a mode of 
bf of thin-walled members. One way to enhance the local buckling strength | 
© economically is to add longitudinal stiffeners to the flat plate components of 
the thin-walled members. The term “‘longitudinal’’ implies that the stiffener 

i ‘a is parallel to the direction of the in-plane compressive stress. Longitudinal 
stiffeners are widely used in thin- walled components of structural shapes — 

_— cold- formed from light gage structural steel sheets. This investigation deals with — 


the structural behavior and strength of such longitudinally stiffened compression | 

The cross sections of typical cold-formed structural : shapes are composed : 
of stiffened and unstiffened plate elements. The former element is supported = - 
along” both edges, such 2 as s the flange of the hat ‘Section shown 


other, ‘such as in the flange of ‘the channel ; section also shown in Fig. (6). 
Due to the lack of longitudinal support at one edge of the unstiffened element, 

its critical buckling stress and ultimate strength are considerably smaller than - 
those of a stiffened element of similar dimensions and material properties. It 
is, therefore, advantageous to add a longitudinal stiffener in the form of - 
_ Straight lip or L shape to the free edge of the unstiffened element, thereby - \. 
‘transforming its behavior into that of a stiffened element. Such a longitudinal — 
‘stiffener is commonly referred to as an edge stiffener. The post- buckling 
distribution and buckling mode of unstiffened, | ‘stiffened, and edge-stiffened 
‘elements are shown in Fig. 


‘Univ., Hollister Hall, Ithaca, N.Y. ‘14853. 
a >Prof. Emeritus, Cornell Univ., Ithaca, N.Y. 14853. 
Note.—Discussion open until July 1, 1981. To extend the closing date one month, 
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This paper is part of the Journal of the Structural Division, Proceedings of the American 
Society of Civil Engineers, ©ASCE, Vol. 107, » No. ST2, 1981. ISSN 0044-— 
8001/81 /0002-0329/$01.00. 
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¥ In this paper, a einige fener requirement is presented for edge stiffened c compression 


elements. The requirement is given in the form of an adequate flexural a 
_ for the stiffener. The behavior of plate elements with various size edge stiffeners © 


_ was studied analytically and experimentally; however, the edge- stiffener — 


BIG. 


y= 


strengths are are formulated for elements with either adequate or less than adequate — 
size stiffeners. It is implied that the ultimate strength of a plate element with 
van adequate edge stiffener is is equal 1 to that of a stiffened element. nent. aus a 


‘Critical Buckling Behavior.—Two interrelated yet fundamentally different 


P buckling modes characterize the behavior of edge- ~stiffened elements. One is 


n conjunction wit stullener requirement, methods [for predicting ultimate 
San 
| | 
| 


stiffener in a direction to of the it is s intended 
_ to support. For this mode, instability of the stiffen = 
ad simultaneous local buckling of the flange. The second mode is s the local plate 
a buckling mode, where instability is initiated by local buckling of the flat plate _ 
elements, as though they are simply supported along their longitudinal edges. a 
Both of these buckling modes are shown schematically i in in Fig. 2%» 
relanvsly 
Buc ing Mode Local Plate Buckling | 


FIG. 2.—Critical Buckling Modes: 


Local Minimum 


tt There is no single critical buckling equation that characterizes both stiffener 
_— local plate buckling modes; consequently, two independent critical buckling 
analyses are necessary to depict the behavior of edge-stiffened pitt ° 
elements. For the stiffener buckling mode and for the local plate buckling t mode, 7 | 
the critical buckling equations are given in Refs. 1 and 9, 
derivations of these equations are adequately documented in the literature and 7 


| 
| 
= <3 
&g 
| 


_ The solution of the determinantal equation for both wesiiied modes can cnt 


; given as a relationship between the buckling coefficient k,, and the aspect ratio 


k,, and takes the form of a “garland” shape curve for. both ‘stiffener 
local — equations. Several « curves are shown for both critical buckling — 


12 
Ed. 


FIG. 5. —Minimum Critical Coefficients. 


- modes in Fig. 3, where each curve represents a different number of buckling ; 

_ waves. There are two local minima shown in Fig. 3(a). For relatively small — 
_ size stiffeners, the secondary local minimum is smaller as shown in the figure. | 
For large stiffeners, the primary local minimum would be smaller. Several 

- quantitative buckling coefficient curves shown in Fig. 4 show this point. For - 
clarity, only curves for a single half wave length are shown in the figure. al ' 


_ curve represents a different edge stiffener size for the asp: cha eo section 


| FIG. 4.—Quantitatis Stiffened Element 

4 ail kling Mo w/t 2456) 
— 
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THIN WALLED MEMBERS 

For both stiffener buckling and local plate buckling equations, the minimum iu 
buckling stress coefficients are plotted in Fig. 5 for various combinations of - 

‘stiffener size D, and flange width w. The variable D, is the flat width d plus — 

7 _ the corner radius at the juncture. From this figure, the influence of the edge a 
"stiffener on the critical buckling strength of the flange can be categorized as 
1. The flexural rigidity of relatively small stiffeners (D,/w less than about — 

0.12 in the figure) is not sufficient to support the plate; consequently, the —e2 &§ 
buckles in the stiffener buckling n mode. Here the secondary local minimum in 

Fig. 3(a) governs the minimum buckling stress. Uns jancture 

2. For moderately large D,/wratios (D, /w between 0.12 and 0.4), simultaneous 
" 5 stiffener and local plate buckling initiate buckling of the assembly. The primary — 


local minimum buckling coefficient in Fig. (a) approaches four 


( micro inches per inch) 


2 
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-Longitudi 


md Behavior: Finite Element Analysi 


however, buckling is initiated by simultaneous local plate and stiffener buckling. ; 
3. For D,/w ratios larger than about 0.4, buckling is initiated solely by local 
“plate buckling. Here the minimum buckling coefficient is governed by the local © 
plate buckling equation. Paces, was prevenied This wns 
_ For the latter case, local instability of the edge stiffener interacts with the | 
. be stiffened flange and initiates a premature local buckling of that element. 
It is noteworthy that appropriately shaped edge stiffeners other than straight 
lip stiffeners are not in general prone to local buckling interaction, and for 
uch stiffeners, excessively large stiffener dimensions do not affect adversely a 


| 
; as the stiffener rigidity incre num buckling — 


Post-Buckling Behavior.—A geometric finite element analysis was 
employed in an investigation of the post-buckling behavior of edge- stiffened 
assemblies. The finite element solution developed in Ref. 8 combines a one-step _ 
‘ incremental (tangent stiffness) approach with one iteration in the Newton-Raphson 
Mode for the solution of the algebraic equations, 
5 Results are given in Fig. 6 for an analysis of several edge-st stiffened aa. 
having flange and web flat width to thickness ratios (w/t) equal to 97.6 and 
mn stiffener flat width to thickness ratios (d/t) ranging between 2 and 18.2. For 


FABLE 1.—Cross- ‘Sectional Dimensions and Test Results = 


Test 
designation 

E (w/t) 


E-21.4-0.0 
E-21.4-1.33 


-—E-23.9-0. 
E-23.9-2.87 
 -E-23.9-5.26 
 E-45.6-5.00 
E-45.6-8. 
E-45.6-10. 
E-45.6-15. 
E-45.6-20.1 
E-45.6-24. 
E-97.6-15.0 
E-97.6-30.0 
_E-133.-9.89 
133.-15.0 
E-133.-18.2 
133.-25.0 | 
E-133.-32.0 
-E-133.-38.0 
4 E-45.6A-16.5 
E-45.6A-22.1 


E-45.6A-27.4 


4 Note: | in. = 25.4 mm; 1 ksi = 6.9 MN/m?; 


inches 


(7) 


(19.4 
20.0 
24.5 


behavior is the deflection of the stiffener- -flange juncture in a « Gneten perpen- 


_ dicular to the edge-stiffened flange component. Ideally, the flange can be treated 

as a stiffened compression element if the longitudinal stiffener is sufficiently — 2 
rigid to prevent normal deflections at this juncture, 
Two observations regarding the behavior of edge-stiffened assemblies can 


be drawn from this analysis: 


ined 


w, in d,in | tin in r,in L, in in | square 
riches inches inches | inches |inches | inches inch 4 
60 | 0.00 | 0.0747| 2.40 | 0.036 | 0.090 | 10.0 40.3 | 
. 60 | 0.10 | 0.0747} 2.40 | 0.036 | 0.090 | 10.0 | 03 -- 
_ 60 | 0.50 | 0.0747| 2.40 | 0.036 | 0.090 | 10.0 40.3 
; 2.50 | 0.00 | 0.104} 3.50 | 0.036 | 0.198 | 18.0 43.9 | 47.6 
2.50 | 0.05 | 0.1046 | 3.50 | 0.036 | 0.198 | 18.0 |M45.0 | 47.6 7 
2.50 | 0.30 | 0.1046} 3.50 | 0.036 | 0.198 | 18.0 
; 2.50 | 0.55 | 0.1046 | 3.50 | Qos6 | 0.198 | 18.0 | 60.0 | 47.6 
3.40 | 0.37 | 0.0745 | 2.55 0.113 | 180 | 410 | 
3.40 | 0.60 | 0.0745 | 2.55 0.113 | 18.0 46.6 | $2.7 
7 3.40 | 0.75 | 0.0745] 2.55 _ 0.113 | 180 | 465 | $2.7 
3.40 0.0745 | 2.55 0.113 | 18.0 39 
3.40 | 1.50 | 0.0745 255 0.113 | 180 | 49.0 | 34.0 
: 3.40 | 1.85 | 0.0745) 2.55 0.113 | 18.0 | 47.0 | 55.7 
6.44 0.066 | 3.75 0.097 | 30.0 36.6 
| 6.44 | 0.99 | 0.066 | 3.75 0.097 | 30.0 | 34.0 | 36.6 
6.44 | 1.20 | 0.066 | 3.75 | 0.097 | 30.0 | 380 | 366 ; 
| 6.44 0.066 | 3.75 0.097 | 30.0 | 380 | 36.6 
8.77 0.65 | 0.066 | 3.75 0.097 | 30.0 31.0 | 330 
8.77 | 0.99 | 0.066 | 3.75 0.097 | 30.0 | 36.0 | 330 
8.77 | 1.20 | 0.066 0.097 | 300 | 35.5 | 33.0 
8.77 | 1.65 | 0.066 | 3.75 |MMMMM | 0.097 | 30.0 | 350 | 333 
8.77 | 2.11 | 0.066 | 3.75 | — | 0.097 | 30.0 | 39.4 | 347 
8.77 | 2.50 | 0.066 | 3.75 | — | 0.097 | 300 | 39.0 | 33.0 a 
3.40 | 0.66 | 0.0745) 255 | — | 0.113 | 18.0 | 46.5 54300 7 
3.40 255 | — | 0.113 | 180 | 49.7 | 5430 
3.40 0.0745 | 2.55 0.113 | 18.0 | 49.6 | $43 
3.40 | 1.65 | 0.0745] 2.55 0.113 | 180 | 50s | 54.3 
340 | 204 | 0.07s| 255 a 
2.04 2.55 | — 0.113 18.0 49.0 | 543° 7 
4 


l. “The not remain in the post- buckling 
_ range but deflects normal to the plane of the flange. This post- ‘buckling behavior © 
is not restricted to relatively ‘small stiffeners but in general occurs for stiffeners — 
of all s sizes. Apparently for stiffeners that are above a certain required minimum 
y or characteristic rigidity, the stiffener- -flange juncture deflects at or very near 
the critical stress corresponding to k » = 4. For those that are below this 
% characteristic rigidity, the juncture deflects at a critical stress corresponding 


2. Ani increase in the flexural rigidity of the stiffener beyond this characteristic 


the flange — 


EXPERIMENTAL 


of 


designation: web, 


> 


inches | inches | i inches | inches | i — inch 7 

| & (10) 
6B-6.76 
E-45.6B-8.84 
E-45.6B-15.0 
—E-45.6B- 0 


y = column and beam specimens. (A stub cobema i is a column whose length © 


_ is sufficiently short to prevent overall buckling but long enough to permit local — 
— buckling of the individual component plates.) The stub columns consisted of si 
y lipped channels and hat | sections, » and the beams of edge- stiffened I ot 


edge stiffener was a straight ‘lip. The cross- "dimensions, material 
properties, and ultimate strength of each specimen are given in Tables | and 


— on the behavior of the ~ohe For this reason, ‘local buckling of the web, which ~ 
_ might adversely influence the behavior of the flange, was prevented. ‘This was a 
achieved by fabricating the specimens such that a double thickness web (Fig. 

7) would be fully effective. Epoxy combined with rivets or screws were used , ; 
reliably at Cornell in Previous oxperenentes work and maintained adequate bond : 
= For test series E-21.4 and E-23.9, the ends of the stub columns were peut = | 
and set against ground bearing plates prior to loading in a hydraulic testing 
machine. the other stub- test series, the ¢ sections were to large 


: : of the stiffener-flange juncture. The fact that the juncture deflects, even for 
j 1as implications regarding the ultimate strength capacity 
___A series of tests were performed to study the critical and post-critical behavior | -_ 
dge-stiff | 
| 
| 
| 2, and their cross-sectional shapes are shown in Fig.7. 
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to be und, we these specimens were welded to end plates. For the beam a 


specimens the test setup is shown in Fig. 8 


region of the beam specimens to establish a1 an | approximate ‘stress ‘distribution 


21.4, E-23.9 


| 


= diaphragm 


vA 


ews. 


across the section and to when local or stiffener buckling, or both, 
— occurred. In addition for test series E-45.6, dial gages measured deflections 
at the stiffener-flange juncture in a direction perpendicular to the edge-stiffened 


q 
4 
= 
Web 
4 
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In Fig. 9, observed critical buckling stresses obtained from 
Strain gage readings are compared with those e obtained from the critical buckling 
analysis. Experimental buckling stresses were obtained by the strain reversal — 


02 


é 


A, =" 
14° 1.6 
FIG. 10.—Post- Buckling Behavior: 
method (7). Due to the scatter in buckling stresses in test series E-133, wi 


pera is not compared with the analytically determined values. The initial 


imperfections which were more prevalent in ‘flat plate elements with larger 
w/t Tatios, are the cause of scatter in this series. Except for this 


rand Analytical Buckling 
x 


series, however, the correlation between 1 the experimental and analytically 
The post-buckling behavior of edge-stiffened elements was experimentally 
F investigated by observing out-of-plane deflections of the flange at its stiffener- 
_ flange juncture. In Fig. 10 the deflections at this juncture for three a 
_ specimens in test series E-45.6 are shown. The larger edge stiffener of specimen _ 
a E-45. 6-15 (45. 6 and | IS are the w/t and d/t ratios) showed a slight tendency 
: to deflect but retarded ‘effectively any post- buckling deflections at its juncture 


with the flange. On the other hand, the smaller stiffeners for specimens E-45.6-5 
TABLE 3.—Cros Cross-Sectional Dimensions and Test Results from Refs. 2 and 2 


teste 


lo, per square i inch P 


= 


Fic. of E Experimental and Analytical Post- 
and E-45.6-8.05 were unable to keep the flange-stiffener juncture from ~~ 
prior to failure in a direction normal to the plane of the flange. The experimental 
post-buckling behavior for test E-45.6- 8. 05 is also compared with the results 
of the finite element analysis in Fig. 11, , and they are in reasonable agreement. 
_ Two other experimental studies are reported in Refs. 2 and 12, and in Tables" 
3 and 4 the results of these studies are presented. The pertinent cross-sectional 
for these specimens is identified in Fig. | 


7 
web/t- kips t,in inches 
| 2.69 | 39.8 | 44.1 0.080 
241 [269] 4.0 | 3.6 0.080 
27.3 | 4.25 | 58.5 42.6 | 27.0 0.080 
“Determined from ultimate stress and cross-sectional dimensions. = = 
= 


In the on an effective width| are 
formulated for predicting ultimate strengths of: (1) Elements that have adequate 

5 stiffener rigidity; and (2) those that do not. In a na section, formulas 


and for assessing an adequate stiffener rigidity are interdependent. ice 


TABLE 4.—Cross-Sectional Dimensions and Test Results from Ref. 12. 


; 


4 


FIG. 12.—Cross- Sectional Geometry o of f Specimens for for Refs. 2 and te 


_ Ingeneral, the procedure for predicting ultimate strength of members containing | 
edge-stiffened elements is to treat the section as an assembly of stiffened and 
—unstiffened compression elements and then to predict the ultimate strength of 
each component element by an effective width approach. 
_ Earlier research by Winter (11) provided an experimental modification of 
ven Karman’s relationship (10) for effective width of stiffened compression — 

a elements. Later Winter (13) revised that equation slightly for design use in Ref. 

i 6.T he following i isa a generalization of this more recent | effective width equation: — 


at 
__Utrmate Strenctu Consipers 
| —_—__‘ for establishing required stiffener rigidities will be given. It will become apparent [i 
Wit rocedures for predictine ultimate speth 
he 
fitie 
Test square| kip- | L,in | b,,in | b,,in| in in | d,. in in| LS, in 
series | w/t| d/t| inch | inches |inches| inches |inches| inches | inches| inches | inches| inches _ 
4) | (8) | (9) | (10) | (11) | (12) (13) 
_ The 22.9| 3.03) 33.1 | 437 | 138 | 7.96 | 7.88 | 7.91 |0.1471] 076 | — | 560 8 § 
= BB} 28.9} 7.35] 30.2 | 428 | 138 | 3.95 | 3.93 | 3.97 | 0.0599] 0.53, | 56.0 
a, 8 32.0] 5.12) 36.2 | 339 | 138 | 7.95 | 8.01 | 7.95 0.1109] 0.77 | — | 560 § 
38.3)10.3 | 30.2 | 49.5 | 138 | 5.11 | 5.09 | 3.98 | 0.0603] 0.78 | — | 56.0 
7 — 33.5} 8.33] 37.9 | 146 | 96 | 4.50 | 6.10 | 7.95 |0.0600] 0.97 | 1.01 | 36.0 i. : 
5.39) 35.1 201 138 | 6.03 | 6.00 | 5.92 |0.1082| 0.77 | — | 56.0 
23.6| 5.33] 36.2 | 305 | 138 | 5.98 | 6.02 | 7.94 | 0.1091] 0.75 
| 19.2] s.oo} 35.1 | 107 | 138 | 497 | 497 | 3.94 |0.1070| 0.74 | — | soo 
d 
Were = 0.95¢ YJ —— | 1-0.209— yf —— 
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k, to 4 4, Eq. 1 is the ‘used in the American 


Iron and Steel Institute (AISI) specification for stiffened elements. wi wit: Pry 

By. More recently, Kalyanaraman (4) at Cornell has shown that Eq. 1 with an- 
appropriate buckling coefficient can be used to predict satisfactorily the effective — 

_ widths of unstiffened elements as well. Therefore, in the following, effective vd 


widths of both stiffened and “unstiffened elements will be this 


_ Depending on the degree of support, the edge- stiffened ‘element i is Son cell 


as wr ty stiffened, partially stiffened, or unstiffened. These categories are 


Ef of the flange equals ape of an identical flange which is stiffened by : 
webs on both longitudinal edges§ ‘ei 
For these elements, the buckling stress coefficient nt for the fi flange a 


| 


‘The letter equation for (k,,),, is based on experimental evidence (Fig. 9). It 
accounts for the interaction “of local buckling between stiffener and flange, - 
which was shown in Ref. 3 to adversely _ affect the ultimate strength of the — 
nem. If, however, the plate is adequately stiffened by an appropriately 
shaped edge stiffener other than a straight lip, (k,,),, is given by Eq. 3 for 
‘The contribution of the edge stiffener. to the ultimate wah. of the assembly ; 
- also determined by the effective width approach. For an adequately stiffened 


element, the buckling: stress coefficient for the straight lip stiffener can be 2 


which is the theoretical value { for an ‘unstiffened element simply su supported - 
one edge. If the stiffener is other than a straight lip, it is assumed to be fully — 
effective. The effective widths of the flange and stiffener are th then n determined + 
using Eq. | with the appropriate buckling coefficient. 
a Partially Stiffened Element.—An edge- stiffened c compression element is partially 
2 stiffened if the flexural rigidity of the stiffener is insufficient to permit oo 


q 
“an 
deauately Stiffened pment —An edoe-ctiffened clement ic adequately ctif- 


‘¢ _ The ultimate strength of these elements is bounded by the ultimate strength , 
of a stiffened element and an unstiffened element of the same dimensions. _ 
Similarly the buckling stress coefficient of a partially stiffened element (k,,),. _ 
4 is also bounded by the buckling stress coefficients of an adequately stiffened 1 
element and an unstiffened element The buckling coefficient f for 
the adequately stiffened element is given by Egs. 3 3 or Sas and that for the unstiffened ’ 
element is between 0.425 and 1. 27, depending upon the rotational Soon 
provided at the supported longitudinal For instance, the theoretical 


In order to design partially stiffened elements, | it to ‘provide 
an equation that predicts (k,,),, without recourse to a detailed critical buckling — 
_ analysis. The following expression is a close fit (Fig. 13) of the Avorn 
buckling stress from the critical analysis: 
which «k, buckling stress coefficient for pai partially stiffened 
i elements; /, ae moment of inertia of stiffener about its centroidal axis a 
to the flange element; J* = characteristic moment of inertia 


G. 13.—Variation of (k, _ Versus Stiffener Rigidity 


_ at which critical buckling of the assembly is initiated by simultaneous stiffener 
and” local plate buckling; and (k,,),, and (k., Vos are as previously given. The 
equation | was obtained by fitting a quadratic expression between theoretical : 
values for (k..) ps at J,/1} equal to 0 and 1. The variable is the adequate 

_ stiffener rigidity based on a criterion of critical buckling; however, when the _ 

Mange w/t ratio exceeds some limiting value, adequacy should be based on 

an ultimate strength « criterion. Since the adequate stiffener rigidity for the latter _ 
criterion, (1,).dequate> © , these criteria are not 


identical Consequently, it would. not be appropriate simply to replace I* with 


| 
| 
| 
| 
} 
| 
i 
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For this reason, the e is restated in terms of adequate 2nd an exponent 
— 1/n, in which n will be chosen such that the following equation is a close 
approximation of the theoretical buckling coefficients 


a, will be defined explicitly in the 


oh variable (k..)ps is then used in Eq. | to determine the effective width of 
partially stiffened elements. The contribution of the partial stiffener to the ultimate _ 
strength of the assembly is treated as follows. 
The magnitude of the deflection of partial edge stiffeners, as shown in Fig. 
depends largely on the flexural rigidity of the stiffener. Likewise, the 
’ effectiveness of the stiffener to resist load should also depend on its flexural __ 
- rigidity. For design purposes, it is suggested that the cross-sectional area of ~ 
= partial stiffener be reduced to reflect its reduced effectiveness i = | Tesisting ; 


This can be done according to: the simple linear 


4 


in which (A, Jerr = effective cross- area a of it behaves 


is determined Eq. 1); and d, = minimum required stiffener “igidity 

necessary to adequately the plate and will be defined in the 

subsequent section > 
a A comparison of this linear expression with experimental results (3) has shown " 
it to approximate the effective cross- sectional area of partial 
él Unstiffened Element. —An unstiffened compression n element can be considered 

: “~ a limiting or degenerate case of a partially stiffened element. The equation 
proposed to predict the effective width of partially stiffened elements must, 
_ therefore, provide a transition between the effective width equations for stiffened 
and unstiffened elements. It has been shown in Ref. 4 that Eq. 1 gives” good © 

| _ predictions of effective widths of unstiffened elements as well as stiffened = 
‘ elements. Consequently, the procedure outlined here predicts effective widths _ 7 

-_ of adequately stiffened, partially stiffened, and unstiffened elements by the 

- same equation (with an appropriate value for k,,), thereby providing a smooth 


‘transition between the sti stiffened and the unstiffened ele elements. 
ft An adequate edge stiffener, referred to in the foregoing, is broadly defined 
as a stiffener of sufficient rigidity so that the plate that it is intended to support 
‘ behaves as a “‘stiffened compression element.’’ This definition could be based 
on one of two following criteria. The first is the critical buckling criterion (CBC) a 
> for which the adequate stiffener rigidity is defined as that rigidity at which 
_ instability is initiated by simultaneous stiffener and local plate buckling. The — 
second criterion is the ultimate strength criterion (USC) for which the adequate 
| stiffener rigidity i is omnes as thes minimum | rigidity which makes the ultimate — 


| 
| | 


of the stiffened flange equal to that of an identical 1 flange 
_ CBC.—To obtain a stiffener requirement based on a CBC, the stiffener buckling — 
equation was solved for those stiffener dimensions for which the primary and — 
secondary local minimums [in Fig. 3(a)] were equal to 4. This represents that 4 
- stiffener dimension at which stiffener buckling and local plate buckling occur 
simultaneously for all aspect ratios This requirement is shown graphically 
in Fig. 14 as (d/t),.,, the required stiffener depth nondimensionalized with 
respect to its thickness ¢. For comparison the current AISI specification 6) 
san egal requirement is also shown in that figure. As expected, they are in 
_ good agreement | since this requirement was also based in part on a CBC. The | 
i specification limits its requirement to edge-stiffened elements having width | 
_ to thickness ratios (w/t) less than or equal to 60. An extrapolation of the | 
_ AISI specification | requirement for w/t greater than 60 is es showe is in — 14 
The th 


fa 


| 


Oy Dependin g upon the materia al y 


4. For such cases, an adequate stiffener rigidity is that rigidity at which stiffener - 
buckling and material occur ‘simultaneously. buckling stress coeffi- 


k= (13) 
_- The family of curves shown in Fig. 14 is the adequate stiffener rigidity for 
this range, where each curve represents a different yield stress. Since the flange 


160 

| 

| 


4 Temains s fully effective up yielding in this it does not not ‘buckle, 


the flange buckles before reaching a, “and, therefore, this 1 range is referred 


1 USC.—In the fully effective range, the CBC and USC are synonymous since — 
critical buckling (in the stiffener buckling mode) and the ultimate load at yielding 
occur simultaneously. In the following, a design equation is given for this range 7 
which is based on the CBC. Since the two criteria are synonymous in the — 
fully effective range, the design equation should also be valid for the USC. hs ann 


dad 


at . 15.—Hypothetical Stiffener Requirements: USC 
To make the proposed stiffener requirement general and applicable to any 
__ stiffener geometry, it is expressed in terms of the stiffener centroidal moment ay. 
(of inertia, nondimensionalized with respect to the plate thickness, i.e., J,/t*. 
_ Also, the requirement will be expressed asa a function of the w/t! ratio normalized 
with respect to (w/t),, in which (w/t), 
is fully effective as an adequately stiffened element. By so doing, the stiffener 


requirement in this range can be expressed as one — than a family 


| 
Se 
| 
| 
| 


‘which A = = = 166; Be = 0.46) 461; and ie 2" 


ee Eq. 14 was eon in Ref. 3 to be a good approximation of the reamtiied 
_ stiffener dimensions ‘obtained by the critical a Buckling analysis and i is valid for > 


n unstiffened plate element | remains fully effective. The values of 

and B were based on(k,),,= 9.85. 
* The format for this equation is believed to be appropriate since both the 
= width or yield stress, or both, are required to delimit the post- buckling 
behavior of the flange, which in turn dictates the degree of support required 
at the stiffened edge. 
_ For the post- -buckling range, it will be shown wee stiffenet ‘requirement, 
based on a CBC (Fig. 14), provides insufficient rigidity to predict ultimate strengths _ 
of edge-stiffened assemblies conservatively. (Predicted ultimate are 


approximately. to > the requirement based 08 on ‘the Cc BC. For comparison, “the CBC 
and AISI requirements are also shown in that figure, 2s” 
_ These simple ee stiffener requirements were chosen for two reasons: 


can be obtained that leads to conservative 


Comparison oF Previcteo AND 


itt Rg 


MEMBERS 
| 
| 
q 
: he stiffener requirements, identified as I-Vill in Fig. 15, are used to compare 
™ ‘predicted and experimental ultimate strengths. The requirement for which the 
| " predicted and experimental ultimate strengths are in best agreement will then 
be suggested as a possible design requirement. 
_ _ The requirements as shown in that figure have a common origin coinciding _ 
with the CBC requirement at (w/t)/(w/t),, equal to one and are simply arbitrarily 
_drawn_straio' of increasi emollient of these 
| ay | 
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in which Pad = _— ultimate ‘Strength determined by the procedures 


- outlined for adequately, partially, or unstiffened elements (whichever is appro- 
priate); (P.;,), = experimental ultimate strength; and (P,,,),, and (P,,,),, = predicted 
= ultimate strengths for adequately stiffened and unstiffened elements and are 
} determined by the width approach Previously outlined. 


5 —Comparison of Predicted and Experimental Ultimate in Fully 


Partially 


R,/R 
or 


0. ome | 
0.746 
0.746 


(105 


1.02, 


1.28 


0.878 


0.942 


To. 596 
49 = 0.719, 


(0.478 


28 Q 
(1.68 


Arithmetic mean 1. 1; standard deviation = 0. 259, 


0.577 


12300 


qe 
| 0.00 | 0.00 | 0.615 | 0.227 
214 | 133 0.615 | 0.532 
| 669 | 178 | 
0.00 | 0.00 0.215 
| 287 0.751 
29 | 22 | 0! 
| 269 | 167 | 0.651 
| 5.80 | 0.631 | 0.903 | 0.631 
28.9 | 7.35 
38.3 | 103 | 1.78 | 0952 | 0898 | | 0.898 
19.2 | 5.00 asa | 
_ These ratios normalize the predicted and experimental ultimate strengths such _ 
| that they provide a measure of the degree to which the edge stiffener supports 
: the flange. That is, if the ratios equal one, the ultimate strength of the edge . 
stiffened flange equals that of an adequately stiffened flange; and if they equal 4 
zero, the ultimate strength equals that of an unstiffened flange. 
Inthe following, the experimental and predicted ultimate strengths are comp 
__ using the previously defined ratios. The tests in the fully effective range 


considered first, followed by those in oie -buckling range. one In tits 
— Effective Range. —For the fully effective range, , the _ experimental and 


in n Fig. 16. In this range the adequate stiffener ‘rigidity is based on ‘the CBC 

4), and (k,,),,is determined by Eq.9 

le arithmetic means of the ratios R,/R for partially stiffened elements [/, 
and for adequately stiffened elements [(/J,) = (/, 

@ 1.03 and 1.11, respectively, with standard deviations of 0.305 and 0.259, 

respectively. In Table 4 it is seen that the -unconservatively predicted cases 
are those of other references and could be due to differences in the details — 
of testing and specimens. 

- Post-Buckling Range.—For the post-buckling range, (k,, )ps is determined by 

me and is a function of m and (J, hones Several values of n (1-5) and 4 
8 stiffener requirements (Fig. 15) will be assumed i in the ia comparisons - 


of R and R, have been extended to determine both a satisfactory  edge- -stiffener i. 
Site and a good approximation of (k,,),, for that requirement. 
_ As in the comparison of ultimate strengths for the fully effective range, the EM 
ratios Ro/Ra are re compared. To identify the value oul n in these comparisons, 7 
determine the combination of buckling coefficient and stif fener: requirement 
which best predicts the experimental ultimate strengths, it is arbitrarily assumed a 
that the predicted ultimate strength agrees satisfactorily with the experimental 


The number of tests (both: and partially. etiffened) that are within 


| 


—— 


"4 
this interval are given in Table 6 for each said of n and stiffener nl 
_ From this comparison, it is apparent that ultimate strengths are better predicted 4 
- for requirements V, VI, and VII, since the largest number of tests satisfies . 
the aforementioned inequality for these requirements. These requirements and 
their respective values of nm are identified in Table 6. There is no apparent = 
advantage in selecting requirement VI or VII over V; therefore, it appears 
ate requirement V (with n equal to 4) is ‘suitable for design, since it is the 
- optimum of the three with respect to weight or material cost. The requirement 
(19) 


the proposed ‘stiffener requirement and "design for these elements. 


6. —Number of Ultimate Strengths Accurately Predicted: Post- Buckling Ri 


‘Stiffener requirement r 
= ks 


VII 


ier 


largest of predicted strengths. 


_ This can be demonstrated by the arithmetic mean and standard deviation we _ 
: ‘the ratios R,/R, for only the partially stiffened elements. From Table 6 the 
best value of n ‘for requirement V is 4, and in Table 7 the arithmetic mean : 
and standard deviation for the ratio R,/R, are 1.03% and 14.7%. These are 


quite satisfactory and show that experimental ultimate strengths of partially a 
stiffened elements a are well Predicted by the proposed design procedures. vs 


= 
| 
Since there are no provisions for the design of partially stiffened elements _ 
J 
q 


THIN WALLED MEMBER 


and R ratios, graphically compesed in Fig. 17, 
; substantiate the choice | of of requirement v based « on the Statistical evidence. It - 


ie clearly inadequate for edge- stiffened elements with large (w/t)/(w/t), ratios. = 
As shown in the figure, the predicted ultimate strengths based on requirement __ 
and the theoretical variation of buckling coefficient = 2) are infair agreement _ 
wi the experimental results for test series E-45.6, E-45.6A, and E-45.6B but a 


TABLE dt of and Ultimate Strengths in Post- 


Buckling Range: Stiffener Requirement V 


— 
Stiffened; 


=(, 


0.507 2s 1.23 
0.877} 1.53 | 1.29 


“Fa 


E-45.6A 
E-45.6A_ 
E-45.6A 
E-45.6A 
E-45.6B 
E-45.6B 


“Arithmetic mean = 1.46; standard deviation = 0.426. a see 
° Arithmetic mean = 1.15; standard deviation = 0.185. 

Arithmetic mean = 1.07; standard deviation 0.155. 
‘amen mean = 1.03; standard deviation = 0.147. 


* Arithmetic mean = 1.01; standard deviation = 0. 143. 
‘Arithmetic mean = 0.936; standard deviation = 0.138. Atoouaey 


are in rather poor agreement for test series E-97. 6 and E-133. This i is not surprising - 
_ since the former tests have (w/t)/(w/t), ratios of about 1.5 which is not far 
_ from the transition between the fully effective and post-buckling ranges. They 
are, therefore, not affected significantly by changes ir in the assumed stiffener i. 
requirement. Significant deviations of predicted from experimental ultimate 
strengths, etieats eames for the latter test series in which (w/t)/(w/t), = 


5.00} 0.167 | 1.50 1.01 | 0.917) 0.869) 
| 8.05] 0.503 1.23 | 1.19 | 1.17 
| 0.860 10.810] 0.917] 0.893] 0.877] 0.877] 0.877 
| 9.94) 0.478 2.67 | 0.712) 1.32 | 1.08 | 1.00] 0.971) 0.952} 
| 9.89] 0.366 | 3.46 | 0.581] 1.37 0.909} 0.870] 0.840 
E-133 15.0 | 0.960 0.974 1.00} 0.990] 0.990) q 
«E133 18.2 | 140 | |0.908; 0.908 
8.87] 0.348 | 1.52 | 0.768) 1.84 | 1 1.28 | 1.22 | 1.18 | 
| 0.919) 1.47 | 1.30] 1.24] 1.22] 120) 
* 
6.76) 0.316 | 1.58 | 0.752) 1.89 | 1.36 1.20] 
(8.84 0.600 0.721] 1.07 | 0.952] 0.926] 0.885] 0.885 
| 20.0} 3.12 om) | | | | | esn 
| 


4 
e 


out 
tee.) 
540 60 80, 204 06 08 | 


02 04 06 08 1 40 02 04 06 08 
/ 


= 
67 and 3.46, Tespectively. For these tests, the predicted ultimate strengths 


the stiffener requirement, as recommended before, is improvement of the 


effective width approach presented herein has been shown to Predict 

re either 
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adequately or partially ‘stiffened. For those flange widths ‘that 1 remain fully 
: _ effective (as an adequately stiffened element), the CBC and USC are synonomous; 


(hes, a stiffener requirement based on a CBC leads to satisfactory a st 


are required tc to , provide adequacy for the latter cr criterion compared to the former. © 
s Thus, for design it is suggested that in the fully effective range the adequate 
stiffener rigidity be determined by Eq. 14 and in the post-buckling range by 
Eq. 19. For calculation of ultimate strengths, Eq. | with an appropriate buckling | 
coefficient predicts conservatively the effective widths of adequately stiffened 
‘Partially stiffened, and unstiffened elements. 
_ One aspect that was not fully investigated in this "study was local buckling 
_ interaction between web and flange elements. Both the analytical and ane — 
_ work were limited to those assemblies in which local instability is initiated 
in either the flange or the stiffener. The analytical model was chosen so that be 
_ web buckling would not initiate local instability of the assembly, and influence 
the post-buckling behavior of the flange, and thereby, affect a parameter study __ 
to assess stiffener adequacy. However, the results of a preliminary study recently jo 
conducted by the second writer (5) indicated that stiffener requirement V and 
the effective width approach developed here may also lead to conservative 
predictions of ultimate ‘Strengths for edge-stiffened assemblies where local 
- instability. is initiated in the web. At the present time, this subject is being ; 
studied at Cornell University, 


7 This research was sponsored at Cornell University by the AISI. The valuable 
cooperation of the Sheet Committees, and the cognizant task group is ‘gratefully 


acknowledged. 
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cross-sectional area of edge stiffener; 
= - effective cross-sectional area of adequate : size edge stiffe 
effective cross-sectional area of partial edge stif fener; 
‘stiffener flat depth plus inside corner radius; heh 
= ‘flat width of edge stiffener; 
effective flat depth of edge stiffener; 
= Young’ s modulus, 30,000 ksi (210,000 NM/m 
moment of inertia of edge stiffener taken about am -centroidal q 
stiffener centroidal moment of inertia at which 5 
Tequired minimum stiffener of inertia necessary to 


k = stress coefficient; 


buckling stress coefficient for edge stiffener; 
k buckling stress coefficient for flange; 
= buckling stress coefficient for flange that is en stiffened 


buckling stress coefficient for flange that is stiffened 


= buckling stress coefficient for flange that is 
= height of mn Sepported of ‘beam 


compressive load of stub coleman speci- 
(Pardon = predicted ultimate of stiffened element; 
= predicted ultimate strength of adequately, partially, or unstiff. 
P q t 
ened element (whichever is appropriate); 
= experimental ultimate strength, bending moment, or compressive 


(Pa 


= 
4 
| 
i 


ult ) us 


predicted of f unstiffened elemen element; 
= normalized predicted ultimate strength; 
inside corner radius at juncture of component ‘elements of test 
flat width of flange; 
= limiting flat width of flange below which it i is - fully effective 
as adequately stiffened element; 
- = limiting flat width of flange below which it is fully nee 
= flat width of wes 


Sor 
yield stress, in kips per square inch (1 ksi = 


OF she 


the 


tor 


| 
| 
| 
= stress, in kips per square inch (I ksi=6.9NM/m?); 
| 
| 


| 
é 
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usual for analyzing the buckling of 
_ Structural members assume that the cross sections of the member remain > 
undistorted. However, when a member is thin walled, - distortion may occur 
during buckling, and consequently the resistance of the _member to buckling 
= _ may be reduced (9). This distortional buckling occurs in a mode in which lateral _ 7 
_ deflections and twists are combined with general distortions of the cross-sectional oy 
shape. More generally, distortional buckling modes have been shown (8) to x 
_ be intermediate between the local buckling modes which occur in very thin-walled 
members, and the lateral buckling modes les predicted for members with one i 


_ walls by ignoring any distortion. 


One method of analyzing distortional buckling i is to determine the lowest 
eigenvalue of the governing partial-differential equations of equilibrium of the — 4 
plate elements which form the member. However, these differential equations 4 
are usually too complex to solve, and in most cases certain simplifying assumptions ~ 
have to be made. In one such study, Goldberg, Bogdanoff, and Glauz (7) analyzed 
the critical loads for simply ‘supported thin-walled beams of arbitrary cross s 
section. They reduced the partial-differential equations to a system of eight 

_ first-order simultaneous ordinary differential equations, and the lowest eigenvalue — 

was obtained by an iterative numerical computer solution, 
x Rajasekaran and Murray (12) presented a finite element analysis of the coupled 
local and lateral buckling of wide-flange beam-columns. Their analysis  superim- 
posed the displacements due to out- -of-plane flexure of the plate segments on 
the displacements which occur in undistorted members. Although this study 
reduced the complexity arising from the close finite element subdivision, it 

; was still necessary to use eleven cross- -sectional degrees of freedom at each 
node to specify the element displacements. A more recent study by Akay, E 
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- several divisions of plate elements are used over the depth of the section. a 
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Jehneen, and Will (1) of the lateral buckling of I-beams idealizes the flanges 

as beam elements, while plate elements are used to represent the web. This — 
formulation makes it possible to predict local buckling of the web, provided - 


Unfortunately, this again leads to a comparatively large number of degrees” 

_ The finite strip method (6) is a on of the finite element ‘method, — 
‘which in certain cases can greatly reduce the number of degrees of freedom — 
y of an element or strip, as only one set of cross-sectional degrees of freedom 

are eee This method ha has been used in a recent application (8) to demonstrate | 


#) sd gam 
bos 


——— Element Actions in the yz Plane 
the interaction between the local, distortional, and lateral buckling modes of 
simply supported beams in uniform bending. However, this method cannot be 
so simply applied to other support or loading conditions. rat 
_ In this paper, an approximate finite element method is presented for analyzing 7 
“the distortional buckling of doubly symmetric I-section beam columns. The 
"stiffness and are developed in terms of six nodal displacements 
of the cross section, , and relatively few longitudinal elements are required 7 
achieve satisfactory convergence of the solution. Thus this method permits the a 


| 
| 
| 
= 


economical computer study of members" under \ various conditions of loading, 

In order to achieve this computational economy, it is to 
that the flanges remain undistorted. Because of this, local buckling loads cannot _ 
be accurately predicted, and must be computed independently (8). This approach ~ 
ies been shown to be satisfactory in a recent study of the distortional buckling © 


f simply supported I-beams i in uniform rm bending 


| 


re recto Deformations in the Plane of the Cross- Section. 


4) 


» Flange al a6 te) bae {May Sab 


FIG. 2.—Element Displacements 
Following the. of the finite element method, the results of computer 

studies of its convergence and accuracy are presented. The method is then 

4 p ad 


used to demonstrate the effect of different end conditions on the distortional 
| 


weer 

General. —The finite clement formulation in this paper 
_distortional buckling is similar to that described by Barsoum and Gallagher 

: (2) for the lateral buckling of beam columns with rigid cross sections. This 

method requires alengthwise subdivision of the beam column into one-dimensional Lhe | ; 

ni elements, for which the stiffness equations for buckling are developed. | 

_ The axis system of each element is shown in Fig. 1. The z-axis is the longitudinal 


centroidal axis, s, while the x and y axes are the — and minor principal axes, 


357° 
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‘The actions are shown in Fig. “1, and are re applied initially in 
plane of symmetry yz. The initial action set is P,, M,, V,, P,, M,, V,. These > . 
actions are increased proportionally by the load factor \ until buckling occurs. “= 
Longitudinal Variation of Displacements.. —For the analysis in this paper, the 

ctl equal rigid 
“flanges ‘connected by a a flexible v web which may distort in the | plane of the — 
_ cross section. The displacements of a cross section of the beam column yall a 
are shown in Fig. 2. _ When the beam column buckles, its flanges deflect u, 
and u,, rotate u/and and twist where (’) indicates differentiation 
with respect to 2 The lateral displacements ‘of the flanges are assumed to be : 


ag 


a, L+a,LZ+0,LZ* +a,LZ’; 
a,L+a,LZ+a, +a 


while 1 the flange twists may —" expressed a as oh 
hd which rer represent linear departures from the cubic variations of the average twist © 
= 1 and | 3 may be represented i in matrix format by 7 


r 


"expresses the degrees of freedom of the nodes | and 2 of the Sener a 
using Eq. 4, the vector {8} can be expressed in terms of the coefficients {a} 


‘The matrix is | given in Appendix I. 
Cross-Sectional Variation of Displacements.—Since ‘the flanges are assumed 
rigid, the deflection of any flange point x ‘is 


OF 


it is assumed that the web becomes a cubic curve, so that the 
of any web point y can be written as oat 


a | 


; wy=0 


coef: ficients {a,,} be obtained by substituting Eq. 12 into the identities 


and N,, and {a,,} are given in 


OY J F Fat 


‘The matrix [C,, is in “Appendix L ‘Thus by Eqs “15, 4 
and 8 into Eq. 


the vector of (8) by 


IC. 
[C] '-(8) 


ww 


‘is given in Appendix 


4 _ Element Stiffness Matrix.—During buckling, the beam- column element deflects ‘ 


twists. The i Ui the strain during can be 


] dz 


hich U,, 


dx dz 


_ where U,, = the flange in-plane bending strain energy; Ur, = 
out- heating and and twisting str strain energy; and U,, = 


U,= the strain energy 


dxdz 


_DISTORTIONAL BUCKLING = #359 & 
(u 
| 
| 
2 
(19) 
« 
| py 


2—— 


can ‘obtained from E Ea. 1, in 


while do webvgeneralzed stress vector is given by 

2 The stiffness matrix [k] ofa beam element can be derived by 


Eqs. 19 and 26 apply for all (8), leads to the stiffness matrix 


: 
The matrix tk. ] is given in Appendix IV of Ref. 5. 
_ Element Stability Matrix.—The beam-column element is subjected to bending 
moments, shear forces, and axial forces as shown in Fig. 1. The decrease in 
energy of these actions ‘during buckling is givenby = 
j 2] 


~ 
‘ 


T 


where 7 = the work - associated with the in- mn flange deformations u; 


a 
q ; 
Veq = Whe WOK associated wi € out-ol-piane ange Cerorma 1005 
element is derived from 
| 


- decrease in potential energy V by a procedure analogous to that « described for 
stiffness matrix in the previous section, and is given by Woke 
sy T wiles with the kody 34) 


The matrix Is. is given explicitly in Appendix IV of Ref. Ro! pisne 
Stiffness and Stability Matrices for Beam-Column.— 


—tThe finite element 
_ tional buckling analysis of a complete | beam-column: can be represented in matrix 


ug in ony [K] a9 —y = the stiffness and stability matrices of the complete Jf 
beam-column, respectively; \ = aload factor; and = the vector of out- 
no nodal displaceme 


te) 


tiffness and stability matrices [k] and [g], respectively. Since each element 


is aligned with the beam- column, only a simple consideration of equilibrium — 
— 
compatibility at the nodes is required. 


values of which yield a nontrivial solution for in Eq. 35 are the 


| 
a 
ite Element Solutions | 
global matrices and [GU] can be assembled from the individual element 
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this paper, the sigue problem is viadived using the routines employed i in 
7 Ref. 8, which are based on the direct methods presented by Bishop et al. igi > 
_ The accuracy of the buckling solutions obtained by the method presented 
in this paper can be demonstrated by a comparison with known solutions for — 
sage supported I-section members subjected to either a uniform moment M, 
a uniform moment gradient M/I, or an axial load P. A member of length / 
= 2 meters and of the cross section shown in Fig. 3 was used for all of ""* 


Local buckling 
cutoff 


= — Classical Rigid Web- q 


Solution (Ref.13) 


0.0001 


The finite : strip of Ref. 8 ‘the comparisons sons to be made for 


the buckling of a beam in uniform The 


y be used for the flexural buckling for a there is no distortion 


| 


63 


of the web to approach infinity. 
a The variations in the errors of the finite element distortional buckling solutions. 


The ¢ errors are largest for the beam under high ‘moment gradient, when greater 7 


_ difficulty is encountered in using cubic polynomials to approximate the maneen 
curvature buckling mode. 


FIG. 5.—Degrees-c of- Freedom at Beam Ends 
where the aa stress F, is plotted against the beam length 1. In the finite 
a strip analysis, both the compression flange and web were divided into four 
- strips, and the tension flange was divided into two strips, with the transverse 
displacements of the strips being represented by cubic polynomials. On the 
other hand, the finite element method described in this paper assumes — 
_ the flanges remain undistorted, and that the web can be represented by a single 
element whose out-of-plane displacements follow a cubic curve. The 
a agreement between the two solutions in Fig. 4 for the range of beam lengths 
re corresponding to distortional and lateral buckling justifies these two assumptions. 
Errects of Eno Conoimions on Distortiona: Buckunc of I-Beams 
& The finite element distortional buckling analysis has been used to study am 
_ effects of different end conditions on the elastic buckling of single-span I- beams = 
a The analyses were carried out for the eight conditions of end restraint shown ; 
in Fig. 5. These end conditions range from total cross-sectional restraint (A) 


Tabulated or closed-form distortional buck 
: in 2 beams under moment gradient, and so the stiff 
in this paper were tested by comparing the Ve 
; i solutions for beams with rigid webs (2,13). In this case, distortion of the cross _ 
urther Comparison Iinite clement Duckling solution with the finite 
strip solution for a beam with ‘significant is shown in Fig. 4 
ho 
= 


and E* correspond to reflections about the 
-axis - of the end ‘restraints B, D, and E, respectively. Each particular beam a 
Z assumed to have identical restraints at both ends, and the flanges of all 
= are free to rotate in their own planes [Fig.2(b)}). sts 


The differences between the classical (rigid-web) buckling stress F,, and the | 
actual bu buckling stress F, for the eight end conditions is shown in Fig. 6, where — 
the percentage reduction from the classical buckling stress F. ob is oe against 


Fop/Fp- 1) 


( 


‘rror 4 


Pe 


tut 
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6. 6- -Effect c of E End Conditions on Elastic Critical Buckling 9 Stress 


increase as the beam slenderness decreases, and the effects of web > slenderness- 
: become more important. The reduction in the buckling stress increases as the 
- degree of end Support decreases. The reductions for beams where lateral 
2 _ displacement of one flange is allowed (D,E) are significant, because that flange — 
is restrained only by the stiffness of the web. The buckling resistances of beams 
with end conditions B, D, and E (for which the compression flange supports 
are incomplete) are lower than the corresponding resistances for beams Bt, 
D*, and E* whose compression flanges are completely supported. iieam, 
these differences reduce as the beam slenderness increases, and there are no 


differences | for long beams. again ra burs 


= 


5 a Curve A shows the effect of web flexibility alone, while the other curves” 
q 
| 


DISTORTIONAL BU 


ment of the compression flange (D,E) may vary with the member length, as — ‘i 
_ Shown in Fig. 7 | for end condition E E. For short-length beams where the flange » 


in-plane stiffness is high, buckling occurs in an antisymmetric mode in which ts 
the compressive —— deflects laterally as a near-rigid body [Fig. 7(a)]. The 


by the large ja between the top and flange twist angles 


ElFo 70 


FIG. 7. —Beam Buckling Modes for End Condition "a 
. For longer beams, ane antisymmetric ‘mode would require 


this, the beam prefers to buckle in the symmetric mode shown in Fig. 7(b) _ 

_in which comparatively little strain energy is stored in the web. It can be seen o 
that the end displacements and twists in this’ symmetric mode are ‘similar 
those required by total cross-sectional restraint (A), and so it can be concluded ; 
that the buckling modes for with end restraint will approach 


* >— buckling modes of beams with end restraints which allow lateral displace- _ ; 

@ | 
Ye 
vin 

| 


those of beams with total ond restraint as shown by the oie stresses indicated 
A finite element method of analyzing the elastic distortional buckling of I-beams f 
_ with flexible webs is presented. Because the flanges are assumed to be rigid, | = 4 
the local buckling behavior of short-length beams is not predicted by this finite . 
element analysis. Local buckling must therefore be checked Pail a 


a design chart developed elsewhere 9). 


for particular cases of loading and end 
ae The effects of web distortion on flexural-torsional buckling are increased 
for short and’ intermediate length beams where end displacements or rotations — 
of t the flanges are allowed. For such beams, the buckling stress. | decreases as 
the number of cross-sectional degrees of freedom at the beam ends increases. 
This reduction in the buckling stress is greatest - short- length deep beams | 
with high ratios of flange-to-web thickness. a 


_ This paper forms pa bility | 
- structures being carried out in the Department of Civil Engineering of oe 
_ University of Sydney. The calculations were performed in the C. A. Hawkins 
_ Computing Laboratory on a PRIME 400 minicomputer. Funds to purchase this s 
system were provided by the Civil Engineering Graduates Association. The 
writers are appreciative of the advice given by G. J. Hancock of the University 


of Sydney, and of his permission to use his computes — routines. — ae 


Longin Vertation of Displacements. .—In Eq. 


= 
| 
) . 
are 
[M] = —M,| —M,| M, 


“In 163 and 


—HZ + HYZ + HY*Z- -HY*Z 


Cross-Sectional Variation of Displacements.—In Eq. 12_ 
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y ‘The following are used in paper: 


Ay area of flange and total cross section, respectively; 


element dimension matrix for web strains; 
coefficient matrices associated with displacement functions 
for web and total cross section, respectively; on a Ee 


Young’s modulus of ‘elasticity 
buckling stress; aid ove 
= flexural-torsional buckling stress for rigid- reection beam; 
cubic function of z in web-displacement curve; 
beam-column stability matrix; 
4 


‘element stability matrix; 


distance between flange centroids; 


= major axis second moment of area; 


= minor axis second moment of area of flange; ia bong ce 
beam-column stiffness matrix; Se 4 
element stiffness m matrix; 
L of element; 
length of beam colum 
= displacement function matrix for flanges; 
[Nw], [N] displacement function matrices for web; 
Nye. A = axial forces on element; 4a 
P, = critical load for column buckling the y-axis; 
= total strain energy of beam-column element; aan 
= contributions to total strain energy; 
—_ displacements of top flange, bottom flange, and web; 


M.. **Theory of Elastic Stabilitv.”’ cGraw-Hill 
ve has 
the tag 
.& 
lg 
— 


= total work done during buckling of clement; 
contributions to total work done during buckling, = — ae 

end shearforcesonelement; = 

—_ ad top flange, bottom flange and web in y 

= 2/L; 


> 


7 


4 


vector of flange displacement coefficients; 
vector of web displacement coefficients; 
= vector - of nodal displacements for beam column; 
"vector of element nodal displacements; _ 


{o} stress vector for 
t = shearstressinweb;and 
‘Fotations of top and bottom flanges about 
‘giig 

doe = 


= 
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DISTRIBUTION AND DovusB_e SKIN WALL 
By Raymond Ting M.ASCE 
sottudisialt beck sit ai goiblind 


double ckin insulated light gage metal wall systems have been 


widely used in buildings for years. The structural components of this type of 

wall system consist of an interior panel referred to as the liner, an exterior 
yy referred to as the facia, and a series of transverse subgirts connecting I 

_ the facia and the liner together. _ The structure is an indeterminate structure 

> with compatible deflections existing 2 at the subgirt locations. The stress distribution | 
2 the fastener load distribution depend on the wind load distribution between _ 
the facia and the liner. It is essential to determine the wind load distribution | 

_ between the two skins for the design of a double skin wall system. Unfortunately, 

; there has been no available reference which provides the specific guidelines 

for determining the wind load distribution b between n the tv two vo skins. The existing 
“design: data “published by various panel manufacturers are based on either of 

_ the following two assumptions. The first assumes 100% of the wind load action 
on the facia. The second assumes that compatible deflection between the two a 

_ skins exist at any point along the span. The second approach is equivalent | 
.* the assumption that the wind load distribution between the two skins is 


irrational. Most of the design engineers have not been aware of the 
in the basic wind load distribution assumption. The purpose of this paper is 
to make the reader aware of the design implications of the wind load distribution 

. assumption and to highlight the need of wind engineering research i ‘in this area. "= 

_ In this paper, the stress" distribution characteristics | in ‘relation to ‘the se 
“load distribution assumptions : are illustrated for a two equal span ‘structure and 

; the importance of the wind load distribution assumption is quantitatively illustrated — 
by an example structural sensitivity analysis on a typical double skin wall panel _ 
construction. The general trends of wind load distribution are explained by _ 

7 the general characteristics of the wind behavior and the wall construction 7 
— Principal, Advanced Engineering Services, 100 North Meadows Drive, Wexford, Pa. 4 
-Note.—Discussion open until July 1, 1981. To extend the closing date one . month, 
fl written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on February 26, 

_ 1980. This paper is part of the Journal of the Structural Division, Proceedings of the 7 

peas Society of Civil Engineers ©ASCE, ‘Vol. 107, No. ST2, February, | bt ISSN 
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working stress design principles are no 
ft for a double skin wall system. The justification of the design adequacy according» a 
: to the working stress design principles hinges on the basic design assumption 
of the wind load distribution between the facia and the liner. ad ae 
‘The current wind engineering research (1) has been devoted to the wind load — 
‘distributions in relation to the building shape and the location of the wall ‘relative 
to the building. There is no available literature that determines the wind load 
distribution between the facia and the liner at a particular wall location of 
- building. In the current design practices, the wind load distribution | ae 
| the facia and the liner has been handled by the engineering judgement of the 
design engineer or the panel manufacturer. Unfortunately, there have been 
diversified opinions in this engineering judgement and in many cases the design — 


_ adequacy of the wall system becomes a major dispute in the court for wall 


In lieu of the lack of supporting data for the determination of the wind load 
"distribution, the wall system can be designed for the ‘most critical wind load . 

_ distribution. However, _ this approach | would lead to a total confusion in the 

= §6€©6=- design philosophy due | to the different design criteria. For example, a 100% 

j load distribution on the liner might be the most critical condition for the liner 

_ stress, while a 100% negative load distribution on the facia would ey 
_ be the most critical condition for the facia to the subgirt connection design. 
In n many cases, the assumption of the most critical load could inflict a a severe 

economical penalty due to the structural sensitivity ora practical design p problem. 4 
For example, if 100% of the negative wind load at the building corner is assumed y 
to act on the facia, in many cases the regular connection method using metal 

_ screws could become inadequate due to the limited number of available flutes 
in the facia panel for the screw fastening and a a costly alternative fastening — 


A typical double skin wall panel construction is described as follows. The 
liner panel is normally a tray- -shaped_ profile with a wide flat plate and two ae 
Gam roll-formed sidejoint lips. on the two sides to engage the adjacent 
panels. The cover width of each liner panel normally ranges from 12 in. (304.8 
mm)-36 in. (914.4 mm). The depth of the sidejoint lips normally ranges from 
q in. (25.4 mm)-2 in. (50.8 mm). A typical liner profile is shown in Fig. 1(b). | 7 
_ The facia panel profiles varies ‘greatly with the ‘structural requirements ¢ and 7 

architectural appearances. The majority y of the facia profiles has a a roll- formed ea 
fluted section with a depth ranging from | in. (25.4 mm)-4.5 in. (114.3 mm), | 


faa a rib spacing ranging from 6 in. (152.4 mm)-12 in. (304.8 mm). A typical 


facia profile is shown in Fig. 1(a). The commercial panels are - commonly 
manufactured from steel or aluminum sheet coils of thickness ranging from 
0. 02 in. (0.508 mm)-0.06 in . (1.524 mm). The subgirt is normally a hat- -shaped — 
or a z- -shaped ‘member with a depth ranging from 0.5 in. (12.7 mm)-1.5 in. 
Oo 1 mm) and is commonly made from steel sheet of thickness ranging from a 

(0. 048 in. (1.219 mm)-0. O75 in. (1. 905 subgirt profiles are shown 


4 
i 
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ESKINWALL 


10). ‘The construction are commonly used in the 


field: (1) Fasten the liner ite ‘to the supporting girts using metal screws; 
_ (2) place the insulating material in the liner trays; (3) fasten the flange of the - 
Subgirt to the liner lips using metal screws; and (4) fasten 


1 of Double Skin Wall 


he 


i FIG. 2. —Typica | Double Skin Double le Span Wall Sytem 


7 airtight. The relative airtightness is an important factor in the re realistic evaluation 
> of the wind load distribution. A typical double skin wall construction is shown 
sin Fig. 2, except that the sealant and the insulation are not shown for clarification — _ 


the liner construction is reasonably airtight, while the facia construction is not S 


q 
| 
| 
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‘The structural sensitivity with the respect to the wind load distribution is — 
r affected by both the boundary conditions and the relative panel properties between | 
the facia and the liner. In order to obtain a quantitative feeling, a two-continuous 
equal span structure represented in Fig. 2 is selected for the following study 
- with the considerations of the midspan deflection, the facia moment over the 
center support, and the liner moment over the center support. The following 
aq 1. The center girt i is rigid and idealized as a pinned support pene “3 
2. The end girt supports are rigid and idealized as roller supports. _ “ihe ; 
i 3. The midspan subgirt is rigid and thus the compatibility in deflection exists 


between the facia and the liner at the 
4. The subgirts over the girt supports are rigid and the — lips over the | 
girt supports are rigid. Thus, the boundary conditions of the facia free- 
the same as those of the liner free-body. 
5. The reactions on the subgirts are idealized as point loads « on both ‘the f 
f facia and the liner free-bodies. 
_ 6. The wind loads on the facia and the — are static and uniformly distributed. 


4 beam theory i ignoring the shear deflection, ‘the facia moment | over r the support, 
_ the liner moment over the support, and the midspan deflection can be readily 
found from indeterminate structural The resulting are 


ed as follows. fom wind oa i ven 


The f facia | load factor is 


in which R, = = facia load distribution factor. BP) 7. 
The eer load distribution factor is given by 
‘The bending stiffness ratio (facia to liner) is given by om 


which R, = bending stiffness ratio (facia to liner); E,= = Young’s: modulus 


The nreceding Theoretical and The Conventional linear elastic 
| | 
i 


of facia material; J, = effective n moment of iesatia of facia profile; E, = = Young’ s 
- - modulus of liner meterial: and J, = effective moment of inertia of liner profile 


‘The reaction at the midspan subgirth is given by 


The facia moment over the center support me the conventional sign convention * 


in which M, = facia moment over the center support with the conventional — 


sign convention; and Cy S facia moment coefficient at the comter — 


i The. liner moment over the center support with the oon sign convention 


in which M, a liner moment over er the center support “a the , conventional 7 
« sign convention; and C,,, = liner moment coefficient at the center support 


3 calculated 
ar 


= 


facia deflection 
midepen facia deflection coefficient calculated from 
ox 1927 oars compe 
Substituting Eq. 6 into Eq. 12 and reducing, 


on. 


_ Eq. 13 indicates that the aye deflection is independent of the load Gatciont 


q 
i 
in which P = reaction force at the midspan subgirt, positive when the midspan a 
ra 
| 


In normal applications, the facia stiffness is normally higher than the 


liner stiffness. Therefore, the sensitivity analysis | is conducted for a stiffness 
ratio, R,, of 2-10 in increments of 2. The facia moment coefficient at the — 
center is plotted against the facia load distribution 


Absolute value of facia moment coeffic 


o 


‘Load Distribution Factor (Two-Equal 


nt coefficient, 


lve of liner mome 
o 


Absolute val 


Facia distribution factor, R. 
FIG. 4.—Seneitivity of Liner Moment over Center Support with Respect to Facia_ 
Load Distzibution Factor (Two Equal Spans) 


= 

upport wi espect to Facia . a 

4 
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Fig. liner moment cosfficient at ‘the center support, is plotted 
= 


against the facia load distribution factor, R,, i in Fig. 4. The maximum variations 
— n the fa facia moment from a condition of 100% load on liner (ie., R = 0) - 
to a ‘condition of 100% load on facia (i.e., R, ©” are presented in Y Table - 

: - tk which the relative percentage change of facia moment is calculated with 


_ Fespect to the load distribution of R, = = 0. Similarly, the maximum variations 


of Facia Moment with Respect to to Distribution 


b Facia Moment Coefficient, Cx, 


Relative change 
of facia moment, 
Xx Stiffness ratio, R, R,=O0 | @ percentage 


+18. 


TABLE 2. Moment with Respect to Load ed Distribution 


Liner Moment Coefficient, C,,. moment, 


0.02143 
al 0.01531 


in the liner moment are presented in Table 2. As it can be seen from ‘Tables _ 
1 and = the calculated allowable load according to the working stress design 
principles w will be s significantly affected by the assumption of the load distribution 
between the facia andthe liner, a "01x 


In order to ‘0 illustrate the » impact of the load distribution on on the design of ; 
; ' -adouble skin wall panel system, the following example design analysis i is provided 
a typical two-equal span structure conforming to the conditions shown in 


Fig. 2. The he following s structural parameters are assumed to be known: 
= 0.043 in.*/ft (5. 868 om cm 


Stiffness r 


= = 0.032 cu in. (1. 719 om 


= E,, = 29.5 x 10° psi (203. 55 x 10° 


: in which S = = effective section m modulus of t the facia at he he center support; 
and S, = = effective section modulus of the liner at the center - support. own 


For illustrative purposes the simplified design criteria. will 


The m midspan ¢ deflection shall not exceed 11/180. 
3 2. The maximum facia stress at the center support shall not exceed 26,6 667 
psi (184,002 kN/m?’). 


3. The maximum liner stress at the center support shall ‘not exceed 26,667 


Case 1 Design: 100% Load on Liner (R, = 0). —Using the atecsmentionsd 
‘the following cal calculations are given: 


5 x 10° x 0. 29.5 x 10° x 0.175 


=4.07 . 


3 x (0.4587) 
Cay = 


Solving fora allowable wind load governed criterion no. | (midspan 


5 x 10° x 0.175 


in which W = allowable wind governed d by design criterion 
(midspan deflection) = 47.6 psf (2,280 N/m? 
_ Solving for allowable wind load governed by the — criterion no. 2 (facia 


0.086 x W, x 10? x 12 =0.2x 26 667 .. 


in which W, = allowable wind load governed by the design | criterion no. 2 - 

(facia stress) = 51.7 psf (2,475 N/m’). tt 

__ Solving for allowable wind load governed by design criterion no. 3 (liner 


| 

| 
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DOUBLE SKIN WALL Galt 379 

0. 039 x W, x 10? x 12 = 0. .032 x 26, 2.10% 


which } W, = allowable wind load governed by the criterion no. 
wim 
~ Comparing the preceding three allowable loads, the f final allowable load for 


* ‘the design is 18.2 psf (873 N/m’) governed | by the liner stress§ _ 


Case 2 Design: 100% Load on Facia (R, fi), .—The following calculations 


asia 


25 


in which W, = 42.8 psf (2,050 N/m’). ob 


0.0211 W, x 10? x 12 = 0.032 x 26,667 . 


in which W, = 33.7 psf(1,614N/m*), 
—~ Comparing the preceding three | allowable loads, the final allowable load is 
psf (1,614 N/m’) governed by the liner stress. al al 
_ Comparing the preceding two designs, the allowable load of the case 2 design 
(100% load on facia) is 85% higher than that of the case | design (100% load -_ 
on liner). This large discrepancy between the two designs has ot 
engineering implications in the following manner. If the actual wind load 
distribution is close to the case | design, the | case 2 design would be unsafe. a 
If the actual wind load distribution is" close ‘to the case 2 design, the case 
1 design would be unnecessarily overconservative. Therefore, the basic under- 
a standing of wind load distribution is essential for a sound economical design — 
_ of a double skin wa wall panel system wet! be boot vel 
208 


The design of any wall system 1 must t consider two wind load directions, the 


_ The positive wind load is primarily produced by the momentum change of the | 

7 “moving air mass. For a double skin wall system, the facia provides the first 
barrier to the moving air mass. Thus, it will directly receive the positive wind 
load. However, due to the fact that the facia construction is normally not airtight, 7 
a certain degree of pressure buildup will happen within the space between the 
~ facia and the liner under a persistent wind condition. Therefore, it can be 


= the load acting outwardly that i is ; normally referred to as a 4 negative wind load. 4 


} 
tape 
— 
“t jo t 
)(same as case design). 
| 


be acting on the facia with the balance of the positive wind load acting on 

liner. The negative wind load is primarily produced by the 
pA - gradient across | the ¥ wall. _ For a double skin wall system, , the following three ; 
atmospheric pressure regions can be Region 1—inside th the 


garface. 
negative wind load on the facia is produced the pressure 

difference between region 2 2 and region 3. The negative wind load on the liner x 
is produced by the pressure difference between region | and region 2. In a _ 
no wind condition, the atmospheric pressures in regions - 1, 2, and 3 are the 
‘o At the instant of a pressure drop occurring outside the facia surface, 

- an outward force will start to act on the facia due to the poeneupe difference 
‘between region 2 and region 3. As the pressure in region 3 decreases with 

- time, the pressure in region 2 will also decrease with time due to the air leakage ‘ 

_ through the end lap joints and the sidejoints of the facia | panels. Since region _ 
«3 is an open field, the rate of pressure drop in region 3 would be greater 
than that in region 2. Therefore, the net outward force acting on the facia 
_ continue to increase until a steady minimum pressure in region 3 is reached z 

‘If the minimum pressure in region 3 steadily persists indefinitely, the net outward _ 

force acting on the facia will start to diminish due to the continuous air leakage 

from region 2 to region 3. At the same time, the , outward force acting on the 
liner increases due to the increase of the pressure difference between region 

- Land region 2. If an ideal steady wind condition i is assumed, the total negative a 
wind load will eventually act on the liner alone. However, the ideal steady 

_ wind condition could never exist in nature and a realistic evaluation of the = 
- wind load distribution for design purp purposes becomes a challenging: task for 
Bprediis in the field of wind engineering. The design negative wind load 


example, the maximum negative wind load on the facia would be considered 
for the facia to the subgirt connection design while the maximum negative wind — 
~y on the liner would be considered for the liner stress consideration. : 


= is further complicated by different design considerations. For 


addition, the _ wind load distribution is affected by the wall opening area of | 
‘The preceding examinations pinpointed the general expected trends of wind 
load distribution characteristics. The quantitative wind load distribution for design 
has yet to be developed a wind 
research. 
Conciusions 
Conc.usions 
The following conclusions may be drawn from this study: _ ibs ae: 
‘od . The assumption of wind load distribution between the facia and the liner 
a is found to be very critical to the design of a double skin wall system. ems: 
2. The suitable wind load distribution for design purposes has not been defined — 
‘in 1 the industry, causing 4 a great oye of design confusion and legal — 
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Appenoix —Rerenences 
» 


Strectures, A58.1-1972, American National Standards Institute, Inc., New York, 
2. Cold-Formed § Steel ed., American and Institute, 


"The following symbols are 1 used in this paper: ae as 

reaction coefficient at midspan subgirt; 

Cy = facia moment coefficient at center support; 
& = liner moment coefficient at center support; 
Young’s modulus of facia material; 
= Young’ s modulus of liner material; 
effective moment of inertia of facia panel; 

I, = effective moment of inertia of liner panel; 
ot L= span of wall structure; 
facia moment at the center support; 


= liner moment at the center support; toy Pa 
P = reaction force at the midspan subgirt; THe 


the ratio of facia bending stiffness and liner bending : stiffness; acre 


= line load distribution factor; yr 


effective section modulus of facia at center support; 
‘S, = effective section modulus of liner at center ‘support; oe 


he 


w= total uniform wind load; Sterne ‘are 


- = uniform wind load on facia; 
= uniform windloadonliner; = 
W, = allowable wind load governed by midspan deflection; 
= allowable wind load governed by facia stress; 
W, = = allowable wind load governed by liner stress; 


| 
| 
= 
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OPOLOGY AND OpTIMALITY ¢ OF CERTAIN TRUSSES 


ow 


of topological ‘optimality of “structures | are more complex and 
~ challenging than the currently dominant types of problems concerning optimizing 
- node coordinates and member sizes. From the weight-saving aspect, optimizing 
the topology can also be more rewarding. In discrete structures, topological — 
- variation can be modeled as a mixed integer programming problem, the solution — 
of which could be quite complicated. In the current state- of- the-art, optimal 
4 connectivity is handled by introducing all acceptable ‘members connecting the 
_ nodes of the structures and eliminating them subsequently by means of a linear 
 peegrmauting or other automatic proportioning techniques (2, 6). The topological 
nature and its relation to optimality are very rarely analyzed. 
= this paper, some known results of Maxwell’ 's (4) and Michell’ s (5) theorems 
are briefly reviewed. A special characteristic o of a certain class of common 
trusses is then presented as an extension. The condition n necessary for trusses 
to have this characteristic is then derived, and its implications on topology, 
optimality, and stress-ratio design are examined. Simple and more complicated 7 
 Bacxcrounn 


SS As early as 1904, he dealt with problems of total optimal layout of structures 
and provided sufficient conditions for optimality, which are still useful today, 


- elated to this paper, will be examined herein. Refs. 1, 3, and Ta are yan 
of modern statements and interpretations of Michell’s theorems. = ~_| 
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SS Michell’s trusses are dunaeleal by a given set of forces, P, in equilibrium, © 
acting on a truss at given coordinates, Q. For such trusses, Maxwell (4) showed 


= 20,2, = =C 


7c 


q “which = the lengths and ‘the in n the ‘member of ‘the 
truss; and C = a quantity nieseaae of the configuration or the origin of © 


Starting from this result, and allowing the stresses in all members to reach 


an allowable constant stress, o,, in tension, and in ‘compression, Michell 


in 1 which Var ond Ra the material volumes of tension and compression members, 
respectively, such that the total + By 
V=—(c a, +9, Me, V, + 0. 


and recognizing { from Eqs. | and 2 

in which the tilde denotes transposition; and a an = constants 


From Eq. 5, ‘Michell concluded that the truss that has the ke least m ties 
vl It should | also be noted that P, by definition, includes is ; applied loads, 
Po, a as well as the reactions, P,. If Q is partitioned into ordinary coordinates, 2. 
which may be P,, and support coordinates Q,, then Eq. | may 
=0,P,+0,P.=L 


For a structure with members exclusively in tension or in compression, 

“the equality Z = +C follows from Eqs. 2 and 4; thus for constants, P and = 

the volume (in Eq. 5) is constant regardless sal the configuration and the 


location of the unloaded nodes. 


| 
4 


mined from statics), P is constant and by Eq. 7, C is constant regardless of 
the topology and the degree of indeterminacy, provided that the coordinates, — 

Q, remain fixed. Statically | trusses are, of course, a subset of this 
3. The volume of trusses with constant, C (such as determinate or only 
externally determinate trusses, as described in the aforementioned), which are 4 
_ also fully-stressed, can be minimized by either minimizing the volume of tension 4 
members, V,, or the volume of compression members, V. , Since these are — 


_ In addition to the types of structures with special features, such as those — 


| stir in the preceding sections based on Michell’s theory, there exists. a ‘class 


of trusses for which the sum, S, as as defined in the following, remains constant, * 
regardless of member proportions, provided that the geometry of the truss remains — 


in “ia which Zand C are defined in Eqs. 4: and 7, respectively. 
, a Before examining the conditions under \ which Eq. 8 is satisfied, its general 4 


be fully stressed, he sum, Ss, remains s constant in 1 such 4 a truss, even "when 


it is made indeterminate by adding redundants. 


¢. aenpening |C| = |2 L,F,| as an objective is consistent with minimizing 


a 3. Consider the stress-ratio » design | method, in which th the cross ‘section, Ay 


7 which F, = the force in member i; o, = either oC, or o.; and the superscript — 


a: = the iteration or the cycle number. This method is known to i 
toam minimum weight design for for a "single loading (8). If Eq. 8 hold 8 holds, then itt 
. follows: from the aforementioned that S remains constant in n every iteration. 

_ Since Z monotonically decreases (8), the quantity |C| increases until it reaches | 
5 ‘the optimal, statically determinate configuration. In that case, of course, both — 
_ Cand Z become the constants C, and respectively, 

4. it rd —_ the definition of C and Z that if Eq. 8 is satisfied, ‘then 


= — = constant; for. all tension members if C > . 


| 
4 


Fe 


F.=—= constant; members ifC < 


of the forces in tension or 
compression, and L, F, and L. F. denote the sum & F, L, for 
members of tension or compression, respectively. Therefore, if the member 


a is designed using the simple safety condition A , = |F,|/o,, a major implication _ 


of Eq. 8 is that throughout the process of i improving the design, e.g., by using 


, the stress-ratio method, the ee (or V.), remains constant, while Ve 


since in fixed configurations: the sign of C i is determined by 


of material for tension members or for compression _ members, , depending upon 
: whether the structure in question is characterized by a tension or compression 
pth. 
i. On comparing this class of structures exclusively in tension or compression, 
_ for which the required volume of material i is constant, as examined in the preceding ~ 
section, this case may be regarded as an extension of the notions of Michell’s 


theory to general trusses. The for hold is derived 
inthe next section, 


Two o of the laws that govern truss are: Equilibrium among the 


applied loads, P,, and the member forces, F; and (2) Compatibility between — 
the node displacements, . 5, and the member deformations, A. 4. These ¢ requirements ie 


may | be expressed respectively as: ‘oft yt 


and N 


—— The matrix, N, incorporates the geometry of the truss as well as the member 
incidences. Its rows correspond to members and its columns to nodal degrees 
of freedom. If augmented by the matrix, N,, compatiity with ©., the 
- coordinates may be related directly to the member lengthss ts 
L= [NIN,] (13) 
The set of members may be partitioned into a basic determinate subset, 7, § 
_ and a complementary redundant subset, R. The vectors, F and A, and the matrix, © 


| F, 


© 


| or I . (105) 
q 
| 
1 
| 


uantities ey and Ss ” be computed from Eqs. 4, 4, 7, and ene 


in 1 which sgn on the sign of C; and J = a 


with +1 along the diagonal depending on the signs of F. 
In what follows the sets, T and R, will be assumed to be optimal Saale 


will be used to denote the forces F, when F, = 0. Thus =e 


on The and S, corresponding to F, by Eq. 16 


4 
the su it can an be shown that 


= S, + Fe{[(sgn C,) 1+ Jp] Lp —N,Nz'[(Gen C,)1+J,] L;}. 
it 
the forces, F,, ce carry the same signs as F.; the a 
E: toduadents does not change the sign of F,. This can be shown to ral true 


under a single loading condition, provided that the basic structure, i.e., the 
4 set T, is selected to be optimal. The proof hinges on the monotonic convergence 
a of the stress-ratio design process (8) during which the cross- sectional | areas 
of a set of redundant bars, R, decrease continuously, defining a complementary, 


hy + L;= = a 


contains of tension members and I L 


4g 


r 4 
| 


a 
of compression members s only; 2 Zeros replace the lengths of compression members - ; 
in L* and the lengths of tension members in so that L=L* 
~ 3 In order that S remains constant as required by Eq. 8, the second term in 
Ls - the right-hand side of Eq. 20 must vanish with nonzero F,. - Therefore, heed 


= 


= loading determines which 0 one | satisfied. 


INTERPRETATION 

2 23 ar are evidently conditions ¢ on m the ery of the given t truss 


_ The deformation of the ae must conform with this displacement for 


_ Compare now Eqs. 2B and 26. ». Eq. 23a, € e.g., , may be i i 

that the set of tension members in the truss can all be stretched Cgc 
to their lengths, while the compression members conform with the resulting 
nodal displacements without any change in length. A similar interpretation af 
‘ Eq. 235 can be made by interchanging roles of tension and compression. a 
: interpretation can be seen to be aeeis with the implications of wad 8, 


introduced in the preceding secti oe, 


Exampues 
illustrate these concepts, “three are presented, emphasizing 
Example 1.—Fig. 1 shows a five-bar truss. The data and calculations for 
this truss are shown in Fig. 2. Bars 1-4 constitute the optimal, statically determinate 
truss, for which the forces, F,, and the qualities Z,, C,, and S, are calculated. * 
Consider now the addition of bar 5 as a redundant, shown with a dashed line 


Fig. 1. It can be ‘shown that both Eqs. are satisfied. This is ‘the: 


carry the sar same fe sign. "aes all member cross ‘sections to be equal, the forces, 
F, Z, and C are recalculated. Maxwell’s Theorem (see Eq. 1), can be easily 
- verified, and S§ is shown to be equal to S,. If now a different determinate 
7 set if formed by deleting, e.g. » bar 3 (the he rows of N must b be ‘Tepartitioned 


: t in cases in which the redundants 
in Eqs. 23 are mutually exclusive. 
&,, Of a determimate set of members: 
= 
@ 


accordingly) none of Eqs. 23 are satisfied. The test fails i 


Eg. 


1/72 
eto. 


Pte 


ae 


L, = 
s=z+ =12=s ~ 4 


FiG. 2.—Data and Calculations for Five-Bar Truss 

_ used with different allowable ‘stresses in tension and in compression. Eq. 9 © 

is: used in the iterative . design. The calculations | and data are shown in ‘Fig. 


4. The expression for the volume in Eq . 5 is ‘shown to hold, since in every x 
iteration the cross section is redesigned to be fully- stenned - The results of “ 


this case because 
_ Example 2.—A three-bar truss is nstrate that 
| remains constant eve method is 
| | | 


sh 


ai 


[en 
[ara oan | 


4.—Data and for ares: Bar’ Truss 
a 


be 

FIG. 5.—Topological Implications of Eq. 23: (a) Optimal Determinate Truss; (b) 


az + BC -27/118V2, 
| 


La 


, th the first two cycles are shown i in n Fig. 4. If ow process is continued, the cross i 

section area of bar 3 will continue to decrease. Bars 1 and 2 constitute the ps 

_ optimal determinate truss for which S, = 6 as for every cycle of the stress-ratio a 
Example 3.—The topological implications of Eq. 23 are shown in Fig. 5; 
a. data for this truss are shown in Fig. 2. The optimal determinate.truss includes 


yorvadhed an sori) We prom 


coroor 


759.7 


FIG. 6. Two Cycles of Iterative Deel n of Eight-Bar Truss 


bars 1-6. When two are added, 7 and 8, they are stressed 
in tension and compression, respectively. In this example, Eq. 23 is satisfied. 
Fig 5(5) shows the compression members to be all shortened to three quarters x 
of their lengths, and the tension members (double lines) all fi it with no —— 


a > 
— 
‘on 
— 


1981 
porven compression are eiemiaeiel and that Eq. 23a does not hold in this case. 
Finally, the forces resulting in the first two cycles of ‘@e if rative design a 
are shown in Fig. 6 to demonstrate that FL L, i.e., the sum of the products 
of compression forces by their respective lengths, remains constant while FE _ 
LL, the tension counterpart, keeps decreasing. It can also be shown that if any 
two redundants other than bars 7 and 8, are removed, then no part of Eq. 
7 
‘In belied trusses, it is more difficult to trace the behavior described. However 


in is relatively easy to check whether Eq. 8 is satisfied 


q 


| 
7 


AG. Trusses: (aT Truss Eq. 
_ For ss the truss in Fig. 7(a) satisfies the condition exactly, while the | 
one in Fig. 1b) does so only approximately. In the second » case, the sum, 
S, approaches a constant very rapidly if the stress ratio design p process is used. 
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Aprenoix il. — Notation 


= identity matrix; 


= diagonal sign matrix; 


Pp = external forces on ‘truss; ker 
P, = reactions, subset of P; 
applied loads, subset of P; 


subset of Q, compatible with Po 
volume of truss members; 
_member deformations; and 
"joint 
Subscripts 


number of truss members; 
= subset of redundant members; 
T= determinate subset of members; 


0 = optimal determinate truss. 


following symbols are used in this paper: 
A = member cross-sectional area; 


q | 
| 
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STRUCTURAL DIVISION 
RE 
_ By Awadh B. Agrawal,' Leslie G. Jaeger,” and Aftab A. Mufti,” A. M. ASCE b 
| 3 In recent years many investigators have dealt with the response analysis of 
: O structures subjected to strong earthquake motions. Scaled models of reinforced 
concrete frames (21), shear walls (14,17,24 25), and shear wall- frame ‘systems: 


a ( il, 12, 16,26) studies, analytical ‘models based « on “modified EI” procedure we were 
developed to compute the inelastic response of the structures under earthquake ~ 
loads. Takayanagi and Schnobrich (24) utilized a variation of the ‘ “modified _ 
- frame” method to simulate the dynamic behavior of coupled shear walls. 

‘‘modified procedure and its variations are well suited to the analysis 
of plane and framed shear walls if one. “gets a pure cantilever deflection for 

_ the structure. However, given a sufficiently stiff coupling by beams or rs a 
or both, a system of horizontal shears is set up and deflections become smaller 
_ with increasing amounts of backbend (15) being induced into them. Consequently, 

- treating the planar reinforced concrete structures such as shear walls to be 
_ina state of plane stress will permit a more realistic assessment of their structural 
“behavior under "earthquake induced ground motions. A method based on finite 

element plane stress analysis will, therefore, be suitable to investigate general : 


a finite element method has to be computationally efficient and reliable. This 
can be achieved by using efficient solution algorithms, experimental tests, and 
reasonable engineering criteria as applied to reinforced concrete. 


sf wall systems. It is to be noted that a computer program based upon 


4 stress finite element approach in determining the nonlinear response of reinforced 
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a ‘concrete structures subject to earthquake motions. s. The finite element formulation 

_ for conducting the nonlinear dynamic analysis requires two types of mathematical — 7 
modeling: (1) A material model to represent the inelastic behavior of reinforced __ 


a under reversible loads; and (2) an analytical model to trace the oe 


force-displacement characteristics of the structure subject to ground excitations. 
The material model (1) treats an element of reinforced concrete in biaxial — 
stress states, and provides for the cracking and crushing , of concrete and yielding = 
of steel reinforcement. The formation of a second set of cracks under — 
reversals in addition to the existing ones, as well as the reopening and closing _ 
of previously formed cracks, are also incorporated in the model. The derivation 
of the material constitutive relationships is based on the idealization of the — 
‘stress- -strain curves of steel reinforcement and concrete in compression. eer 
The analytical model for conducting the nonlinear dynamic analysis is based — “4 
‘on a lumped-mass approach, and utilizes the Newmark’s “B- method’’ to perform 
_ step-by-step integration of the equations of motion. At the end of each 
_ time interval the model checks for the cracking, yielding and crushing of concrete, 
opening and closing of the previously formed cracks, and yielding of steel 
reinforcement. These effects, resulting in nonlinear response of the structure, : 
vedio bag aT anotiom bows: mood (21) 


q 


FIG. 4 —Uniexial Stress- Strain Cyclic Loading 

are incorporated in the a analysis procedure by mo modifying the elasticity n matrices 

of the affected elements, and computing the release or corrective loads (28). 
The next time interval considers these loads and the incremental — 

accelerations, to continue the analysis, RAS tg 

_ The finite element formulation is employed to trace the response of a — : 
-storyed shear wall tested by Hsu (14) under simulated El Centro earthquake 

of May 1940. The usefulness of the material and analytical models is substantiated 

by comparing the results of the analysis | with the experimental data. bili a 


a 
= 


SHEARWALL 


as a — elastic and brittle material in tension, and elastoplastic in uniaxial | 
compression. A simple bilinear stress-strain curve, reflecting the ctusithediening, 
_ is employed for steel reinforcement. The hysteretic response curves for the 
constituent materials are presented in Fig. 1. 
5 F _ Plain concrete is assumed to be a homogeneous and isotropic material prior 
to its cracking or yielding. In a state of ’ plane stress i= elasticity or material 
property 1 matrix is written as: lwo we 28; 
in which E. and v = the modulus of elasticity and Poisson’s ratio of ae 7 


laine! 
dos 
to 


* Cracking and Crack Closing.- —The effects of cracking and crack closing in 
1 the material model are considered according to the present state-of-the-art 
_ (7,9,13,23,27). The concrete has very limited capacity to resist tension, and 
 istherefore allowed to crack if and when a principal stress exceeds the permissible . 
Ee stress of LA set of cracks i is formed in an element normal to the principal i 
iA stress direction as shown in Fig. 2. The cracking is modeled by reducing the 
modulus of elasticity E. to zero in the principal stress direction. The material | 


i 
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the transformation matrix diane the material coordinate system 


xy -y to the element coordinate system X-Y; and ‘‘a’’ = the shear transfer facto ¥ 

account for the aggregate interlock and action on rough surfaces 
opened cracks. Although this factor decreases with the opening of the crack, ad 
and increases with the relative displacements parallel to the crack, an approximate 
constant value of a = 0.5 (for normal weight concretes), seems to be sufficient 
Al-Mahaidi and Nilson (6) have the effect of crack 
defining the parameter as: Gow try 4 a: 


| 


as 


of the cracking, This is based on the assumption 
€, = 5e,, at the time of the failure of the structure subject to monotonic loads 
and moreover does not take into account the increase in the shear capacity 
of the cracked concrete due to the shear-slip along the crack faces. In the 
‘present investigation, therefore, ‘‘a’’ is assumed to be a ‘constant quantity at 
In addition to the single set of open cracks considered in the foregoing, the oa 
_material model provides for five other crack configurations. The six possible 
crack configurations (9, 27) for an element of reinforced concrete are: (1) No 
cracks; (2) one set of open cracks; (3) one set of closed cracks; (4) one a 
of open and one set of closed cracks; (5) two sets of closed cracks; and (6) 
two sets of open cracks. No distinction is made between the material properties — 
of uncracked concrete and concrete which has cracked and has all of its =— 
closed. Similarly, the material properties remain the same for concrete which — 
has just one set of open cracks, and the one in which one set of cracks | are 
closed and the other set of cracks are open. 
_A set of open cracks, say located in the x direction (see Fig. 2), will close 


_ parallel to eracks) exceed the strains &, acting normal to cracks. In other words 
these cracks, crack closing will occur if or the crack width 


= €x sin? cos’ - 


an 


It It should be noted that the term “‘crack width” 

since it is expressed in terms of strains. eal — oad 
=. A set of closed cracks reopen as and when a hl acting norr 


cracks become tensile. The concrete having open cracks will crack once m 
giving rise toa second set of cracks, if o{ > of . The new set of cracks i 


__in which [7] . 
| 
q 
- J 
nal 
in a Ggirection which is normal to the existing set (/). OWever, in d 


SHEAR WALL 


the cracking of uncracked concrete is to determine the direction 
_ the new set of cracks. The concrete having two sets of open cracks is assumed — 
to be capable of transferring shear only, and a is reduced to one half for 
a two sets of cracks in Eq. 2. Since cracks are formed at different _ 
bp angles, the property matrix for each set of cracks is transformed separately — 
__ to the element coordinate system. The two are then added up to give the material 
i. "property matrix [D.] ., for an element of concrete having two sets of open 
Yielding Crushing of Concrete. —The von Mises yield criterion, shown 
in Fig. 3, is utilized to check the yielding and plasticity of concrete in biaxial _ 
The concrete starts yielding if > 0, and F(o) is given 
in which of = the uniaxial stress of concrete (see ‘Fig. — “initial 


ote sn. 3.—Fa Envelopes for Concrete 
the yielding of concrete in ‘the reggae procedure. The concrete is assumed — 
have been crushed if strain a a value of 0.003. 
Itmay that an alternate approach based on ‘the concept of ‘ equivalent 
<9) uniaxial strain’”” has been proposed by Darwin and Pecknold (9,10) to model ° 
the behavior of concrete in plane stress states. By ren of this concept, biaxial: q . 


model for concrete is obtained under biaxial compression. The use 
of von Mises criteria has advantage because the elastoplastic theory could be 

applied for the computational purposes. Keeping in mind that the order of accuracy 
f in both approaches is comparable (3), and the cost of computation at present 


is therefore, idealization for concrete in in compression and 


i 
| 
| 


: ‘the use of von Mises yield criterion under biaxial compression, in the opinion — 


: material, and smeared or ‘uniformly laid over an n element of concrete (9,13,23,27). a 
For the reinforcing bars inclined at an angle y to the X-axis of the element, a 
the material property matrix [D,] is written as: 
= 0 a (9) 
in which E, = the ‘heii of caatieny of steel; and p p= the reinforcing hy , 
Composite Material Property.— —The composite material property matrix [D, 
or an element of reinforced concrete is obtained by adding the material property | 
- matrices [D.] and [D,], respectively. A null matrix is used to replace [D, a 
> the concrete is crushed; and E,,, the modulus of steel in the strain hardening e 


range, seplaces E, in matrix [D,] ] once enters inelastic range. 


Move. FOR ANALysis 

Finite’ Element. plane stress (22) with three degrees-of- 
- freedom at each node is employed to discretize the structure for nonlinear © 
reinforced concrete analysis. The degrees-of-freedom are the translations u and 
oy in the X and Y directions respectively, and the inplane rotation 6, = dv/ax 


of the node. The presence of the inplane rotational degree- -of-freedom permit 
- the idealization as line elements of coupling beams and frame members 2) i 
in coupled shear walls (5,18) and shear wall-frame systems. Ce 
_ Formulation and Solution of Equations of Motion.—For the structural system, 
the formulation of the equations of motion is based on the concept of lumping 
the masses at appropriate nodes. The finite elements, therefore, utilized fund 
- the system are considered as “massless: elements. The degrees- -of- 


freedom associated with zero-mass terms, i.e., dependent degrees-of-freedom, 
- are condensed out (1) before performing the dynamic analysis. The motion 
: of the simplified model in terms of the independent degrees-of-freedom is governed © 
by differential equation of the type: » > om betas: 
(MCF) + + = 
in which [M] = , diagonal ‘mass matrix; [C] = damping matrix; {r,} = vector — 
of base accelerations; and {r}, {7}, and {7} = vectors of relative displacements, 
velocities and accelerations of the masses, The stiffness matrix. 


if 


stiffness matrix “TKI. which eonneapent to the dependent and independent 
_ degrees-of-freedom, respectively. A closed form solution of Eq. 11 is not possible 
because of: (1) The need of changing or updating the stiffness matrices [K] 
& and [Ke ] to reflect the inelastic behavior of the structure as the loading | - 


| 
| 
| 
| 
i 
1 
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WALL 


#1 
in which {A7,} = vector of base and {Ar}, 
and {47} = the vectors of incremental relative displacements, velocities, and — 


accelerations of th the 1 masses, 22 


time increment or step At. The value B = 0.25 corresponds to the ding sa 
thai the acceleration remains constant within a time step, and an average value 
of the accelerations between the nth and (n + 1)th time steps is used to evaluate e 
the required quantities for (n + 1)th step. 
_ The solution of the set of second order differential equations representing 
the motion of the simplified model, Eq. 13, is reduced to a sequence of ieieiien ; 
matrix operations by assuming that between two time references the material 
remain unchanged f for the purposes of “computing the incremental 
displacements in the structure. This assumption is identical to the one utilized 
performing the “iterative incremental nonlinear analysis’ of reinforced 
concrete structures ‘subject to monotonically increasing and reversed cyclic 
fe loading (6,9, 13,23,27) where within each iteration the structural system is idealized | 
_ to behave in a linear elastic manner. Once the incremental displacement vector 


nd. is known, the incremental displacement vector {Ar} associated with the 


Damping Matrix. _The damping matrix [C] is usually assumed proportional | 


to either the mass matrix 0 or r the aeiness matrix , or a li linear combination of 


= the constants to be determined based on the natural - 


In this investigation, “the damping matrix: [C] is assumed "proportional on 


— 


‘The nonlinear is traced by solving the 


and distributing the unsupported stresses to the adjacent elements 


__ Taking the steel reinforcement to be elastic and concrete to be uncracked _ 


Summing the element stiffness matrices in an appropriate | ‘manner, 
A the structural stiffness matrix [K] for the complete structure is assembled. 4 
ee Loads are — at the appropriate nodes to form the mass matrix [M a, 


r progresses; and (2) the ran ke induced base accelera- _ tq 
j 
| 
1‘ 
the preselected damping factors for th ure, 
| 
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and then the is mee the matrix ix [K* The fir first i increment 
s base acceleration i is applied tc to the structure to initiate e the analysis procedure. 


a a check i is s made for the state of materials fo or r individual elements. Whenever _ 
q the material constitutive relationships are violated, the element’s material property 
and stiffness matrices are regenerated or updated. Examples of such — 
7 _are the cracking and crushing of concrete, opening and closing of the previously - 
sac cracks, and yielding of steel reinforcement. The plastic behavior a 


concrete does not fall into this category (29). The matrix updating is followed 


by adjustment of stresses in the affected elements, and calculation of the release 
‘ or unsupported stresses. These are then converted into release or corrective - 
loads (28, 29), and are applied to the structure as a set of ‘‘static loads’’ (24). 


Wall thickness = 1" 
Reinforcement: 
No. 
es 1" c/c 
(44) Horizontal No.8 
wires 2" c/c 


_6x3."=18. 
& 


al 
AG. 4— Shear Wall D-4(14) (lin. = 254mm) 
‘Their effect, bakit ' is consi time interval. This means thet 
a for each time interval, except the first, incremental nodal displacements are 
_ produced due to two effects: (1) The increment of input base accelerations; * 
, ¢ (2) the release loads originated in the preceeding time interval. Lastly, 

: if required (because of the updating of the material property and stiffness 
matrices), the structural stiffness [K] is reassembled, and the analysis procedure | 
continued for the remaining part of the earthquake record. 
a In this investigation, the stiffness matrix [K] is assembled as an upper triangular — 

banded matrix. The decomposition procedure used for obtaining the matrix [K*] al 

from the matrix [K] is also utilized to determine the incremental nodal 
- displacements due to corrective loads. These and other details concerning the 

a of release stresses ane corrective loads 2 are given elsewhere (1,4). ia 


1 


— 
> 
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The usefulness of the proposed material mode! and the analysis procedure ; 
_ in simulating the required behavior is demonstrated by their application to a 


; 4 three storyed shear wall tested by Hsu (14). In the experimental test, the specimens Z | 
* were subject to simulated N-S component of the El Centro earthquake of May - 
4 18, 1940. The time axis of the original earthquake record was compressed by 


wil 


‘oes 1.—Material Data for Shear Wall D-4 (1 4) 


win 


| Shear Wall D-4 Model Ast and Mode Hote 
bens FIG. ified Model for a ove (M1) bee 
_ The geometry, cross section, and the finite element idealization of the selected | ‘ 


specimen, designated as D-4 in the experimental test (14), are presented in _ 
_ Fig. 4. The pertinent material data for the shear wall is given in Table 1, and 
includes modification ir into bilinear form of the stress: stress-str: strain be of 


q and were inc! walls was 
| par | per |» | per | per | 
Test square square | square square square | Verti- | = 
4 {3.8 x 10°| 4,720 | 409 |29x 7.03 x 10¢ | 53,500 | 53,500 


wire reinforcement as given in 1 (14). The wall is discretized into 120 rectangular be 
*; elements. The finite element grid i is varied between ‘storys t to o improve the ¢ estimate 


__ The simplified model for the purposes of dynamic analysis is shown in Fig. 
5. Also marked in this figure are the nenneal frequencies and mode shapes — 


FIG. Sibi ‘Base Shear and Base Moments for Shear Wall D-4 4a diel 
of the uncracked elastic system. The | computed value of the lowest i 
_ frequency of the shear wall is 11.84 cps against a measured value of 8.5 —_ 
in the test. For a similar shear wall, specimen D-3 not modeled analytically 
_ here, a natural frequency of 10.5 cps was obtained in the test (14). Otani (19, 2) 
: and Hsu (14) have attributed this difference between the measured and competed 
frequencies to the existence of the fine cracks in the test specimens — 
are developed during its handling, and due to shrinkage of the concrete. Se woah “| P 
-.. The base accelerations of the El Centro earthquake record, with maximum 
_ amplitude of 1.05 g and time axis being compressed by a factor of five, were 


| 

| 


ames at a constant interval, AT, of 0. 002 sec, ‘and were applied to ‘me 
_ Structure in the X direction. However, for numerical integration, a time step : 
__ At, of 0.0005 sec was employed so as to ensure the stability of the — 


“process. The nonlinear response of the shear wall has been traced here for 
_ the first 1.5 sec of the earthquake. This is the duration within which aot 
‘Tesponse parameters attain their respective maximum values during the first 


“a 


AIG. 7 —toanies Relative Displacements for Shear Wall D-4 (1 ‘in. = = 25. 4 mm) : 
enn run (14). This approach of analyzing the structure only for the most intense 
portion of the earthquake record instead of its entire duration has also been | 
adopted by other investigators (24), wW 
- The analysis has been performed v with three different combinations of ol 
shear transfer factor a and the damping coefficient B, of the first mode. The 


combinations are: ou 0.5 and By = 0.0 02; Qa = 0.1 and 


the 


| The general observations about the wave forms are: 
.- The base moment wave form, shown in Fig. 6, is exclusively governed © Eo 

the first mode component, and is very smooth. 


lirst mode con 


3rd level accelerations 


FIG. 8. —Computed Response Accelerations for Shear Wall st D- 4 
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2. Thee displacement w wave ve forms at are very / smooth; at all levels, and are dominated 
f by the first mode component (see Fig. 7). abe 
_ 3. The acceleration wave forms presented in Fig. 8 contain more of higher 
- components, especially at the first and second levels. As the loading © 
"progresses, higher mode components become more perceptible even at the third 
level wave form, in which first mode components have prevailed previously. _ o. 


mbinations were 
4 
Me Nase wave [orm as ine loading 


the | test. . The nye values of the Tesponse parameters of this analysis are are 


Magnitude of 
be 
_ Maximum base acceler- 
Duration of earthquake 


Response accelerations: 


Second level 


First level 


= 


displacements: 
(Double amplitude) 
‘Thirdlevel 
Second level 
First level 


= base shear 


Response base 


Natural frequency 
before test run 
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frequency during 


test/analysis | 


vate 


a 
Measured 


in. 
kip-in. 


DPB 


0.76 in. | 


Test 


ao 


B, = 0.015 


0.80 in. 
0.48 in. 
0.19 in, 
«$38 kip 
183 kip-in. 
84 cps 
2.32cps| 2.48 cps 


aires 1 in. = 25.4 mm; | kip = 4.45 kN; 1 kip-in. = 

ab 


| yielding of steel ante 
reinforcement 
y 
vv A. With: 


Factor 
Damping 


wel Shear Factor a = 0.1 
Damping 0.015 of 


j e response _wave forms presented here are similar to those reported in 
] ent § 

Calculated |_a=01 az01 | a=05 
B, = 0.02 8 =002 | 6,=002 
-  &§ 

10g | 105g g ¥ (105g «10658 
2.03 2.28g | 228g | 245g 
173g | 14g | 139g | 135g | 140g 
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a= - 0. and B, = 0. 015, are » to weer values 
_ than those computed by Hsu (14) using the ‘‘modified EI’’ approach. The proposed — 
- analysis, however, over-estimates response base shear and under-estimates first 
level response accelerations by wide margins as compared to the test. * 
The analytical results for the shear wall | obtained with shear transfer factor — 
- 0.1 and B, = = 0.02 do not differ significantly from those obtained with — oa 
= 0.5 and B, = 0. 02. This implies that in nonlinear finite element —_ 
of reinforced concrete structures the actual value of the parameter a is not 
as important as the fact that the cracked concrete has been allowed to retain 
_ The analytical and experimental crack patterns for the shear wall are compared = 
in Fig. 9. It should be noted that the experimental - crack patterns have been 
“4 i reported in the test only after the third test run was completed. However, 
_ a good match is observed between the two, both in the location and the direction 
of the cracks. The analysis indicate that at all levels of the shear wall, the 
first cracks to form appeared to be the flexural cracks. The cracks then propagated 
_ from lower to higher levels, and from outer sides towards the inside of the — 
shear wall. Due to reversals of the loading, the second set o of cracks +. ue 
formed in many elements before the first set of cracks could get a chance — 
to close. This effect is more pronounced in the crack patterns presented in 
Fig. 9 where the analysis was performed with a = 0.5. The inelastic behavior 
of the steel reinforcement is confined primarily to > the first level elements afi 


principal objective of this investigation was to the use of the 
_ plane stress finite element method in the nonlinear analysis of reinforced concrete 
Structures to earthquake loads. 
— On the basis of the results presented herein, it is shown that the proposed _ 
‘material and analytical models, as well as the analysis procedure, are capable 2 
_ of tracing the response of planar reinforced concrete structures under simulated 
or earthquake induced ground motions. As compared | to the “modified EI”’ 
approach, the results given here are better in terms of response displacements, 
more or less the same for response accelerations at higher levels and base» 
- Rone and not so good in terms of base shear. The proposed model could — 
be improved further by providing a crack-width dependent (6) value of the | 
“shear transfer factor to account for the aggregate interlock and dowel actions 


across open cracks in concrete. This, however, needs experimental 


under cyclic and dynamic loadin 
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Cis constants used to define damping 
damping matrix; nt 


folowing are used in paper: 


aterial property matrix for concrete; “e as 
material property matrix for cracked concrete; : 
material property matrix for reinforced concrete; oo 

material property matrix for steel reinforcement; 
modulus of elasticity for concrete; 


4 © 
modulus for steel in strain hardening range; 


o) = biaxial stress field for concrete; 
ads = acceleration due to gravity; 
(Kul, (Ka ‘partitions of structural stiffness matrix [K]; 

 [K*] = structure stiffness matrix used in dynamic 
= 
= ground accelerations vector; 
= grour 1 v a 
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vectors sof relative dis velocities a and 


freedom; Yer LRN NAL THE 


“transf ormation matrix; 


material coordinate system with with -axis aligned 


ORSON A 8, aa coefficient for first mode in terms of 


ar 


= time interval; 


= time increment or step; 
{Ar} (Ar) (ar) = vectors of incremental relative 
velocities and accelerations of de- 
vector of incremental displacements of of depen- 
incremental ground accelerations vector; 


equivalent uniaxial strain for crushing of 


normal and shear strain components; 


= strains normal to open ¢ crack; 


v wing aa 
¥ 


strains normal to. open ‘crack due to Poisson’s 
effect 


th 
reinforcin an le w.r.t. X- had ‘ 


Poisson’s ratio for concrete; 
= angle of crack in concrete w.r.t. X- 
tensile stress for concrete; Opes DY. 
ove 4a stresses acting parallel to cracks; 
normal and shear stress components ts for con- 
= uniaxial yield stress of concrete; comin 
yield stress for steel reinforcement; and 


= natural frequencies. 
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_ The problem of the instability of a steel ‘stanchion subjected to m to moments — 
about oth its bending axes has attracted the attention of many investigators 
in the recent past. Birnstiel and Michalos (1) conducted ultimate load tests 
on biaxially loading H-columns. A rigorous method of analysis of such pT 
was proposed by Harstead, Birnstiel, and Leu (2). However, this analysis was 
restricted to columns bent symmetrically about the midheight. Furthermore, 
the axial thrust was applied at fixed eccentricities, SO that end 1 moments were — 
developed in proportion to the applied thrust. Milner and Gent (4) ‘studied the | 
problem of columns with end restraints and with initially applied moments, 
but the analysis was again restricted to columns bent symmetrically about the 
midheight. Another method, with the same limitation, was developed by Santath- _ 
adaporn and Chen (5), the analysis being based on a tangent stiffness matri 
_ method. The effects of material yielding, residual stresses, end warping r restraint, 
end bending restraint, and initial imperfection were included. A more general 
-method in terms of nonlinear material properties and unequal end moments © 
was developed by the first writer (8), primarily for steel-concrete composite _ 
7 columns. All twisting and warping effects were ignored, but the analysis was ad 
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The solution was based on the Newton-Raphson iteration technique, and propor- 
tional or nonproportional loading paths were included. Another method applicable 
to steel H-sections with generalized end conditions of loading and flexural = 


_'Lect., Dept. of Civ. Engrg., The City Univ., oat, Square, London EC cv 4 


extended to similar columns with unequal biaxial flexural end ‘restraints (9). 
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6 and included the effects of twisting and warping of the cross section. The _ 
method was based on the use of successive over-relaxation for the solution = 

of nonlinear equations. A very much similar method was given by Vinnakota 


and Aoshima (7) with the limitation that elastic, , perfectly 

It would thus appear that rigorous methods are now to 

: _ calculate the ultimate loads of columns with generalized end loading and flexural 

restraints, i.e., with biaxial loading and biaxial end restraints, unequal at the 
=~ ends. The “ultimate load” could be either the axial thrust — or any one 


stress-strain characteristics, initial imperfections, and residual edi patterns = 
_ are allowable. The method developed by the second writer (6) included the ia 
_ effects of twisting and warping when applied to steel H-sections or I-sections. 
_ The method developed by the first writer (9) on the other hand ignores torsional © 
5 effects altogether, as it was developed primarily for torsionally stiff sections, 
in particular composite columns. . Individually, both methods have been verified 
— ageinet experimental results on columns of their r representative types. However, 
4 > no comparisons have so far been made for ultimate loads obtained by ignoring | 
torsional effects against those obtained by including them. Clearly, the procedure — 
_ ignoring the torsional effects will give speedier results. At the same time, ignoring 
; torsion could in some cases lead to nonconservative designs. Moreover, in view a 
the new design methods being developed at "present [see, e.g., the work 
of Wood (10)] it would also be of interest to investigate any difference in — 
failure modes predicted with or without the consideration of torsional effects. ie 
This paper attempts to identify cases where it would be imperative to consider 
_ torsional effects in estimating the failure loads or failure modes, or both, =—=—> 


Failure Loads Including Torsional | Effects .—The deflected culdiguaiive of an 
H-column under generalized end loading is shown in Fig. 1. The forces ig 
_attheends including restraining moments are shown in Fig. 2 where, for simplicity, 
= the minor axis restraints and deflections are shown. The deflected shape 
of the column in equilibrium with an applied ert can be obtained by the 


jon is: a 


in which . M, and M, = = the external m major or and minor axis moments; b= 

the angle of twist; x. = the major axis flexural stiffness calculated including — 
the effects of material nonlinearity; the total deflection; and = the 


— 
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in which M,, and M,, . = the moments about ‘the x-axis at the bottom and 
top of the column after any correction due to restraining moments have been = 
accounted for. A similar expression for M, can also be written, 


“ay 


-$M, + + M, = LY, 


in which Y, ; = the minor axis flexural stiffness calculated in a manner similar 


fects: 
The third equation of equilibrium relates to = ‘twisting ¢ effects a 


= | 


— (Cu + Cy)- —(Cy- Cr) (4) 

in which C, «= ae of the comp compression ‘flange alone; Cy = the si stiffness 
tension Senge alone; =F = the K, = the St. Venant’ 


ig q 
| 

a 

2.—Forces and Deflections Shown for Min 


COMPUTED ULTIMATE LOADS 
ry 


: = the term relating to the Wagner « effect and; D= ell 
a ‘The three foregoing | differential equations vary from point to point along 

the column length due to inelasticity. A numerical procedure (6) involving — 

successive over-relaxation has been used to solve the equations at a selected 
5 number of discrete points along the length of the column. 2 bl 

_ Failure Loads without Torsional Effects—The method is based on the assump-— 
4 tion that the orientation of the cross section along the length | of t lu 

remains s unchanged inspite of the biaxial deflections. The to 


solved 1 may be written as: 


of the setlted thrwst, measured with respect to the initially straight column 
axes, and corrected for any restraining end moments. 


- The net curvatures are compounded and the bideneation f the ‘neutral axis 
eed Assuming plane section remains plain, the strain distribution is 


established. From the stress-strain properties of the materials the internal stresses, — 


and thus moments, are calculated and compared with the external moments. 
The solution for the equilibrium deflected shape is obtained by using the 
Newton-Raphson iteration technique for a system of nonlinear equations (9). 
4 In both of the methods outlined previously, equilibrium « deflected shapes are 

a obtained for i increasing values of the loads following a given load path until, 
a for a given load, no convergence is obtained. Provided that the load increment — 
P is small, the highest load for which an equilibrium solution is found may be 

taken as the load for the column. 
Seco Pa 8 
cope OF Present Stupy 

parameters are likely to ‘influence twining warping in 


restraints. The calculated ultimate loads of torsionally stiff ¢ cross sections 
g not likely to be influenced by allowing for torsion and warping in the column. 
; Open sections like the I-section or H-section may, however, under certain 
conditions show instability to torsional and warping deflections. Thus, for the 
a present study an 8 x 8 x 48 Universal Column section was adopted. Being 
near the lower end of the range of practical sections, the cross section should — 
exhibit maximum slenderness effects even at normal story heights. Ree 
i In all, eight cases have been analyzed, with three different column lengths. — 
The three lengths of 90 in. (2,286 mm), 123 in. (3,124 mm), and 153 in. (3,886 
mm) correspond to weak axis slenderness ratios of 43, 59, and 74, respectively, — 
7 _ Fepresenting stocky, medium, and long columns. Six of me columns analyzed 
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were p pin- ended w while two two of ‘the longer columns had minor axis 
Table 1 lists the pertinent data. Fig. 3 shows various experimentally obtained 
stress-strain curves used in the analysis and referred to in Table 1. hi wt 
The columns analyzed are identified by labels, the first character representing, 
_ the length (S for 90 in., M for 123 in., and L for 153 in.) and the next three = 
digits representing the axial thrust on the column in ton forces. Restraint on _ 
the column is further identified by a trailing character (W for weak restraint, _ 
.and H for high restraint). Absence of a trailing character denotes no restraint. _ 


A 
Stress-Strain 
Curve 


(See Fig + 
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IG. —Strese- Strain Curves Used in Analysis 
about both the bending axes were considered. The initially 


~~ shape is assumed to be sinusoidal in nature, the maximum values 
being those given in Table 1. Small minor axis moments were applied, aimed 


he cee typically applied unintentional end moments. The columns wer 


| 
| 
Column | L,in | tons | uo, in | | force | | | 
inches | force | inches | inches | inches | B, B, | Flange Web 
$250 | 250 | 0.023 | 0.045 | 01s} 10] 2 | 4 
230 | 0.020 0.030 4.0 10.15] 1.0 
M200 123: | 200 | 0.015 | 0.070 0.15} 10} 2 |. 
L200 «| «153 | 200 | 0.000 | 0.077 | 40 jois} 10} 1 | 
LITO | 170 | 0.000 | 0.077 | 4.0 1.0} 1 
Lisow 150 | 0.015 | 0.110 | 00 |0.15 | 4 
| 
e 


analyzed for a major axis moment gradient of 0. 1s (ratio of smaller major 

g axis end moment to the larger value). This ratio is maintained throughout the cag 
path. The actual loading path adopted consisted of first applying the - 
o> axial load on the column, and then increasing the major axis moments __ 

rn TABLE 2. —Computed Ultimate Moments, in tons force- inches 
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Column label Torsion (Ref. 6) 9) 


635.06 
| 3 | 
= 


ltonf=9.96kN. | 
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4 
to The ‘ultimate load of the column the maximum value 
g the larger major axis end moment the column can support in equilibrium. ™. 
Resutts 3 3 
major axis ultimate moments obtained with and without 


torsional effects are ‘shown in Table 2. Results are also plotted in Figs. 411 
— in which the major | and minor axis deflections and twists are plotted for the 
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column satis as well as the major and minor axis end rotations. Continuous 
lines in these figures relate to the behavior when torsional effects are included, 


FIG. 5.—Rotations and Twist for Column $250 -in. = 0. KN. 


tonf-in. 


500 
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MOMENT 


° - 
ne. 6.—Deflections for or Columns M200 M230 a tonf- f-in. = 0.2531 0. 2531 | -m, in. 


; high load on the cross section. It may be seen from Figs. 4 and 5 that the > 
* response obtained by ignoring the torsional effects is essentially the same as 
~ that obtained by including them. Except for very small differences in minor 


. i and the broken lines relate to the behavior when twisting and warping are ignored. § 
7 
7 == 
* 
! 
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COMPUTED ULTIMATE LOADS 


axis deflections, the response is nearly identical. Cubiiagteliits, the twists 


recorded are also small. The difference in the two ultimate loads is only 10.2%, 
calculated as a of the load obtained by ignoring torsional 


TORSION 
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we 


i 
- deg. 


= 0.2531 kN- 
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4 


—Deflection for Columns L L140, 1170 ‘and L200 (1 in. = 0.2531 ‘we 
- effects. Strictly speaking, the more exact value based on torsion calculations — 
should have been used for calculating the percentage of errors. However, the | 
program including the torsional effects sometimes stopped a little 
short oft the —— value, whereas th the computer program i ignoring torsional ef! fects 


| 
ie 
g 
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7.—Rotations and Twist for Columns M200 
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, Figs. 6 and 7. For these columns, the behavior obtained on the basis of ceding 
torsional effects was somewhat different from that obtained by including them. 
; ‘ _ Divergence was noted in the minor axis deflections almost from the start, , but = ‘4 
le the major axis deflections were nearly identical. The ultimate loads obtained a 
were also very close, the difference being 2.8% for the M200 column, and 
_ 5.5% for the M230 column. A noteworthy feature was the substantial amount — 
Z of twist obtained for both onnae. Clearly, _ for such columns although the 


500 
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i: — L140, L170 and L200. —Three co columns umns of length 153 in. @B, 886 
mm) were analyzed. The axial loads were 140 tonf (1,395.0 kN), 170 tonf (1,693.9 _ 
kN), and 200 tonf (1,992.8 kN), respectively. The results for the three columns 
are shown in Figs. 8 and 9. The notable feature of the results obtained is 
that a marked divergence in behavior occurs when torsional effects are ignored, 
both in "terms of the ultimate moments and the failure modes. Inclusion of 

torsional. effects results in a reduction in the ultimate moments of 16. 9%, 23.8%, 
and 23.9%, respectively, for the three columns. Major axis deflections, however, — - 
show divergence at an early stage and, at collapse, the differences in the values x 

_ obtained are substantial. Perhaps the most significant aspect is the occurrence 

of very large angle of twist near A comparison with earlier cases 


wi 


was able to trace the response beyond the peak and for this reason the more 
| hused as the 
24 
length 123 in. 
(3,124 mm) were analyzed. One of the columns had an axial load of 200 —_ 
the may be subject toerror. > 
2 
| 
FIG. for Columns L140, L170, and = 0.2531 
» 
| 


ae Effect of Minor Axis Flexural End Res aints. :—The: — for all of t the columns 2% 
m previously, when including the effects of twisting and warping, 

ats, obtained assuming that the column was free to twist and warp at both end 


] columns are restrained : in SOME way both flexure 


GNORED 


DEFLECTION 


| 


(ROTATION - 


a. 


Pe In this section, two columns were analyzed with only minor axis flexural — 


2 restraints. Both columns were 153 in. (3,886 mm) long and carried an axial , * 
7 Risa) = 150 tonf (1,594.6 kN). One column (L150W) had a flexural end restraint = 


— 
— 


of stiffness 12, 000 tonf-i in. /rad 037.0 kN-m /rad), and the other (L150H), 
of 20,000 tonf-in./rad (5,061.7 kN-m/rad). Results are shown in Figs. 10 and © 
11. The most significant observation is that there appears to be a negligible 
(less than 1%) difference ‘in the ultimate moments obtained between the case 
when torsion is included, , and and the case when torsion is not included. | Major = 
axis deflections and rotations are virtually identical. It is only near collapse, 
that a slight divergence is obtained for the minor axis deflections. In spite — 
- of the closeness of the deflections obtained for the two cases, the values of < 
_midheight twists calculated were of significant magnitude. Thus, even though - 
; the twist does not appear to influence the failure load, it could govern the 
failure mode. It is important | to ‘Tealize ‘that in n the present c cases, _ only flexural 


Sev Furthermore, the beam stiffnesses of 12,000 tonf-in./rad (3,037. Oo 

kN-m/rad) and 20,000 tonf-in. /rad (5,061.7 kN-m/rad) adopted in the analysis 
were not of a very high order in relation to the column stiffness. The conclusion - 

= be that in the presence of flexural (or torsional, or both) restraints, a 

very good estimate of the behavior can be obtained | without considering the 

__ torsional effects of twisting and warping. In such cases, weak axis deflections 

_ are likely to be somewhat underestimated, but the estimate of failure loads — 

— will be fairly accurate. In view of the fact that the inclusion of twisting and 
warping effects imposes a severe penalty on computational effort in terms of — 

= computer time © (average | ratio A 5), this result can have an important bearing = 


| 


he 


occurred. o~ minor axis moments were of the order of unintentional end 
- moments, usually specified in codes of practice. It is considered that the load 


path chosen is likely to show maximum interaction between the flexural and 


— Ithas been shown that in many cases for steel I-columns or H-columns subjected 
to biaxial bending, a good estimate of the ultimate load can be obtained without a 
considering twisting and warping of the cross section. The conditions under 7 
which this applies are: (1) Short and medium length columns with or without — 
flexural and torsional end restraints; and (2) long columns with ee and an 

‘For long columns with the minor axis slenderness ratio greater than int 7 
* 60, the failure mode is likely to be accompanied with predominant twisting 


and warping. For such columns, neglecting these effects could lead to an erroneous © _ 


estimate of the load deflection response. Also for such columns, the ultimate — 
loads calculated could be in nonconservative error by as much as as 


4 is, a need for further systematic study of this problem. 


\g The work described in this paper was carried out in the Department o of Civil — 


‘research leading to their doctoral dissertations. 0. 
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_The folowing symbols ar 


"stiffness of the compression and tension flanges respectively; 


- distance between the center lines of flanges; 8 
~ biaxial eccentricities of effective line of action | of "applied thrust" 
with respect to initially straight column axes; 
= shear modulus; 
torsional stiffness relating to Wagner ‘effect; 
St. Venant’s torsion constant; onstant; 


M,,,M » = "moments at bottom of column; 
moments at top of column; 


M,,M, major and minor axis moments; 


biaxial deflections; 


= majorand minor axis flexural stiffnesses calculated d by inch incteding 
effects of material 4 


e used in this 
age 
“ef 
| 
| 
G 
s 


x pap Rite coordinate axes for cross es in its s initially ‘straight he 


= distance along column height; 
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yet, to the point where they warrant publication as a Proceedings paper in a Journal, the a 
Baba: of Technical Notes was extort by the Board of Direction on October 16-18, 
1967, under the following guidelines 


. An culajaat manuscript and two copies are to be submitted to the Manager | of Technical — 
and Professional Publications, ASCE, 345 East 47th Street, New York, N.Y., 


10017, , along 
with a request by the author that it be considered as a Technical Note. 


2. The two copies will be sent to an appropriate Technical Division or Council for review. — 
3. If the Division or Council approves the contribution for publication, it shall be returned ; 


to Society Headquarters with appropriate comments, 
4, The technical publications staff will prepare the material for use in the earliest possible 


sue of the Journal, after proper coordination with the author, 2% 
‘a 5. Each Technical Note is not to exceed 4 pages in the Journal. As | an approximation, 


full manuscript page of text, tables, or figures is the equivalent of one- eee a —_— 


The Technical Notes will be grouped in a special section of each Journal Sr 
. Information retrieval abstracts and key words will be unnecessary for Technical Notes. 


The final date on which a Discussion should reach the Society is given as a footnote 
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FLExuRAL CRACK WIDTH IN PRESTRESSED 
hes (ha fore ConcrETE BE AMS haw 
arotd By Stephen W. Meier,’ A. M. ASCE and Peter Gergely,” M. 
Partially prestressed concrete beams, those which are designed f for a limited — 
devel of tensile stress, have | several advantages over their fully prestressed 
counterparts—there is a better control of camber, a Savings in prestressing” 
. a : steel and end anchorages, and a more economical utilization of mild steel. Partially 
Prestressed concrete beams do, however, have several disadvantages: (1) Reduc- 

; tion in the ultimate strength; (2) increased deflection under overload conditions; Jf 
and (3) cracking. This note is based on a recent study (9) of flexural cracking © 
ia partially prestressed concrete beams which included th the examination of all 

Prediction of maximum crack widths for prestressed concrete beams has 
received far less attention than that for reinforced concrete. Experimental data 
related to prestressed concrete cracking are extremely limited and often collected — 
: as secondary information, and yet the complexity of the crack width calculations © 
are greatly increased over reinforced concrete by the presence of additio 1 


_ The width of cracks in a concrete beam « depends on the strain in the member, ° 


of the member. The majority of previous investigations approach crack width — 
prediction using either the 1 net steel stress method 1 3, 5,6, 10, 11) or the nominal 


on the net change in the steel stress from the decompression moment 
sponding to zero strain on the tension face) to the load stage at which the 
crack width is to be determined. The nominal concrete method is based on — 
the fictitious tensile stress in the concrete at the extreme fiber, assuming an 
uncracked, homogeneous section, 

onl The limited test t data available have been treated empirically (9) by examining» 

trial prediction equations composed of measured and calculated variables. 
‘Design Engr., Chicago Bridge and Iron Co., Chicago, 

Re. ?Prof. of Struct. Engrg., Cornell Univ., Hollister Hall, Ithaca, N.Y. 14853. eLokienas 
_ Note.—Discussion open until July 1, 1981. To extend the closing date one month, 

\ sd a written request must be filed with the Manager of Technical and Professional Publications, oe 
ASCE. Manuscript was submitted for review for possible publication on March 20, 1980. 
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Geometric and physical variables may accuracy 
while load- -dependent variables are to uncertainties of prestressing 


‘the: location of the neutral a axis were an iterative 

equilibrium solution of a modified reinforced concrete beam similar to that 
by Bennett and Chandrasekhar (4). 

A total of 48 trial prediction equations were selected to compare the sensitivity — 
— applicability of the two different approaches, to examine the variables, 5 

a... to study new and previously | proposed peng The evaluation of coeffi- z 


Data were Mk from seven different sources (1,2,3,5,6,8,11) containing _ 


‘rectangular, I- or T-sections with deformed bars, prestressing strands or wires, a 
and varying geometric and physical properties. The data \ data were divided into tw« two 


—Number of Deta Points in DataSet 


Complete ; Limited 
(2) 


Modified strands 


major categories: deformed reinforcing bars and strands, depending on the type = 

_ of steel closest to the tensile face, because of their differing bond characteristics. 

_ Additional refinements to eliminate bias were reauired since the Nawy investiga- : 
- tion (11) had a significantly greater number of points than other studies and 
_ it represented beams with similar geometric and physical characteristics (the = 
data are identified as “modified strands’). A number of data points were 
eliminated because the crack widths were less than 0.003 in. (0.076 mm); such 
small values are very erratic and therefore cannot be predicted with acceptable _ 7 
accuracy. Data points corresponding to stresses greater than 2,500-psi (17-MPa) _ 
nominal tensile stress were also removed because such high Stresses are 
_ impractical. The remaining data set is identified as “‘limited’’ in Table 1. 
yr Histograms of the number of data points were constructed for each of the 
: major variables for each data set to indicate their range and possible bias present 

in a data group. A reliable statistical analysis must examine the ranges of all 

possible factors or variables in the data because a variable cannot be included _ 

in a prediction equation if it is not properly reflected in the data. Likewise, 

if a variable assumes only a few values in many beams, its significance may 

be erroneously reflected in a statistical evaluation. Some of the equations proposed \ 


- (10) correlate well with the data from which they were sone but are poor 


in predicting crack widths from other 18) 


q 
| 
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The trial prediction equations were chosen to investigate the effect of the 3 
_ major variables and to compare previously recommended equations. Five equa- 
“tions were variations of the Bennett and Chandrasekhar "equation (3), which | 
has the form w = kef .,, in 1 which = a constant; c= the net cover; and 
Sa = the nominal tensile stress at the tensile face. A « comparison of these : 
_ equations supports the findings of other investigations that the crack width 
is proportional to the clear cover. The coefficients k obtained (9) were 2.5 
Fr for rebars and 3.7 for strands, which are 83% and 74%, respectively, of those a 
The e addition of the effective area p per bar in the equation shows little 
_ improvement i in the prediction for the data set. The use of the cube or square 
root of the effective area per bar does, however, seem to have an advantageous | " 
effect on the reliability of the prediction. This supports the findings of the § 
_ Gergely-Lutz equation (7) for reinforced concrete - When a an equation v with the ca 


aa ‘the scatter is somewhat larger ‘than in several other equations because = 

of the inclusion of V d. instead of d. which was found to be more applicable. Oo 

lair 


20 


FIG. 4 of Ratio of Measure 

A comparison of equations based ¢ on the ne net - steel stress approach with those | 
containing the nominal concrete stress indicates a significant increase in scatter 


In light of the extensive statistical analyses, the following equations ns are proposed 
4 fort the prediction of of ma: maximum flexural crack width: few 


in 5 which w = maximum crack width; C,, C, = bond coefficients which he 
on the type steel nearest to the tension face; for rebars C, = 12, Cy 
for strands C, = 16, C,= 12 23 = nominal ‘concrete tensile ‘strain a at tl 
tensile face, equal to f,,, in which E. = the modulus of elasticity of concrete; 
d. = minimum concrete cover to centroid of steel nearest the tensile face; 
= effective concrete area per t bar as defined in the American Conties 


Eq. 1 


. 3 but Eq. 2 shows somewhat better accuracy om: wide beams ¥ with large horizontal 


| 
q 


spacing. The use of the tensile ‘strain €., . makes Eg. 1 dimensionally correct. ne 
equations poetics the of the maximum crack widths. The scatter 

1 which shows the histogram of the 
ratio of measured and peter values for prestressed beams containing pris 
deformed bars. The prediction is not as good for beams with strands and without | 
_ deformed bars. The accuracy of the other 46 equations was found to be inferior : 
to that of Eqs. 1 and 2, even though some were 


‘The maximum crack width in concrete beams has been shown to be related — 
to the strain and curvature in the member. The relationship between the crack 

— width and the Strain in the member is usually modeled by the stress or strain 

“in the steel, | or by the nominal tensile stress in the concrete, modified by a 
multiplier to’ account for the bond and physical characteristics of the section. 

_ Examination of cracking behavior relative to previously established factors does 

- not justify complex prediction expressions based on the limited data available. _ 

_ Prediction equations based on the nominal concrete stress has the advantage 

_ of avoiding the iterative calculations required to obtain the steel stress. Compari- a 

son of the results obtained here for the two approaches: suggests that the net 7 

steel stress method has little advantage over the simpler nominal stress method. — 
Proper design and detailing practices can minimize the cracking problems — 

in partially prestressed and reinforced concrete members. Distributing the steel — 

and placing deformed bars, which have superior bond characteristics, near the — 

— tensile face can distribute the cracks and reduce the maximum crack width. __ 


Further experimental investigations, covering an expanded range of the major 
variables with a larger data base, are “necessary before significant advances _ 
can be made in flexural crack width prediction for partially prestressed concrete 

beams. The two equations proposed should suffice until more data becomes 
available for the estimation of flexural crack width. 
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MULTISPAN ConTINuOoUS BEAM 
COEFFICIENTS 


Gerard C. ‘Pardoen,' M. AS ASCE 


The objective of this modeling procedure is to isolate a single span of a 
-multispan continuous by beam while providing for the rotational and translational = 
rigidities: of the Spans flanking this’ single span. Consider, for example, the 
mth span of the ‘multispan continuous beam depicted in Fig. 1. In determining 


“te rotational and translational rigidities of the spans flanking this mth span, 


the following assumptions are made: 

© The flanking spans are Supported by rigid rollers within the interior of 

the continuous spans and | are ‘supported at the end by ‘either a roller [as in 

Fig. 2(a)] or a clamped support [as in Fig. 2(b)). 1... ; 

2. Each individual span is prismatic and can be characterized by a different — 

length, L,, and a constant bending rigidity, E/,. The effect of shear deformation 

lcmmding the usual translational and rotational degrees-of-freedom at each | 

end of a beam, the 4 x 4 stiffness matrix for the unsupported mth span can 


‘Asst. Prof., Civ. Engrg., Univ. of California at Irvine, Irvine, Calif. 92717. ae 
_ Note.—Discussion open until July 1, 1981. To extend the closing date one month, 
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= flanking. spans, sespectively. The stiffoess cosfficients he mth span me Kw 


associated with the flanking spans depend upon several factors, among t them 


Continuous Bea 


being the number of individual span spans, s, the bending rigidity and length of each’ 7 
__ individual span, and the support condition at the end of the flanking span. . 
_ These coefficients are derived by sequentially imposing unit translations and — 
rotations at each end of the mth span and denoting the reaction forces needed ~ 
to maintain the flanking span in the appropriate displaced 


Eno Rower Support on Eno CLampeo SuPPorT 

-“ Consider, e.g., n individual spans flanking to the right of ‘span m in Fig. 

2, in which the individual spans are supported by rigid rollers within the interior => 


of the continuous span and at the end by a roller. Assume each 
span has a different length, L, (i al Ae 2. .. n), but also has, at least for 


symmetrical 0 0 
4 
an 


the time being, a constant bending rigidity, EI 
can be derived from matrix techniques and shown to be 


"Expressions for the effective length t terms 4 sppeating in 2 are given 
in Table 1 for one to five individual flanking spans" with the same ‘bending 
| + 4 


q Translational Stiffness Coefficient Kaa 

Keo 


“rigidity, EI,. If, however, any the has a constant bending 
rigidity other than EI, over the span length, then one needs only to normalize | 
all ill bending rigidities t to El,. Specifically, this can n be accomplished most | easily — 
by by r replacing all terms L, appearing in the definitions of L* and L* with L, Et, 
es practical purposes, one would rarely need to account for more than five 7 
_ individual flanking spans because the computed values of K,,, K,,, and Kg, o 
are relatively unaffected by the presence of even the fourth and fifth individual 
* spans. This conclusion i is based upon the curves shown i in Fig. 3(a) fora multispan: : 
continuous beam of ‘equal span length | and bending | rigidity, EI. These curves, . 
. which depict the influence of the number of individual flanking s Spans in computi 
K4,, K 49, and Ky, reflect the tabular valuesin Table2, 
ye Similarly, the stiffness coefficients K,,, K 49, Ko for the n flanking. 
a camped ent end 700 can be shown to of of 


ficients 
4 
im 
‘ 
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TABLE 1.—Stifiness Coefficients of Flanking ‘Spans with Roller End; Ks, 2 
(BEI, + 4L3)/GLF + 4L$); = + 4L3)/GL? + 
and K,, = (3EI,/L, + 4L3)/GL3+4L3) 


| flanking 


spans | 


IL, L, L L,(G/4)L,L,L, + G/4)L, al 
L2+L,L + (9/16) L,L2+ (3/4)L,L,L, 


5 
L,+(3/4)L2L, + B/4)L3L, + (3/4) L2L 


} 


c 
= 
w 
= 
a 
a 
c 
£ 
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Stiffness Coefficients - Clamped Support at End aie 


cee 


= 
4 
{ 
| 


TECHNICAL 


L, 


Expressions fort the effective terms * Ls ‘appearing in Eq. . 3 are given 


in — 3 for one to five flanking spans with the same bending rigidity, El, n 


oe ved of 
TABLE 2.—Tabular Value es of Stiffn 


2.286 | 2.385 | 2.392 | 2. .392 | 2.392 | 2.392 | 2.392 | 2.392 
K L?/3EI 1.429 | 1.462 | 1.464 | 1. 
1.143 | 1.154] 1.155 | 1. 


- TABLE 3.—Stiftness Coefficients of Flanking Spans with: Clamped End; K,, = 
/ Li + LE) + 413); K = GEL, + 2L3)/(4L} + 4L3); and 


= (GET, /L, + +403) th | uw & - 


flanking 


neat 


+(G/4)L, L2+(3/4)L,L, Le + (3/4)L, L? i+ G/4)L,L 


+ (3/4) L?2L, + 9/16) L? ibs. +L, 
+ (3/4)L,L,L, 1/L? 
_ Again, if any of the flanking spans has a constant — rigidity other than 
EI, over the span length, then one only need the terms appearing 
_ in the definitions of and L} with L,/E EI,. 
a Again, there is little need to account for more than five flanking spans because _ 
‘the computed values of K,,, K,,, and K,, are relatively unaffected by the 
presence of even the fourth and fifth flanking spans. This conclusion is based — 


ess Coefficients Reflected by Curves in Fig. 3(a) 
7 
| 
| 
i 
q 
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2 
bo? 
5 
upon the Curves shown in Fig. [or a muluspan Continuous OF equal 


‘span — and bending vin ee These curves, which depict the influence 


of the number of flanking spa spans ir in computing the 
tabular values in Table 4. 


Eun 

As an example of the ts consider the. continuous 

beam system shown in Fig. 4. Each span is of length, L, and has bending 

rigidity, EJ, except for span, AB, which has bending rigidity 2EJ. Assuming 

TABLE 4 _—Tabular Values of Stiffness Coefficients Reflected by Curves in Fig. — 


0.598 | 0.598 | 0.598 | 0.598 | 0.598 | 0. 598 7 
1. 000 0. 750 0. 733 0. 732 | 0.732 | 0.732 | 0.732 | 0.732 | 0.732 | 0. 732 A 

1.000 | 0.875 | 0.867 | 0.866 | 0.866 | 0.866 | 0.866 | 0.866 | 0.866 | 0.866 o 
‘ rae 


Span of Interest - Spon 


Left Flanking Seen - Span AC 


that span CD is the span of significant interest, | then span AC is eomsideced 


the left flanking span, and span DG is the right flanking span. The overall 


4x 4 stiffness matrix associated with span CD (before the implementation 
of the boundary conditions) i 


— 
0 -12 -6L 12 -6L 


0 0 93 STIL 


r 
| 
j 
| 


The first 4 x 4 matrix widhidin the rigidity of the left flanking span “AC oan 
is derived from Table 3. Note that the special consideration is made for the 

- gigidiny of span AB. The second 4 x 4 matrix reflects the rigidity of span : 
CB. The third 4 x 4 matrix reflects the rigidity of the right flanking - : 


“a Implementing the displacement | unine conditions at rollers C and D gives 
the matrix equilibrium equation for s CD as 


sults of Eq. 5, 6.=- 
agree with those when o one e considers the span as a Six degree- 


= 
4 The sine and formulation for re representing the translational and rotational 7 
‘Tigidities of flanking spans of a multispan, continuous beam have been presented. : 
_ This modeling procedure, which is particularly helpful if the significant eatin 
occurs within a single span, can drastically reduce the kinematic indeterminacy 
of an analysis problem. Based upon several examples, there appears to be no — 
significant effect upon the translational and rotational rigidities when the numbers f 
of flanking spans exceeds three, compared to the rigidities calculated upon f 
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Discussions 
— may be submitted on any Proceedings paper or technical note published i in 2 1 any 


or on any paper presented a at any Specialty Conference or other meeting, the | Proceedings 


anyone who has significant comments or questions regarding the content of the paper/technical © 
note. Discussions are accepted for a period of 4 months following the date of publication 
a a paper/technical note and they should be sent to the Manager of Technical and Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 
extended by a written request from adiscusser. 
7 The original and three copies of the Discussion should be e submitted on 8-1/2-in. (220-mm) 
11-in. (280-mm) bond Paper, typed double- -spaced with wide margins. The length of 
® manuscript including figures and tables); the editors will delete matter extraneous to ie 
subject under discussion. If a Discussion is over two pages long it will be returned for shortening. 
_ All Discussions will be reviewed by the editors and the Division’s or Council’s Publications — 
- Committees. In some cases, Discussions will be returned to discussers for rewriting, or they — 
may be encouraged to submit a paper or technical note rather than a Discussion. _ ae 
Standards for Discussions are the same as those for Proceedings Papers. A Discussion is 
subject to rejection if it contains matter readily found elsewhere, advocates special interests, : 
is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, 
or is foreign to the purposes of the Society. All Discussions should be written in the third 
person, and the discusser should use the term “‘the writer’’ when referring to himself. The 
author of the original paper/technical note is referred to as “the author. 

Discussions have a specific format. The title of the original paper/technical note appears — 
a the top of the first page with a superscript that corresponds to a footnote indicating the 
month, year, author(s), and number of the original paper/technical note. The discusser’s full 
“ame should be indicated below the title — Discussions herein as an example) together with 


§ page of the manuscript, "along with the Proceedings paper number of the original eoneg ve Deno 
- note, the date and name of the Journal in which it appeared, and the original author’s name. ~ we 
- Note that the discusser’s identification footnote should follow consecutively from the original _ 
_ paper/technical note. If the paper/technical note under discussion contained footnote numbers 
7 | and 2, the first Discussion would begin with footnote 3, and subsequent Discussions would _ 
continue insequence, 
_ Figures supplied by the discusser should be designated by letters, starting with A. This leo 
_ applies separately to tables and references. In referring to a figure, table, or reference that : 
appeared in the original paper/technical note use the same number used in the original. 
It is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 
the inse-esiguoa of ASCE for more detailed information on preparation and submission of — 
manuscripts. 
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irtep Desicn or CorruGaTeD SHEAR 


Closure | by Kinh H. Ha,° M. ASCE, Noor 'El-Hakim,’ sradw 
od) and Ps Paul P. M. ASCE aqile mase 


q _ The writers wish to thank Easly and Chockalingam for their discussions. — js 
- In Ref. 14, Easly has used the term “‘seam slip’’ A, in the definition of the 

seam fastener force F, (Eq. 14, Ref. 14), and the relationship between | the 
was established by the testing of “bolted lap joints.”’ This fact led the 
writers to believe that Easly’s “seam: slip’’ represents the total deformation 
1 in the seam connection. However, as clarified in the discussion, the term “‘seam 

_ Slip’’ was used to denote only one half of the deformation in the seam connection. _ 
_ Since the deformation in the corner end connection which is in line with the | 
‘ seam equals A,, the force in that connection is only half of that force in the 
connection (assuming all fasteners to be identical). This" fact, which 
contradicts Fig. 5 of Ref. 14, has been “extensively v verified by finite element. 
analysis and was also recognized by Davies (9). Thus, Easly’s clarification of 
the term ‘‘seam slip’’ still does not correct the deficiency of his theory. 
The writers Fh with Easly’s contention that the perimeter framing member _ 
should pass through the centers of all the In fact, under such 


a condition, one must have (Fig. 3): she we 


Fig. may give” rise to possible violation of only ‘the moment equilibriam 
; - the interior panels. The reason for the writers’ emphasis on ensuing equilibrium 
of the end panels has already been stated inthe paper, 
a Eq. 37 may also be interpreted as the condition for all panels to deform — 
identically; thus, the writers find it difficult to support Easly’s contention that 
his theory i is applicable to the case in which all fasteners are identical. ip de 
The data in Table 3, prepared by Chockalingam, are interesting; however, 
ee writers feel that still better agreement between the predicted diaphragms’ 
: strengths and measured values would be obtained had actual fasteners’ strengths 
been available. The fact that these diaphragms failed along the seams is further _ 
phrag 
“- “July, 1979, by Kinh H. Ha, Noor El-Hakim, and Paul P. Fazio (Proc. Paper 14681). wl 
 °Assoc. Prof., Centre for Building Studies, Concordia Univ., 1455 de Maisonneuve ~ 
Doctoral Candidate, Centre for Building Studies, Concordia U niv., 1455 de Maisonneuve 
7 Bivd. W., Montréal, Québec, Canada, presently, Sr. Structural Analyst, Pratt and Whitney — 1 
Canada Ltd., Longueuil, Québec. a Univ., 1455 de Pager 
*Prof., Dir., Centre for Building Studies, Concordia Univ., | vd. 
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‘evidence of for t the seam fasteners’ forces. by the writers’ 
f method. Thus, the writers are surprised to see the poor agreement in the predicted | 
and measured seam slip values in Table 4. However, careful eee of these © 
values and t their sources (2) re reveals the points 
1. The seam slips predicted by the writers’ theory | are. the a average | values, 
whereas those by finite element analysis and by experiments are at the — : 
centers where they re 
2. The seam slips in Col. 2 of Table 4 are tcdbniotent with the fastener 
forces given in Ref. 2. By using the fastener stiffness and forces, at the gage 
location, as obtained by finite element analysis (2), the seam slips in diaphragms 
57-2 and 59-4 should be 0.0048 in. and 0.00415 in., respectively. It is worthwhile 
to note that Cols. 2 and 5 of Table 1 show excellent agreement between the 
. writers’ method and finite element analysis for the average seam forces along a 


already been ‘explained by Ammar and Nilson 

~ Concerning Chockalingam’s | s final comment, ~~ eee would like to to point 
out that no distribution of the horizontal end fastener forces within a panel 
is necessary nor implied in the theory, and thus there is no inconsistency in 
this respect while calculating D, and D,. . Finally, i in light of the foregoing facts, 

the writers feel that the “discussers have been unduly overconcerned with a > 
possible violation of the moment equilibrium of the interior panels. waste ryt 
_ Errata.—The following corrections should be made to the original paper: a 


Page 1366, Should send “Bryan (3), and Bryan and El-Dakhakni (4, 5 6). 
Davies subsequently modified some of Bryan’s expressions by assuming an 
improved distribution of the fastener forces (8,9). In general, Bryan’s theory | a 
_ is based on the energy method and’’ instead of ‘‘Brian (3), and Brian “— : 
7 El-Dakhakni (4,5,6). Davies subsequently modified some of Brian’s } expressions © 
: by assuming an improved distribution of the fastener forces (8, 9). | In In general, 
Brian’s theory is based on the « energy method and” 
= Page 1367, paragraph 1, line 1: Should read “Similar to Bryan’s method,’ instead 
Page 1370, paragraph 2, line 2: Should read ‘‘the same as those derived by 
Bryan” instead of ‘‘the same as those derived by Brian” 
Page 1373, Table 2: Fo 2: For r diaphragms 57-2 57-2 and 59 59-4, data in Cols. 4 and 6 and 6 should 7 


| 

| 
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_ The authors are to be congratulated for bringing to notice some of the spurious © 
‘ie cases that can occur i ‘in numerica al models for inelastic - material behavior. The 
purpose of this discussion is to supplement | the original technical note. The 
thors have correctly said that Clough’ s model (2) has probably _ inspired | 
numerous other models. The o- B model (6) developed by the writer is one. 
eae model, most effectively shown in Fig. l(a), gives the gist, and is 
not to be directly extended into details suck, as shown i in Fig. 1(6). Its extension — ; 
oF Fig. 2 2(6) might be valid. The writer’ s material m model (6) was specifically 4 
3 delat to simulate the special stress-strain behavior of mild steel in monotonic 
and reversed loadings. As a particular case, Clough’ s model is obtained. A 
a more general stiffness degradation model without strain hardening was also’ 


developed adopting the Masing technique, 


_ The writer’s paper (6) explains some more spurious cases that can happen 
4 in | material modeling an the | consequences of avoiding « or not avoiding these 
spurious behaviors. For example, Fig. 2(b) explains two situations, one when > 
a 7 unloading occurs up to +0 load level and then reverses, the second when unloading 
occurs up to —0 load level and then reverses. The writer has discussed a similar 
situation. As the authors would have meant, it is not beyond the reach of 


the analyst- programmer to avoid such irregularities. But at what mental effort 
and physical cost vis-a- a-vis the accuracy achieved? 
- To answer this, an important distinction has to be made as to ‘the meth¢ dology 


of solution built-in into the software. A nonlinear analysis starting from he 
- modeling of the basic stress-strain behaviors of the constituent materials is 
relatively expensive, but the overall effects of ignoring the spurious cases of 
type | Fig. _2(6) may not be serious. the other hand, it may 


e against the stress-strain relationship) as, ¢.g., moment-curvature, axial 


force-axial deformation etc. This is evident from Fig. 3(a) and “a .. > 


_ The writer wishes to mention that the authors’ new stiffness degrading model, 


‘Fig. 4, is closely obtained as 7 special case of the writer’s a-B mild steel stress 


strain model (6) by using the stiffness ; degradation Parameter a =la 


‘Santhanam, es “Model for Mild Steel in Inelastic Frame Analysis,”’ Journal 
the Structural Division, ASCE, Vol. 105, No. ST1, Proc. a od 14333, Jan., 1979, 7 


3 “December, 1979, by Rafael Riddell and Nathan M. Newmark (Proc. Paper 15017). _ 


*Scientist, Structural Engrg. Research (Regional) CSIR Madras- 
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STANDARD PRESTRESSED Concrete Be: BEAM im 


Discussion by Bohdan K. Boczkaj’ 


It is a duty of an author to give precise definition of new idea or a concept 
intends to introduce. As a benefit, can understand Paper easier 
s author has additional ‘opportunity to find mistakes 
the case of papers by D. S. Hatcher (6,7) including the o one under discussion. 
J Introduction of the ‘‘Loading Interaction Diagram’’ and we Stress Interaction — 
Diagrams,” one for “‘stresses’’ at transfer, second for “‘stresses”’ at service 
load, without accurate explanation of coordinate systems leads to ambiguous — 
results in relations between all of these diagrams. In the original paper, the 
word stress was, in the writer’s opinion, improperly used in reference toa 
_ force or a ‘moment. For the sake of better understanding the writer decided , 
to keep it in this discussion but in quotation marks. The name “Stress Interaction 
7 7 Diagram,”’ itself, for a diagram in coordinates, P-force and M-moment does 
not help one to understand the method of superimposing it with “Loading | 
Interaction 1 Diagram. ” The “Stress Interaction | Diagram” for transfer gives a 2 
region of prestressing force, P , and bending moment M for which the stresses 
in any point of section are in range (a,f’ or © f'). Similarly, at service load a 
the range is a ‘f!, ce’ f). Both types of diagrams (Figs. 2, 3, 4, 5) are drawn 
in the same coordinate system. It makes the impression that these two diagrams» a 
overlap each other by wide areas and at our disposal is a wide range of prestressing _ 
force and moment, it means also eccentricity. It is easy to 
two diagrams are in two different coordinate systems: P,, , M, = + 
for transfer, and »P,, M,, = —nP,e + M, + M, for service load, and ovedapeing, 
area is very small, as one may see it on Magnel’s diagram and on Fig. 8. 
In the case of a fully stressed section, it is only one point. It is easy to ages 
it, when both diagrams are drawn in one coordinate system, bad Pa as i: 
For the fully stressed section, Eqs. 1 and 3, and also Eos. 2 
identical and only one pair of P and e satisfies Eqs. 1-4, so any constraint 
on eccentricity makes that section nonfeasible. It also applies to the range of 
eccentricity. The chapters called ‘Sections Having Maximum Eccentricity 
Constraints’’ and ‘‘Sections Having Range in Eccentricity’’ in Ref. 6 for this 
aa lead to false results. The method of calculation k, + k,—Eq. 39 in 
Ref. 6—is not always applicable, because it leads, especially for bigger spans, 


 Indescribing ‘ ‘Stress Interaction Diagrams”’ ‘instead of describing certain points, = 


§ seems to the writer, that it is better to describe. ‘the | diagram as an area 
inside lines as given by the following equations: — SSS 


“January, 1980, by David S. Hatcher (Proc. ‘Paper 1510 15106). 


Engr., Dravo Corp., Pa. 


» 8. —Stress Interaction in Diagrams w with Area 


transfer a, 


: 


za af. of to le. 
nd 


meh fhe 

in which = e+M M, 

The shape of fully stressed section on coefficient A = (nc, — a 

— na,). When A > 1, the fully stressed section has a convoidel @ axis in the 

upper half, as on Fig. _ 3(c) in Ref. 6. When A = | section is symmetrical, 

as on Fig. 3(a) in Ref. 6. When A < 1, the section has a centroidal axis in 
= the lower half, as on Fig. 3(b) in Ref. 6. 

Tt is necessary to add to chapter ‘‘Review”’ that Magnel created a ‘graphical 

, method to find eccentricity e and prestressing force P, (or their reciprocal 1 /P, 7 - 

for the section, which was earlier selected and which had tos satisf. y two vo equations: 


pe ake When Me of 
| 
| 


 Egs. 57 and 58 58 are are derived fi trad aiih's 1 a 3 and Eqs. 2 and 4, 4, respectively. 

‘When selected section satisfy both these equations, the line corresponding to 7 
. 6 lies below that corresponding to » Eq. 8 and the line corresponding to 


Eq. 5 lies to right of that corresponding to Eq. 7; there is no need to check 


selection of the eoction , P, and e may be calculated using Magnel’s diegrem 
or by solving one combination (of a possible four) of two equation chosen 
a Egs. 1-4. Eqs. 1 and 4 give P,,,;, and ee The combination of Eqs. p 
2 gives P, =  P, (Eq. 24). The combination of Eqs. 3 and 4 gives 
Py (Eq. 25). The combination of Egs. 2 and 3 gives Fiat and e,,;,- It is of 


- 


| 


P.nax > P, and P,,,, > Pp-, but the section is feasible, as 


éman (max 


€ = Cin: bined Prestressing force and eccentricity has to be in the range P,, 


max’ Cm 


tmax 


_A design procedure outlined in “Noncomposite 


in "Fig. represent A and M for the ‘combination ‘of ‘equations indicated by 
_ Coefficients R, and R, are described as greater than 1 when af. #0 - 
a’f’ #0. It is a mistake, because Me R, > 1 only when a,f’ and a’f’ 

0. When and > 0, R,, Ry < 1, when and = 0, R,, R 


introduction of the concept of generalized kern (14), but it also can be an. 
asa way the prestressing force, or compression resultant moves inside generalized - 
kern after a load is applied. The total distance a resultant can travel without | 
overcoming allowable stresses is R,k, + R rr as shown in Fig. 9. After 
_ ~prestressing, ' when the weight of the beam M,, Starts to act, a resultant, or 


her equations ac linng ac the eccentricity ic within an allowahle range A pal 1 
{ 
g 


compression force is shove eccentricity e on M,/P, to the 

owe limit of the generalized kern. Allowable stresses are not exceeded and _ 
the stress diagram has the biggest value at the bottom. When the value of _ 
prestressing force drops to P,, the resultant moves up for a distance (1 — > 

1) M,/nP,. The last part of the denominator of Eq. 31 in Ref. 6 represents 


this distance, which has to be provided to keep the resultant all the time within © 
the limits of the kern R,k, + R,k,. When M, is applied, the resultant moves 
to the upper limit of the kern ‘ond the stress diagram has the biggest value 4d 


@ Stresses on top. Using this jaa ws pg Eq. 21 in Ref. 6 can be easily 


* The method Sudesh in the paper does not apply to pretensioned sections, 
because in this case in calculation of prestressing force moment M = M, +. a 
Mg has to be taken, instead of M = M, + (1 — 4)Mg, as it is in post-tensioned 
‘Structures. Also eccentricity for a pretensioned section can be fast shed e=R re 
of the end of the beam. . = = 


(14. Kaufman, S. , Mosty sprezone (Prestressed Bridges), ‘bed Poland, 6 Polish). 


EE 


‘The author has procedure for the preliminary 
_ estimation of reinforcement tatio for reinforced concrete columns subjected 

biaxial bending ‘moments. According to the author, the method could be 

* used to ) predict ‘the area of reinforcement in an assumed cross s section, after — 
which a detailed analysis using the available methods could be made. But the © 
method seems to lack a theoretical basis or a proof (as to why Eqs. 8a and 
8b, which are derived from uniaxial bending analysis would give valid rein- 
forcement ratios in a biaxial bending case), and does not indicate the errors — 
inherent in the derivation of the proposed e equations. T The writer wishes to to make 4 


correct form of Eq. 4 is: 


“May, 1980, by Himat T. Solanki (Proc. Paper 15384). 


| 
| 
| 
| 
rgh, Pa. 15261. 


in which all ter terms are as defined in in the paper. 


The author did not indicate as to how Eq. 5 was affected by the ae, 
that the compression steel is located very close to the middle of the compression 
Eq. 9 cannot be derived directly from Eqs. 6 or 8a, as the depth to the 
_ center of the compression block for the example presented in the paper is 
1. 68 in., which is smaller than a minimum effective cover of ai in. caeee 
to this, Eq. 6 has’ to be modified as following: 
A, -—- —h+ gh) 
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e, 
assuming an effective cover of 2 in., the alg the compression steel is” 


Ag) 


| FIG. 5.—Biaxial Bending 


e+e (in which and é, = eccentricities along Y and 1X directions, 
respectively). By drawing an arc from the origin O, with radius e, it could 
be seen from Fig. 5 that when the load is acting at point A, e, = e and e, 
7 0; and ey when the load is p acting at point B, e,=e and e, = 0. 


> ACL and compression s steel strain < 


e — 0.5¢ + 0.5 


Su. 
| 
explained =—s 
through Fig. 5 in which X and Y axes coincide with the major and minor _ 
2 
— 


egils 


The author has. ab thas. the 


a _ other terms are as defined in the paper. 


the is also assumed here that A in or 


e-—-O.5at+ Vv: 


at t 0.85 

in whete 2A, the area of reinforcement on the sides parallel to the a ie 
Eq. 14 is identical to Eq. 8a and Eq. 15 could be made identical to Se 
¥ 8b by multiplying with the factor (1 + a)/(1 + 2a/3). This alternate method 
for deriving the equations: proposed by the author (Eqs. 8a and 8b), shows 
that the e philosophy of the approach i is to ) assume the net moment ) pr a to act 
in both directions X and Y. The writer tends to believe that the major saat 
of this method is that for any combination of eccentricity e, and e,, which 
_ makes the net eccentricity e to fall along the arc AB, Eqs. 8a and 8b predict 
the same reinforcement ratio. If, for the — presented in the ‘paper, M 


er: ft _kips and M, = 357 ft kips or, M . = 200 ft kips and M, = 300 at 


ies ; the re resulting net eccentricity is the same as that given in the paper. ‘Therefore, : 


. for both these cases, the maximum reinforcement ratio given by either Eqs. 
_ 8a or 8b, would be the same as that determined in the paper. This does not 
seem to be correct because, in the former case, the design is nearly identical 
to that of the uniaxial case, whereas in the latter, more reinforcement is needed 
along the longer side tov withstand a comparatively larger moment. The equations — 


7 _ proposed by the author do not ‘seem to take into account of this aspect and 
thus the writer tends to doubt the general usefulness of the approach. peiryy Dh) 
_ In addition, for the moment (P,e) acting in the X or Y directions, one could 
determine the reinforcement ratios (for an assumed cross section) directly from 
the uniaxial diagram without solving Eqs. 8a and 8b or one could 


_ It can also be verified from Eq. 76 that for a = ae 3 ‘ oxeate section), p, = 
1.2 (2A,/A,), and that given by Eq. Ta is = 2A,/A,. Therefore for a square 
cross section with equal moments acting in both directions, — Eqs. 7a and 76 


“predict tw two different reinforcement ratios. 


from Eqs. 8a and 8b depends on the assumed cross- s-sectional dimensions. The 

proposed equations could be valid only for a certain range of moment combinations 

and a. But the major problem in the design of reinforced concrete column 

_ subjected to biaxial bending moments, lies in finding a suitable cross- sectional - 
dimensions and the proposed me method does not address the issue. — bi 
oe 


pillidedre _ gens Wisew « bebi vor 


an Das syllivoq Jad! bejon isan boiw 


*‘Columns,”’ Design Handbook, Vol. 2, American Concrete Institute, Detroit, Mich., ca 


8. Somayaji, S., and Ramakrishnan, V., of ‘Simplified Method of Design of Reinforced 


Concrete Columns Using Equilibrium Diagrams,” ” Journal of Civil | Engineering Design, 
ie 


| 
| 
| 
| 


The attractive title to the original the reader to 
a simplification of biaxial bending. The ‘‘simplification’’ does not look simple. 
_ Procedures recommended elsewhere (3,4,6,7,8) require fewer calculations, less — 
engineer time, and they produce more accurate results. The proposed simplified 


selecting a cross section as if the resultant skew moment were acting about 
the principal axis of the cross ‘section. Eqs. 8a and 8b were derived on the 
basis of assumptions that are nearly appropriate only if large eccentricities exist — 
_ for bending about both principal axes. Even if the equations produced good © 
estimates of steel required ‘‘in opposite faces,’’ no advice is given regarding 
the dual role of | the bars i in corners of sections. dns dod 
— ‘The difficulty encountered by y efforts to apply the proposed method is revealed 
clearly in the example illustrative problem. The proposed procedure led to a 
recommendation for a 4% reinforcement ratio. A cross section with only 2.7% 
_ reinforcement (10 #10 bars, four in each side face and three in each end face) 
can be shown adequate for the specified loading. If the ‘‘simplified”’ Procedure — 
leads to half again as much ste¢ steel as that which is needed, the simplified SS 
needs refinement or rejection. In | its present forms it does iis look promising 
APPENDIX. —Rererences 


Furlong, R. W., ““Concrete Columns Under Biaxially | Eccentric Thrust,” Journal, 


10. Furlong, Richard W., “ Ultimate Sista of ‘Square Columns U Under Biaxially Eccentric 
PDFo OF Wind: Press ON Low-Rise But Butpincs" 


The author hi has provided a a useful analysis ¢ of ‘the 


ell, Jr. 


1980, by Theodore Stathopovlos (Proc. Paper 15409). 


fr. Lect., ept. IV. ngrg ‘ost ice, James niv., OwnsvVvl1 
— 2Sr Lect., Dept. of Ci &S Engrg., Post Office, James Cook Univ., T. ill 


Texas, at Austin, Austin, 78712. 


7 procedure is not simpler nor is it as accurate as the elementary expedient of _ 
| 
| 
using wind tunnel measurements. It is noted that positive skewness and non- — 2 : 
| = 
. 


hin Sete tails were initia for the windward wall pressure fluctuations, in 
og addition to the well- established negative skewness and non-Gaussian tails asso. 

ciated with negative mean pressures. The author has also observed that the 

positive skewness cannot be entirely explained by the skewness of the upwind > 

It should be pointed out that there is a nonlinear relationship between upwind — 4 

_ velocity fluctuations and windward wall pressure fluctuations, apart from the 

_ distortion: of the turbulence due to the flow around the building. The writer 

Pe has recently (18) investigated the probability distribution of dynamic pressure — | 
fluctuations associated with Gaussian velocity fluctuations. It is shown that 7 
the PDF of corresponding pressure coef (t ), is is given byt the following 

4 


the tails, when the turbulence intensity, is around 20%. 

- ~The latter is about the average turbulence intensity at the building height for 
_ the “‘open country”’ and “‘built-up”’ exposures used by the author. This distribution 
showed excellent egretintat with the measured distribution for the pressure | 


fluctuations on the windward wall of a house measured by Ho Holmes and Best | 


+ 


| ea Although the departures from the Gaussian form at the positive tails, in the 
E author’ S measurements, are somewhat greater than those indicated by Eq. 21. 

& for J, = 0.20, it seems likely that the nonlinear pressure-velocity relationship — 
the major “‘other factor’ contributing to the skewness" distribution 


18. J. D., “Non- Gaussian Characteristics of Wind Pressure Factions Journal — 
of Wind Engineering and Industrial Aerodynamics, Vol. 7, No. 1, Jan., 1981, pp. 
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